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A REAL-TIME DIGITAL COMPUTER PROGRAM FOR THE 
SIMULATION OF A SINGLE -ROTOR HELICOPTER 

By Jacob A. Houck, Lucille H. Gibson,* 
and George G. Steinmetz 
Langley Research Center 

SUMMARY 

A computer program (Langley program Cl 152) has been developed for the study of 
a single-rotor helicopter on the Langley Research Center real-time digital simulation 
system. This report includes descriptions of the helicopter equations and data, program 
subroutines (including flow charts and listings), real-time simulation system routines, 
and program operation. Program usage is illustrated by standard check cases and a 
representative flight case. 


INTRODUCTION 

In view of the expanding interest in helicopters at the Langley Research Center, the 
development of a man-in-the-loop real-time helicopter simulation was necessary. An 
existing computer program for a single- rotor helicopter was modified and implemented 
on the Langley Research Center real-time simulation (RTS) system. This program was 
then used to gain experience in the area of nelicopter mathematical modeling (specifically 
in the area of rotor mathematical modeling), and it was used to determine the impact 
(computer memory and timing requirements) of such a simulation program on the comput- 
ing facility at the Langley Research Center. 

The objective of this paper is to present the development of the computer program 
for the RTS system. Equations, aerodynamic data, flow charts, computer listings, and 
operating procedures for the simulation program are given. 

The simulation program, with modifications, is presently being used in support of 
two piloted studies of the Sikorsky S-61 helicopter, which is a five-blade, single-rotor, 
commercial passenger-type helicopter currently being used by New York Airways. One 
study concerns pilot workload, route structures, and air traffic control procedures for 
instrument flight in congested short-haul markets. The second study concerns develop- 
ment and evaluation of various computer- generated displays and instruments for aiding 
pilots during flight, especially during approach and landing. 

Electronic Associates, Inc., Hampton, Va. 



SYMBOLS 


Measurements, calculations, and programing were made in the U.S. Customary 
Units. They are presented herein in the International System of Units (SI) with the 
equivalent values in the U.S. Customary Units given parenthetically. 

Most of the symbols used herein are not standard but were adapted from the 
FORTRAN notation used in the computer program of reference 1. 

A SPD indicated airspeed, kilometers/hour (knots) 

A OS main rotor collective pitch, radians 

A OS£ collective pitch control input, radians 

A 1S main rotor lateral cyclic pitch, radians 

A 1S C lateral cyclic pitch control input, radians 

AlS CS SAS lateral c y clic P itch control input, radians 

ADVPj aerodynamic variable parameter, per meter 2 (per foot 2 ) 

ADVP 2 aerodynamic variable parameter, per second 

a X ,a Y’ a Z H ne ar acceleration along X, Y, and Z body axis, respectively, 
meters/second 2 (feet/second 2 ) 

B lg main rotor longitudinal cyclic pitch, radians 

®1Sq longitudinal cyclic pitch control input-, radians 

B 1S CS SAS longitudinal cyclic pitch control input, radians 

C J SAS constants used in transfer function 1, where i,j * 1,2,3 

c D^ y main rotor drag coefficient at ^ R and Y 

C L ^ y main rotor lift coefficient at and Y 
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Q 


'Xi 


w 


X , 


C Yi 


C Y 


C Z 


C TR 

VS 

F 


J HS 


J W 


'm. 


' m Fr 


n t 


n T 


D 




D 


W 


power coefficient 
torque coefficient 
thrust coefficient 
wing drag coefficient 

fuselage drag coefficient as a function of sideslip 

fuselage drag coefficient as a function of angle of attack 

fuselage side-force coefficient 

tail rotor side-force coefficient 

vertical- stabilizer side-force coefficient 

fuselage vertical-force coefficient 

horizontal -stabilizer vertical-force coefficient 

wing vertical-force coefficient 

fuselage rolling- moment coefficient 

fuselage roll-rate damping coefficient, kilogram -meters (slug-feet) 
fuselage pitching- moment coefficient 

fuselage pitch- rate damping coefficient, kilogram -meters (slug-feet) 
fuselage yawing-moment coefficient 

fuselage yaw-rate damping coefficient, kilogram-meters (slug-feet) 
summing variable (main rotor blade drag), meters Vsecond^ (feet^/second^) 
rotor downwash coefficient 
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d x longitudinal distance from center of gravity to main rotor shaft (positive 

rearward of shaft), meters (feet) 

d XrTo longitudinal distance from horizontal stabilizer to main rotor shaft (positive 

Ho 

rearward of shaft), meters (feet) 

dx_ D longitudinal distance from tail rotor to main rotor shaft (positive rearward 

T xv 

of shaft), meters (feet) 

dx,, c longitudinal distance from vertical stabilizer to main rotor shaft (positive 

Vo 

rearward of shaft), meters (feet) 

dx w longitudinal distance from wing to main rotor shaft (positive rearward of 

shaft), meters (feet) 

d z vertical distance from center of gravity to main rotor hub (positive upward 

to hub), meters (feet) 

d Z TR vertical distance from tail rotor to main rotor hub (positive upward to hub), 
meters (feet) 

d z vertical distance from vertical stabilizer to main rotor hub (positive upward 

V b 

to hub), meters (feet) 

d z ^ vertical distance from wing to main rotor hub (positive upward to hub), 

meters (feet) 

E Earth coordinate in east direction, meters (feet) 

E velocity component in east direction, meters/second (feet/second) 

e flapping hinge offset, meters (feet) 

Fmrs main rotor hub moment parameter, kilogram-meters2 (slug-feet^) 

Fxr main rotor force along X body axis, newtons (pounds) 

F Y r main rotor force along Y body axis, newtons (pounds) 
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G eh main rotor ground effect due to height 

G E y main rotor ground effect due to forward velocity 

G RW t vehicle gross weight, kilograms (pounds) 

g acceleration due to gravity, meter s/second^ (feet/second^) 

h Earth vertical coordinate, meters (feet) 

h velocity component in vertical direction, meters/second (feet/second) 

hjy desired velocity component in vertical direction, meters/second 

(feet/second) 

Ig blade moment of inertia, kilogram-meters2 (slug-feet^) 

Igg horizontal- stabilizer built-in incidence angle, radians 

Ig rotor moment of inertia inboard of flapping hinge, kilogram- meter s2 

(slug-feet2) 

Irqt main rotor moment of inertia, kilogram-meters2 (slug-feet^) 

Iyg vertical-stabilizer built-in incidence angle, radians 

I^y wing built-in incidence angle, radians 


*XX’*YY>*ZZ moment of inertia about X, Y, and Z body axis, respectively, 

kilogram- meter s2 (slug-feet^) 


SAS constants used in transfer function 2, where i = 1,2,3 and j = 1,2, 3, 4 



main rotor constants for flap angle calculation 


area constants used to determine main rotor lift and drag forces, 
meters2 (feet^) 
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Ki 


Ei 


constant used in engine torque calculation, kilogram-meters^/second 
(slug-feet^/ second) 


K M >Kq>Kql volume constants used to determine main rotor moment and torque, 

meters^ (feet^) 


K R/F 

K r '»2 


Ki 


W 


IWM 


J MR 


J PDR, 


'I'Y 


J RH 


Mt 


M 


RH 


M 


TIP 


m 


m 


rotor downwash fraction acting on fuselage and wing 
ratio of elevator area to horizontal-stabilizer area 
constant used in calculation of Yrpj^, newtons (pounds) 
wing downwash coefficient 

total aerodynamic rolling moment, newton-meters (pound-feet) 
fuselage rolling moment, newton-meters (pound-feet) 
main rotor lift force, newtons (pounds) 

summing variable (main rotor blade lift), meters4/second2 (feet^/second^) 
main rotor hub rolling moment, newton- meters (pound-feet) 


total aerodynamic pitching moment, newton-meters (pound-feet) 


fuselage pitching moment, newton-meters (pound-feet) 


MpoR^y summing variable (main rotor blade flapping moment), metersVsecond^ 
(feetVsecond2) 


main rotor hub pitching moment, newton-meters (pound-feet) 


main rotor tip Mach number 


B 


vehicle mass, kilograms (slugs) 


blade mass, kilograms (slugs) 
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N 

N 


AZ 


B 


N 
N 

N F 

n rad 

P 

P 

Q e 

q e d 

Q 

^raax 


Q 

Q 

q 

q 

q T 


MR 


pdr* y 


Rt 


r 


Earth coordinate in north direction, meters (feet) 
velocity component in north direction, meters/second (feet/second) 
total aerodynamic yawing moment, newton- meters (pound-feet) 
number of rotor azimuth stations 
number of rotor blades 

fuselage yawing moment, newton-meters (pound-feet) 

number of blade radial stations 

roll rate, radians/second 

roll acceleration, radians/second2 

engine torque, newton- meters (pound-feet) 

desired engine torque, newton- meters (pound-feet) 

maximum engine torque, newton- meters (pound-feet) 

main rotor torque, newton-meters (pound-feet) 

summing variable (main rotor blade torque), metersVsecond^ (feet5/second2) 

pitch rate, radians/second 

pitch acceleration, radians/second2 

lateral dynamic pressure, newtons/meter2 (pounds/foof2) 
vertical dynamic pressure, newtons/meter 2 (pounds/foot^) 
main rotor blade radius, meters (feet) 
yaw rate, radians/second 
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yaw acceleration, radians/second^ 
shaft horsepower, watts (horsepower) 
fuselage lateral-force area, meters2 (feet^) 
fuselage vertical-force area, meters2 (feet^) 
wing longitudinal-force area, meters2 (feet^) 
fuselage lateral-force area, meters2 (feet^) 
vertical-stabilizer force area, meters2 (feet^) 
fuselage vertical-force area, meters2 (feef2) 
horizontal-stabilizer force area, meters^ (feet^) 
wing vertical-force area, meters2 (feef2) 
time, seconds 


blade element velocity perpendicular to blade span axis and to U T at 
and Y, meters/second (feet/second) 

blade element velocity perpendicular to blade span axis and to rotor rotation 
axis at '1 >j^ and Y, meters/second (feet/second) 

linear velocity along X body axis, meters/second (feet/second) 

linear acceleration along X body axis, meters/second2 (feet/second2) 

speed of sound, meter s/second (feet/second) 

total velocity of air through main rotor, meters/second (feet/second) 
tail velocity along Y body axis, meters/second (feet/second) 


fuselage rolling- moment volume parameter, meters^ (feet^) 
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v m fuselage pitching-moment volume parameter, meters^ (feet^) 

F 

V n fuselage yawing- moment volume parameter, meters3 (feet 2 ) 

n F 

v linear velocity along Y body axis, meters/second (feet/second) 

v linear acceleration along Y body axis, meters/second 2 (feet/second 2 ) 

W F main rotor inflow distribution at any radial station 

W Hg tail velocity along Z body axis, meters/second (feet/second) 

W T i 7 main rotor inflow redistribution as a function of velocity, meters/second 

ir 

(feet/second) 


W IWM 

W %Y 

w 

w 


main rotor mean inflow velocity, meters/second (feet/second) 

wing downwash velocity, meter s/second (feet/second) 

main rotor local inflow velocity meters/second (feet/second) 

linear velocity along Z body axis, meters/second (feet/second) 

linear acceleration along Z body axis, meters/second 2 (feet/second 2 ) 


total aerodynamic force along X body axis, newtons (pounds) 

X AOS cockpit collective stick input (positive up), centimeters (inches) 

X cg cockpit cyclic stick longitudinal input (positive forward), centimeters (inches) 

X F fuselage aerodynamic force along X body axis, newtons (pounds) 

X tr cockpit tail rotor pedal input (positive right), centimeters (inches) 

X w wing aerodynamic force along X body axis, newtons (pounds) 

Y main rotor radial distance (measured from hinge), meters (feet) 
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Y a total aerodynamic force along Y body axis, newtons (pounds) 

Y cs cockpit cyclic stick lateral input (positive right), centimeters (inches) 

Y F fuselage aerodynamic force along Y body axis, newtons (pounds) 

Yp E main rotor radial station plus hinge offset, meters (feet) 

Ytr tail rotor force along Y body axis, newtons (pounds) 

y TW blade twist at any radial station, radians 

Yyg vertical-stabilizer force along Y body axis, newtons (pounds) 

Z A total aerodynamic force along Z body axis, newtons (pounds) 

Zp fuselage aerodynamic force along Z body axis, newtons (pounds) 

Zrs horizontal-stabilizer force along Z body axis, newtons (pounds) 

Z w wing aerodynamic force along Z body axis, newtons (pounds) 

ofp fuselage angle of attack, radians or degrees 

ofjjg horizontal- stabilizer angle of attack, radians 

wing angle of attack, radians 
a egg main rotor coning angle, radians 

Q! jgg main rotor longitudinal flap angle, radians 

d-lgg main rotor longitudinal flap rate, radians/second 

main rotor blade angle of attack at <k R and Y, radians 
|3p fuselage angle of sideslip, radians or degrees 

0yg vertical- stabilizer angle of sideslip, radians 
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P 1SS 

^1SS 

At 

6 E 

9 

0 

9 TR 

9 


TRc 


0 


TR CS 


'J'Y 


P 

P 


$ 

<t> 

't' 


main rotor lateral flap angle, radians 

main rotor lateral flap rate, radians/second 

main rotor total blade flap angle at radians 

main rotor total blade flap rate at 'l'^, radians/second 

time interval, seconds 

horizontal- stabilizer deflection angle, radians or degrees 

fuselage pitch angle, radians 

rate of change of pitch angle, radians/second 

tail rotor pitch angle, radians 

tail rotor pitch control input, radians 

SAS tail rotor pitch control input, radians 

main rotor blade pitch angle at ^ ^ and Y, radians 

main rotor tip- speed ratio 

air density, kilogram/meter 3 (slug/feet^) 

time lag, seconds 

fuselage roll angle, radians 

rate of change of roll angle, radians/second 

main rotor blade inflow angle at and Y, radians 

fuselage yaw angle, radians 

rate of change of yaw angle, radians/second 


'I'p main rotor blade azimuth angle, radians 

S2 main rotor angular velocity, radians/second 

tl main rotor angular acceleration, radians/second^ 

Op desired main rotor angular velocity, radians/second 

Dot over symbol denotes time derivative. 

PROBLEM DESCRIPTION 

A description is presented of the development of an all digital simulation program 
for a single-rotor helicopter to be run on the Langley Research Center real-time simula- 
tion (RTS) system. The program was designed to represent the entire operational flight 
envelope from hover through transition to forward flight. 

The equations of motion for a single- rotor helicopter were developed in reference 2. 
They are suitable for analog, digital, and hybrid simulation programs. The equations 
include dynamic modeling of the main rotor and airframe to provide a mathematical model 
which will accept control inputs and calculate the resulting vehicle dynamics and 
orientation. 

In reference 1, the helicopter equations of motion were applied to a specific vehicle. 
The simulation consisted of an all digital airframe mathematical model programed on the 
Electronic Associates, Inc. 8400 digital computer and a stability augmentation system 
(SAS) programed on an analog computer. The digital programing language was 
FORTRAN IV. 

References 1 and 2 have been used to form the basis of the Langley helicopter sim- 
ulation. A complete collection of aircraft equations is contained in appendix A. Several 
modifications and updatings of the computer program have been instituted and are dis- 
cussed herein. The simulation program has also been structured to run on the Langley 
RTS system which is described in reference 3. Operating procedures are discussed in 
appendix A of reference 4, and in the Program Organization and Program Usage sections 
of this paper. 

One of the major modifications to the computer program involved a change in inte- 
gration formulas. Reference 1 uses an Euler formula to integrate the state and rotor 
variables with a solution rate of 25 iterations per second. The Langley program uses a 
second-order Adams- Bashforth predictor integration formula with a solution rate of 
32 iterations per second to solve the state variables. This allows for a more accurate 
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solution and is one of several integration formulas available through the RTS system soft- 
ware. This integration of the state variables occurs only in the operate mode of the sim- 
ulation. There are five sets of rotor variables which are integrated continually, regard- 
less of the simulation computer mode, in order to stabilize the rotor dynamics. They 
are the inflow velocity, the main rotor coning angle, the two components of the main rotor 
flapping angle, and the main rotor angular velocity. These variables are integrated by 
using an Euler formula at a solution rate of 32 iterations per second. 

Other modifications involve function generation and trim condition calculation. The 
function generation system used in reference 1 has been replaced by a technique used for 
several years with the Langley RTS system. This technique allows for rapid generation 
of aerodynamic functions used in real-time simulations. The technique is detailed in 
reference 5 and in appendix B by Lawrence E. Barker, Jr., and Kemper S. Kibler. The 
function data are presented in appendix C . The trim program of reference 1 has been 
replaced because of its slow convergence rate and its poor accuracy. The new trim cir- 
cuit, a modified secant method, is detailed in reference 6 and has been reproduced for 
convenience in appendix D by Gary A. McDaniel. It has been used in several RTS system 
programs and is rapidly becoming one of the standard techniques used with this system. 
The technique is extremely fast and requires only a few seconds to reach trim criteria at 
speeds up to 140 knots in level flight. A forward velocity of 140 knots appears to be the 
upper limit for trim convergence in level flight in this program due to the uncertainty in 
the helicopter aerodynamic data. Trim conditions at higher velocities have been reached 
in descents; and with an improved data package, it is felt that trim conditions at higher 
forward velocities can be reached in level flight. 

A stability augmentation system (SAS) has been implemented in the Langley simu- 
lation. This system is not discussed except to mention that it was necessary to develop 
special integration techniques due to the numerical instability of the standard RTS inte- 
gration formulas when handling time constants of the magnitude involved. These convo- 
lution integration techniques are derived from reference 7. 

One final feature has been added to the simulation program. This feature is the 
ability to conduct preprogramed flight maneuvers from the program control station (fig. 1) 
and allows the analyst to run research and validation tests on the mathematical model 
without having a pilot present to fly the simulation. The operating procedure for this 
feature is detailed in the Program Usage section of this paper. 

With minor changes to the mathematical model (data, control systems equations, 
fuselage equations, etc.), the program described in this paper can be used to simulate a 
variety of single-rotor helicopters. Future studies include air traffic control, control 
systems, navigation systems, and pilot display systems. 
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PROGRAM ORGANIZATION 


Langley program Cl 152 was written in the FORTRAN IV language to be run in a 
real-time mode on the CDC 6600 computer system. The program requires a field length 
of 65 000 octal locations. The computing time depends primarily on the integration 
scheme used. For the normally used one-pass Adams -Bashforth integration scheme, 
the computing time required is approximately 11.25 milliseconds (total available is 
31.25 milliseconds) per iteration. 

This section contains a list of the FORTRAN variables in labeled COMMON and the 
various DSPLAY arrays, a brief description of the real-time system subroutines, and 
individual program subroutine descriptions with corresponding flow charts and listings. 
The following diagram shows the deck configuration needed for execution: 

ft I 

7 

8 

9 


f= 

Function Data 


V 


8 

9 


C1152 Main and 
Subprograms 
Block Data 



/ 


^ 


Real-Time 
Control Cards 

/ 


Deck Configuration 
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Labeled COMMON 


The following list contains the FORTRAN variables appearing in labeled COMMON 


and descriptions of each variable: 


COMMON 

label 

FORTRAN 

variable 

Description 

/CHECK/ 

TABLE (150) 

Array of floating point numbers to be displayed 
and/or changed 


INTEG(09) 

Array of integer numbers to be displayed and/or 
changed 

/INTCOMM/ 

T 

Time, update in subroutine IGRATE1 


DT 

Integration step size used in subroutine IGRATE1 


INT 

Flow control parameter used in subroutine 
IGRATE1 


NEQ 

Number of integrations performed in subroutine 
IGRATE1 


ISC HE ME 

Selects integration scheme in subroutine IGRATE1 


DERINT(2,12) 

Array of integrals and derivatives in subroutine 
IGRATE1 

/INTINT R/ 

INTERN(5,12) 

Temporary storage for elements of subroutine 
IGRATE1 

/REALTIM/ 

ADC(32) 

Array of analog-to-digital converter inputs 


DAC(64) 

Array of digital-to-analog converter outputs 


LDISI(108) 

Array of logical discrete inputs 


LDISO(196) 

Array of logical discrete outputs 
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COMMON 

label 

FORTRAN 

variable 

Description 

/REALTIM/ 

NOPER 



NHOLD 



NRESET 

> Return addresses from subroutine RTMODE 


NTERM 



NPRINT 



NREAD 

J 


/ACCEL/ 

AX 

a X 


AY 

a Y 


AZ 

a Z 

/ADV C/ 

OMEG 

ft 


WIM 

W IM 


VT 

v T 


QV 



QL 



ASPD 

a spd 

/ADV P/ 

ADVP1 

ADVPj 


ADVP2 

advp 2 
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COMMON 

label 

FORTRAN 

variable 

Description 

/ANGSAS/ 

PHIL 



THETAL 



PSIL 

Last values of Euler angles and their derivatives 
' used in SAS calculation 


PHIDL 


THETDL 



PSIDL 


/COEF/ 

CT 

C fp 


CP 

Cp 


CQ 

C Q 


MU 



MTIP 

Mtip 


VS 

v s 

/CONTRL/ 

AOS 

A os 


A1S 

A 1S 


BIS 

B 1S 


THTR 

0 TR 


OMEGD 
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COMMON 

label 


/CNTL C/ 


/CNTL P/ 


FORTRAN 

variable 


Description 


AOSC 

A1SC 

B1SC 

THTRC 

A1SCL 

B1SCL > 

THTRCL 

AOSCO 

A1SC0 

B1SC0 

THTRCO 

J 

XAOSG ^ 

YCSG 

XCSG 

XTRG 

XAOSR 

YCSR 

XCSR 

XTRR j 


A os c 

A 1S C 

B 1S C 

0 TRC 


Storage locations for last values of A ls , B 1S 
i C C 


and 6 , 


THc 


Control constants in A 


OS c > A 1S C > B 1S C 


and 0 - 


TRc 


Control gains 
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COMMON 

label 

FORTRAN 

variable 

Description 

/CNTSAS/ 

A1SCS 


Als cs 


B1SCS 


Bls cs 


THTRCS 


0TR cs 

/CPIC/ 

ACI 




BCI 

> 

Instrument lag parameters 


DC I 




XRCM1 




YRCM1 




XRTM1 




YRTM1 




XPBM1 

> 

Instrument constants 


YPBM1 




CSK4 




CSK6 




CSK7 




F05(18) 



/ENG C/ 

QMAX 


Q 

^max 


QE 


^E 


SHP 


S HP 


COMMON 

label 

FORTRAN 

variable 

Description 

/ENG C/ 

GRWT 

g rwt 


OMEGDT 

n 

/ENG P/ 

IROT 

X ROT 


CKE 1 

K El 

/ERROR/ 

NERR 

Error return parameter 

/EUL CS/ 

COSFI 

cos($) 


SINFI 

sin(<i>) 


COSTH 

cos{0) 


SINTH 

sln{9) 


COSSI 

cos('l') 


SINSI 

sin('I') 

/FCNNAME/ 

CLF1 

C ^F 1 = f (M 


CMF1 

c m Fi = f (« F ) 


CNF1 

Cn F 1 =f (M 


CXI 

C X 1 = f (%) 


CX2 

C X 2 =f(«F) 


CYTR 

C Y tr = f ( u - V TR^TR) 


CYF 

c y f = f (%) 


CZF 

c z F = f ( a f) 
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COMMON 

label 

FORTRAN 

variable 

Description 

/FCNNAME/ 

CXW 

c x w - f ("w) 


CZW 

c z w ■ f (“w) 


CYVS 

C Y V S ' ! (Ns) 


CZHS 

C Z HS " f (“Hs) 


CDSIY 



CLSIY 

= t (“' I ' Y ' UT «) 


DW 

D\V = f( u ) 


GEV 

Gev = 


GEH 

g eh = 

/FNMC/ 

XA 

X A 


YA 

y a 


ZA 

Z A 


LA 

l a 


MA 

m a 


NA 

n a 

/F NMP/ 

DX 

d x 


DXHS 

%s 


DXTR 

d X T R 


DXVS 

dx vs 


COMMON 

label 


FORTRAN 

variable 


Description 


/FNMP/ 


/FTRMPLS/ 


/FUNCS/ 


DZ 

d Z 


DZTR 

d Z TR 


DZVS 

03 

> 

N 

T3 


DXW 

dv 

A W 


DZW 

“ z w 


IR 

*R 


EQLMR "1 



1 

EQAOSS 

) Trim parameters 


EQWIM v 



HDOTD 



F001(28) 

Array of data values for 

c m F 

F002(28) 

Array of data values for 

Cv 

a 2 

F003(14) 

Array of data values for 

C 7 

^F 

F004(14) 

Array of data values for 

c / 

* F 1 

F005(14) 

Array of data values for 

C n 

n Ff 

F006(28) 

Array of data values for 

C Xj 

F007(14) 

Array of data values for 

c y f 

F008(28) 

Array of data values for 

Cx w 

F009(14) 

Array of data values for 

o 

N 
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COMMON 

label 

FORTRAN 

variable 

Description 

/FUNCS/ 

F010(28) 

Array of data values for C v 

Y VS 


F011(35) 

Array of data values for Cg 

HS 


F013(07) 

Array of data values for Dyy 


F014(07) 

Array of data values for G^y 


F015(07) 

Array of data values for Ggpj 


F016(35,10) 

Array of data values for C n 


F017(35,10) 

Array of data values for Ct „ 


F018(7,3,4) 

Array of data values for Cy,^ 


D001(9)-D011(9) 

Arrays containing the generated function val- 
ues and function parameters for F001-F011 


D013(9)-D015(9) 

Arrays containing the generated function values 
and function parameters for F013-F015 


DO 16(18) -D017(18) 

Arrays containing the generated function values 
and function parameters for F016-F017 


D018(27) 

Array containing the generated function value 
and function parameters for F018 

/FUS C/ 

XF 

Xp 


YF 

y f 


ZF 

Z F 


LF 

l f 


MF 
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COMMON 

label 

FORTRAN 

variable 

Description 

/FUS C/ 

NF 

n f 


XW 

x w 


WIWM 

W IWM 


ZW 

Z W 


ALFF 

“F 


BETF 

^F 


ALFW 

“w 

/FUS P/ 

CLF2 

C; 

1 F 

2 


CMF2 

o 

3 

CO 


CNF2 

Cn F 

F 2 


CKRF 

k r/f 


CKWIWM 

Kw IWM 


IW 

% 


VLF 

% 


VMF 

s 

> 


VNF 

< 


SXF1 

Sv 

X Fi 


SXF2 

Sx F 

F 2 


SYF 

Sy 

y f 
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COMMON 

label 

FORTRAN 

variable 

Description 

/FUS P/ 

SZF 

s z F 


SXW 

Sx w 


SZW 

£ 

N 

GO 

/GPARAM/ 

DELT 

At 


RHO 

P 


G 

g 


PIE 

7T 

/LAGSTK/ 

NILAG 

Time lag input parameter 


NSTICK 

Console potentiometer input parameter 

/MODEC/ 

MRESET 

Mode control reset 


MHOLD 

Mode control hold 


MOPER 

Mode control operate 

/MROT C/ 

FXR 

f xr 


FYR 

f yr 


LMR 

l mr 


LRH 

l rh 


MRH 

m rh 


QMR 

q mr 


COMMON 

label 

FORTRAN 

variable 

Description 

/MROT C/ 

AOSS 

a OSS 


ALFSIY 

“*Y 


UTSIY 

U T 

X *Y 


BETSI 



BDTSI 



A1SS 

a iss 


AD1SS 

a iss 


BISS 

^1SS 


BD1SS 

^1SS 

/RCOEFF/ 

CKL(20) 

k l 


CKD(20) 

k d 


CKQ(20) 

«Q 


CKDL(20) 

k dl 


CKQL(20) 

“ql 


CKM(20) 

k m 


Y(20) 

Y 


WF(20) 

W F 

/ROTCON/ 

YPE(20) 

y pe 


YTW(20) 

y tw 
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COMMON 

label 

FORTRAN 

variable 

Description 

/ROTCON/ 

WIF(20) 

Wjf 


BOVN 

Ratio of Ng to N ^2 

/ROT P/ 

TWIST 

Blade twist 


FMRS 

f mrs 


CKAO 

k a 

a O 


CKAB1 

k ab x 


CKAB2 

k ab 2 


CKAB3 

m 

00 


NRAD 

n rad 


NAZ 

n az 


NB 

n b 


RB 



EFH 

e 

/SAS P/ 

A1C1 

c ! 


A1C2 

C 2 


A1K1 

A 


A1K2 

K 2 


A1K3 

K l 


A1K4 

K l 
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COMMON 

label 


/SAS P/ 


/SCALES/ 


FORTRAN 

variable 

Description 

B1C 1 

C 2 

L 1 

B1C2 

C 2 

B1K1 

K 1 

B1K2 

K 2 

K 2 

B1K3 

A 

B1K4 

A 

THC 1 

c 2 

THC2 

c l 

THK1 

A 

THK2 

K 3 

z 

THK3 

A 

THK4 

A 


SF1 
SF2 
SF3 
SF4 } 
SF5 
SF6 
SF7 j 


Cockpit scale factors 
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COMMON 

label 

FORTRAN 

variable 

Description 

/SCPSIR/ 

SINSIR(50) 

Sin pR) 


COSSIR(50) 

cos (*r) 

/SHIP P/ 

IXX 

: xx 


IYY 

l YY 


IZZ 

X ZZ 


MASS 

m 

/STICKS/ 

XAOS 

X AOS 


YCS 

Y cs 


XCS 

x cs 


XTR 

X TR 


DPR 

Degrees per radian conversion factor 


RPD 

Radians per degree conversion factor 

/STINPUT/ 

XI 

Initial input on X^Qg 


X2 

Initial input on Y^g 


X3 

Initial input on X^g 


X4 

Initial input on Xrp^ 


T1 

Time at which to apply stick input 


T2 

Time at which to remove stick input 
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COMMON 

label 

FORTRAN 

variable 

Description 

/TAIL C/ 

ALFHS 

“hs 


BTVS 

Vs 


ZHS 

Z HS 


YVS 

Y VS 


VTR 

V TR 


WHS 

W HS 


YTR 

y tr 


DELE 

6 e 

/TAIL P/ 

CKTR1 

^Rj 


CKTR2 

KTR2 


IHS 

J HS 


IVS 

Vs 


SZHS 

Sy 

Z HS 


SYVS 

Sy 

*VS 

/TRIMAT/ 

XMAT(11,12) 

An (NEQN x NEQP1) matrix of independent trim 



variables 


EMAT(12,13) 

An (NEQP1 x NEQP2) matrix of functional 



evaluations 


TRSUM(12) 

An NEQP1 dimensional vector of error sums 


NEQN 

Number of equations to be trimmed 
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COMMON 

label 

FORTRAN 

variable 

Description 

/TRIMAT/ 

NEQP1 


NEQN + 1 


NEQP2 


NEQN + 2 


XMATG(ll) 

Increments added to XTRIM to start the iteration 

/VALUES/ 

A1CS1 


SAS value of first transfer function on A lg 


A1CS2 


SAS value of second transfer function on A lg 


B1CS1 


SAS value of first transfer function on B lg 


B1CS2 


SAS value of second transfer function on B lg 


THCS1 


SAS value of first transfer function on 0 TR 


THCS2 


SAS value of second transfer function on 0 TR 

/XIC/ 

A1X1 ^ 




A1X2 




A 1X3 




B1X1 




B1X2 

> 

Second transfer function parameters 


B1X3 




THX1 




THX2 




THX3 

y 




COMMON 

label 


FORTRAN 

variable 


Description 


/XIC/ 


AX1TN 

AX2TN 

BX1TN 

BX2TN 

TX1TN 

TX2TN 


First transfer function parameters 


DSPLAY Arrays 

The real-time simulation program uses specified arrays for displaying and/or 
changing the values of desired program variables. The desired program variables as 

defined in these specified arrays are assigned display addresses as shown in the following 
table: c 


DSPLAY 

Display 

Maximum number 

arrays 

address 

of elements 

TABLE (I) 

I 

199 

INTEG(I) 

I + 200 

99 

LOGIC (I) (not used) 

I + 300 

99 

ADC (I) 

I + 400 

99 

D AC (I) 

I + 500 

99 

LDISI(I) 

I + 600 

99 

LDISO(I) 

I + 700 

199 


This type of addressing is called "IN TABLES" addressing. For program variables 
not in the DSPLAY arrays, a type of addressing called "NO TABLES" addressing is used. 
The DSPLAY arrays are listed below with their associated FORTRAN variables and 
descriptions. The array elements are equivalenced to the FORTRAN variables, except 
where equality signs are indicated. 
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TABLE is a floating point number array. 


TABLE 


TABLE 

element 

FORTRAN 

variable 

Description 

1 

= A1SS 

“1SS 

2 

= AD1SS 

d lSS 

3 

= BISS 

^lSS 

4 

= BD1SS 

^1SS 

5 

= BETSI 


6 

= BDTSI 

% 

7 

= E 

E 

8 

= EDOT 

E 

9 

= H 

h 

10 

= HDOT 

h 

11 

= N 

N 

12 

= NDOT 

N 

13 

= OMEG 

n 

14 

= OMEGDT 

si 

15 

= P 

p 

16 

= PDOT 

p 


TABLE 

element 

FORTRAN 

variable 

Description 

17 

= PHI 

$ 

18 

= PHID 

<f> 

19 

= PSI 

* 

20 

= PSIDT 

* 

21 

= Q 

q 

22 

= QDOT 

q 

23 

= R 

r 

24 

= ROOT 

r 

25 

= THETA 

6 

26 

= THETDT 

6 

27 

= U 

u 

28 

= UDOT 

u 

29 

= V 

V 

30 

= VDOT 

V 

31 

= W 

w 

32 

= WDOT 

w 

33 

= AOS 

> 

o 

CO 

34 

= A1S 

A 1S 

35 

= BIS 

B 1S 
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TABLE 

element 


FORTRAN 

variable 


Description 


36 

37 

38 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 

51 

52 

53 

54 

55 


= THTR 


= OMEGD 


AOSO 


A1S0 


B1S0 


THTRO 


PHIO 


THETAO 


U KNOTS 


HDOTDO 


OMEGO 


GRWTO 


0 T r 

t 


Initial conditions 
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TABLE 

element 


FORTRAN 

variable 


Description 


56 

P0 

57 

Q0 

58 

R0 

60 

A1SSTEP 

61 

B1SSTEP 

62 

THTRSTP 

63 

A1SGN 

64 

B1SGN 

65 

THTRGN 

66 

AOSGN 

67 

PGN 

68 

QGN 

69 

RGN 

70 

UGN 

71 

VGN 

72 

WGN 

73 

HGN 

74 

PSIGN 




J 

>1 


> 



Initial conditions 


Control step inputs 


Scale factors for variables on time history recorders 
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TABLE FORTRAN 

element variable 


75 

PHIGN 

76 

THEGN 

77 

A1SSGN 

78 

B1SSGN 

79 

ALFFGN 

80 

BETFGN 

81 

BETSGN 

82 

BDTGN 

83 

DTCON 

84 

X10 

85 

X20 

86 

X30 

87 

X40 

88 

T10 

89 

T20 

90 

SF10 

91 

SF20 

92 

SF30 ^ 


Description 


Scale factors for variables on time history recorders 


Divisor on step size DT 


Control stick inputs, initial conditions 


Time at which to apply input, initial condition 
Time at which to remove input, initial condition 

Cockpit scale factors, initial conditions 
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TABLE 

element 


FORTRAN 

variable 


Description 


93 

SF40 


94 

SF 50 

) Cockp 

95 

SF60 


96 

SF70 


125 

= CT 

C 

126 

= CP 

Cp 

127 

= CQ 

C Q 

128 

= MU 


129 

= MTIP 

m tip 

130 

= AS PD 

a spd 

131 

= SHP 

S HP 

132 

= XAOS 

X AOS 

133 

= YCS 

Y CS 

134 

= XCS 

X CS 

135 

= XTR 

X TR 


Cockpit scale factors, initial conditions 
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INTEG is an integer number array. 


INTEG 


INTEG 

element 

FORTRAN 

variable 

Description 

1 

NSTICKO 

Denotes cockpit or potentiometer inputs 

2 

NT 

Determines print interval 

3 

IRUN 

Run number 

4 

NILAGO 

Denotes lag input 

5 

INTSCME 

Selects integration scheme 

6 

KNOTS 

Selects airspeed 

ADC is an input array from potentiometers or cockpit inputs. 



ADC 

ADC 

element 

FORTRAN 

variable 

Description 

1 

XAOS 

X AOS 

2 

YCS 

Y CS 

3 

XCS 

x cs 

4 

XTR 

X TR 

5 

YCSCOL 

Lateral stick trim adjustment input 

6 

XCSCOL 

Longitudinal stick trim adjustment input 

7 

XTRCOL 

Tail rotor pedal trim adjustment input 


DAC is an output array. Elements 1 through 24 are for the time history recorders 
Elements 25 through 40 are outputs to the cockpit. 

DAC 


DAC 

element 

FORTRAN 

variable 

1 

XAOS *XAOSGN 

2 

AOS *AOSGN 

3 

XCS *XCSGN 

4 

BIS *B1SGN 

5 

YCS *YCSGN 

6 

A IS *A1SGN 

rr 

XTR*XTRGN 

8 

THTR*THTRGN 

9 

P*PGN 

10 

Q*QGN 

11 

R*RGN 

12 

PHI *PHIGN *D PR 

13 

THETA *THEGN*DPR 

14 

PSI*PSIGN*DPR 

15 

ALFF*ALFFGN*DPR 

16 

BETF*BETFGN*DPR 


Description 


Normalized time history recordings 
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DAC 

element 


FORTRAN 

variable 


Description 


17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 


H*HGN 


U*UGN 


V*VGN 


W*WGN 


A1SS*A1SSGN 


BISS *B1SSGN 


BETSI BETSGN 


BDTSI BDTGN 


-COS(PHI) 


-SIN(PHI) 


-COS(THETA) 


-SIN(THETA) 


-COS(PSI) 


-SIN(PSI) 


COS(AA4) 


SIN(AA4) 


-COS(AA5) 


SIN(AA5) 


Normalized time history recordings 


Euler angle drivers 


Altimeter drivers 


Rate of climb meter drivers 


41 



DAC 

element 

FORTRAN 

variable 

Description 

35 

LIM(COS(AA6),- .886,.886) ^ 

Limited two minute turn drivers 

36 

LIM(SIN(AA6) ,- .500, .500) J 


37 

-COS(AA7) 

Bank angle indicator drivers 

38 

-SIN(AA7) J 


39 

LIM(OMEGRPM/ 500 . ,0 .0 , .9996) 

Limited RPM indicator driver 

40 

LIM(AS PD/140., 0.0, 0.9996) 

Limited airspeed indicator driver 


LDISI is a logical discrete input array where the descriptions are for true (.T.) 
values of the switches. 

LDISI 

LDISI Description 

element 

1-16 Data entry keyboard 

17 OPER (OPERATE) mode switch, integration of equations of motion 

18 HOLD mode switch, stops integration of equations of motion and holds all 

variables at present values 

19 RESET mode switch, initialization of variables (conditions for t = 0) 

20 TERM (TERMINATE) mode switch, terminates control at the program 

control console and transfers control to the graphic display unit 

21 CHANGE mode switch, allows changing of variable values 

22 SCAN mode switch, scans through the display addresses in conjunction 

with the subroutine SCANNER 
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LDISI 

element 

23 

27 

29 

30 

31 

32 

33 

37 

38 

41 

42 

43 

44 

45 

46 

47 

48 
55 


Description 

RELEASE mode switch, releases CHANGE or SCAN mode switch 

ERASE mode switch, erases real-time disk file 

IDLE mode switch, uses minimum computing time 

PRINT mode switch, prints output of previously stored variables 

READ mode switch, selects printout storage format 

RELEASE mode switch, releases ERASE, IDLE, PRINT, or READ mode 
switch 

Dynamic check doublet inputs on A^g, Bjg, or 

Trim circuit switch 

Selects second set of printout variables 

Stability augmentation system (SAS) switch 

Cockpit instrument test 

Cockpit mode control 

Cockpit tie-in 

ADC inputs 

AUTOTYPE, a logical variable used with the typewriter 
TYPESW, a logical variable used with the typewriter 
"IN TABLES" addressing 
Cockpit SAS switch 
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LDISO is a logical discrete output array. Elements 61 through 99 are used to turn 
the white indicator lights on and off. Elements 101 through 110 are recorder event mark- 
ers. Elements 111 through 118 are used to turn the red indicator lights on and off. Only 
the following elements are used. 

LDISO 

LDISO Description 

element 

61 Variables are trimmed 

65 SAS is activated _ 

73 a F out of range 

74 /3-p out of range 

75 a w out of range 

76 Q! H S out of range 

77 /3y g out of range 


78 

U'T' 

1 'kY 

out of range 

79 


out of range 

80 

Cv 

y tr 

out of range 

81 

X AOS 

is trimmed 

82 

Y cs 

is trimmed 

83 

x cs 

is trimmed 

84 

X TR 

is trimmed 


102 Timing marker on recorder no. 1 
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LDISO 

element 

Description 

110 

Timing marker on recorder no. 2 

111 

Trim unsuccessful 

112 

Aqs out of range 

113 

A-^g out of range 

114 

Bjg out of range 

115 

$Tr out of range 

116 

Automatic hold 

117 

Negative altitude 

118 

Engine power out of range 


Real-Time System Subroutine Descriptions 
The real-time system subroutines and brief descriptions of their functions are as 


follows: 


Subroutine 


Description 


APRINT 


Provides an alternate entry point to subroutine RTMODE which returns 
control to subroutine RTMODE for further processing 


ARE AD 


Provides an alternate entry point to subroutine RTMODE which returns 
control to subroutine RTMODE for further processing 


ATERM Does final processing and halts execution 


CYCLE Sets up return address from subroutine RECYCLE 

DATABLX Specifies the variable arrays for subroutine DSPLAY 
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Subroutine 


DAYTIM 

HALT 

IGRATE1 

INOUT 

LOSTIME 

NAME CRT 

NM218 

PLAYBAK 

PRINTER 

READOUT 

READY 

RECORD 

RECYCLE 

RELEASE 

RTMODE 
RT ROUTE 


Description 

Provides date and time of day 

Signals completion of real-time portion of program 

Integrates variables in array DERINT(2,J) and stores in array 
DERINT(1,J) (J = 1,2,...,NEQ) 

Sets up arrays for input/ output conversion 

Sets address to which control will be returned in the case of iteration 
time exceeded 

Identifies and initializes usage of the CRT console 

Initializes usage of typewriter 

Plays back data recorded on real-time disk one frame at a time 

Routes MFI disk file to printer 

Specifies variables to be recorded and frequency of recording for the 
real-time disk file 

Signals that the program is ready for real-time operation 

Records variables as specified in subroutine READOUT 

Signals end of a cycle in the operating loop and returns to address 
specified by subroutine CYCLE 

Releases a real-time data file that has been opened by subroutine 
READOUT 

Real-time mode control 

Indicates dynamic printing will occur during the job run 
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Subroutine 


Description 


SCANNER 

SYNCH 


TYPEVAR 

XDSPLAY 


Increments display address 

Issues a dayfile message that indicates the current instruction that was 
processing when synchronization with real time was lost 

Types data displayed on digital decimal display unit 

Initializes data entry keyboard and digital decimal display unit, routes 
program listing 


Program Subroutine Descriptions, Flow Charts, and Listings 

This section contains a brief description of the main program and each of the sub- 
routines. A simplified flow chart and a FORTRAN listing are given at the end of each 
description. The principal aircraft equations used in the simulation program and the 
references for their sources are contained in appendix A. 

HELIC 

Helic is the main program from which all the subroutines are called. HELIC 
contains the COMMON and EQUIVALENCE statements needed for communication with the 
various subroutines and the other FORTRAN declarative statements such as DIMENSION, 
REAL, and LOGICAL. The main parts of the HELIC program are the real-time system 
software, the program constants, the variable initializations, the aircraft trim circuit, 
and the calls to the subroutines. 
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PROGRAM MEL l C (fNPUT.201, OUTPUT-201) 


TYPE STATEMENTS ANO DIMENSIONED VARIABLES 

REAL I HS « I VS * IROT* 1W* I XX . I YY , IZZ* LA, LF , LMR, LRH , MASS* 

X MRH, MU • MA , mF« MTIP. N, NDOT • NA , NT, NO 

LOGICAL AUTOH, AUTOTyP* 1NTABLS* LOIS1, LOISO. LOGIC* 

X TYPESW.VARCHmG, mreset , mhold. moper 

DIMENSION IVARBUF<5), LOGIC(Ol), TIME(2>. MESAGE ( 2 ) • 

X XTRIMtll) , FTP I M ( I 1 I • 0VECU2) 

real time control communication 

INTEGRATION communication 

COMMON 

X /CHECK / TABLE! 150)* INTEGI09) 

X /INTCOMM/ T . DT * I NT * NEO * 

X ISCHEME. DER I NT (2*12) 

X / I NT l NTR/ I NTfRN (5*12) 

X /REALTIM/ ADC ( 32 I * DAC<64), LDISU108). LDTS0(196)« 

X NOPeR • NHOLD , MRESET * NTERM , 

X NPR I NT « NREAD 

SUBROUTINE COMMUNICATION 


X 

/ACCEL/ 

AX 


AY . 

AZ 







X 

/AOV C / 

OMEG 


WIM , 

VT . 

QV 

« 

QL 

* 

ASPD 


X 

/AOV P / 

ADVPl 


ADVP2 








X 

/ANGSAS / 

PHIL 


THETAL • 

PS IL , 

PHI DL 


THETDL 


PS I DL 


X 

/cntl c / 

AOS C 


Al SC , 

B ISC . 

THTRC 


Al SCL 


Bl SCL 

* 

X 


thtdcl 










X 

/CNTL p / 

AOScO 


A 1 SCO , 

B 1 SCO 

THTRCO 


XAOSG 


YCSG 

* 

X 


XCSp 


XTRG , 

XAOSP * 

YCSR 


XC SR 


XTRR 


X 

/CNTSAS / 

Ai SrS 


BlSCS . 

THTRCS 







X 

/COEF / 

CT 


CP ♦ 

CO 

MU 


MT IP 


VS 


X 

/CONTRL/ 

AOS 


Al S 

BIS . 

THTP 


OMEGD 




X 

/CPIC / 

AC I 


BC t * 

DC I * 

XRCM1 


yrcmi 


XRTM 1 

• 

X 


YRTMl 


XPBM1 , 

ypbmi , 

CSK4 


CSKfe 


CSK7 

* 

X 


F05 < i© J 










X 

/ENG C / 

OMAX 


OE . 

SHP . 

GRWT 

* 

OMEGDT 



X 

/ENG P / 

IROT 


CKE1 








X 

/error / 

NERP 










COMMON 












X 

/EUL CS/ 

CO Sf I 


StNFI » 

COSTH . 

S I NTH 


cossi 

, 

SINS! 


X 

/FCNNAME/ 

CLF J 


CMF 1 . 

CNF 1 . 

CXI 


CX2 

* 

C YTR 

* 

X 


CYF 


CZF , 

cxw • 

CZW 


CYVS 

, 

CZHS 

, 

X 


CDS I Y 


CLSIY , 

Dw . 

GEV 


GEH 




X 

/F N M C/ 

XA 


YA . 

7 A • 

LA 


MA 

, 

NA 


X 

/F N M P/ 

DX 


DXHS , 

OXTR , 

DXVS 


DZ 

• 



X 


DZTp 


DZVS . 

DXW « 

DZW 


IR 




X 

/FTRMPLS/ 

EQLMR 


EOAOSS. 

EQWlM , 

HDOTD 






X 

/FUNCS / 

FOOi (20 ) 

F002 ( 20 ) , F003< 

4 ) • F 00 4 (14), 

F 005 (14) 


X 


F006 ( 28 > 

F007 < 1 4 ) » FOO8<20>, F009 ( I 4 ) . 

F01 0(28) 


X 


F01 l (35 > 

F013C07)* FO 1 4 ( 07 ) * F015<07>, 

DO! 8 (275 


X 


DOO i (09 > 

D002 < 09 ) * 0003(09), D004(09>» 

0005 (09 ) 


X 


DOOfi (09)< 

D007<09). 0008(09). 0009(09)* 

DO 1 0(09) 


X 


DO 1 i (09 ) 

DO 13(09). 0014(09). DO 1 5 ( 09 ) , 

0016(18) 


X 


DO I 7 ( 1 0 

) 

FO 1 6 < 35, 10). F 0 1 7 ( 35 « 

10), FO 

10(7,3.4) 

X 

/rus C / 

XF 


YF ♦ 

ZF 

LF 


MF 

, 

NF 

• 

X 


xw 


W I WM 

ZW « 

ALFF 


BE TF 

, 

ALFW 


COMMON 












X 

/FUS P / 

CLF? 


CMF 2 , 

CNF2 * 

CKRF 


CKWI WM 

, 

1 W 

, 

X 


VLF 


YMF . 

VNF « 

S XF I 


SXF2 

, 

S YF 

, 

X 


SZF 


sxw . 

SZW 







X 

/GPARAM/ 

DELT 


RHO . 

G « 

PIE 






X 

/MROT C/ 

FXR 


fyr 

LMR 

LRH 


MRH 

, 

OMR 

, 

X 


AOSS 


alfsi Y, 

UTS1Y 

BETS I 


BOTS 1 

« 

Al SS 

. 

X 


AO 1 eS 


R 1 SS ♦ 

RD1 SS 







X 

/RCOEFF / 

CKL (20 ) 


CK D ( 20 ) 

CKQ (20 ) 

CK DL ( 2 0 > « CKQL 

(20 ) , 


X 


CKM (20 ) 


Y ( 20 ) 

WF ( 20 ) 







X 

/ROTCON / 

YPE (20 ) 


YTW<20 > 

W IF 120 ) . 

BOV N 






X 

/ROT P / 

T W I FT 


F MRS 

CKAO . 

CKABJ 

« 

CKAB2 

• 

CKA83 

, 

X 


NR A 0 


NAZ . 

NB . 

RB 

* 

EEH 




X 

/LAGSTK / 

N1 LAG 


NSTICK 








X 

/MODEC / 

mreset 


MHOLD » 

MOPER 







X 

/SAS P / 

A J Cl 


A l C2 

A1KJ 

A 1 <2 

, 

A] K3 

, 

A 1K4 

, 

X 


BI Cl 


R 1 C2 . 

R 1 K I , 

Bl K2 

• 

PI K3 

, 

B 1 K4 

, 

X 


THC 1 


THC2 

THK 1 , 

THK 2 

, 

THK3 

, 

THK4 


COMMON 












X 

/scales / 

SF 1 

• 

SF2 , 

SF 3 , 

SF 4 

. 

SF5 

, 

SF6 

• 

X 


SF7 










X 

/SCPSIR / 

SINf IR ( 5r 

). COSS I R ( 50 1 







X 

/SHIP P/ 

1 XX 

• 

I YY • 

IZZ , 

MASS 






X 

/STICKS / 

YA O c 

. 

YCS 

xcs 

XTR 

, 

DPR 

, 

RPD 


X 

/ST I NPUT/ 

XI 

• 

X 2 

X 3 • 

X4 

, 

T 1 

. 

T 2 




X 

/T AfL C / 

ALFhS 

. 0 TVS 

ZHS 

YVS . VTR , 


X 


WHS 

. YTR 

DELE 



X 

/TAIL P/ 

CKTRl 

. CKTR2 

INS . 

I VS . SZHS . 

S YVS 

X 

/TRIMAT/ 

XMAt ( 1 1 

♦12). EM AT (12.13). 

TRSUM ( 1 2 > • NEON 


X 


NEQP1 

* NEOP2 , 

XMATG ( t 1 ) 


X 

/VALUES / 

A 1 CFt 

. A 1 CS2 

BICSl • 

SICS2 • thcsi . 

THCS2 

X 

/xi c / 

A1 XJ 

* A 1 X2 

A I X3 . 

01X1 .01X2 . 

B I X3 

X 


THX ) 

. THX2 

THX3 « 

AX1TN . AX2TN • 

ex 1 TN 

X 


9X2 TN 

. TX1TN 

TX2TN 




EQUIVALENCES 


EQUIVALENCE! INTEG ( 1).NST1CK0>. 

X { INTEG ( 3). IRUN >. 

X ( I NTEG ( 5 ) « INTSCME ) » 


t INTEG! 2 ) * NT ). 

(INTEG! 4 ) ♦ N I LAGO ) » 
(INTEG! 6 ) * KNOT S > 


EQUIVALENCE (LDJSl (46 AUTQTYP ) * (LDISI (47) . TYPES W ) . 


X 

(LDISI (40) 

. INTA8LS > 




EQUIVALENCE 

X 

(OERINT ( 1 t 

I ) . P ) . 

(DERINT (2. 

1 ) 

• PDOT ) 

X 

(DERI NT ( 1 , 

2) .0 ) « 

(DERINT (2. 

2 > 

* QDOT ) 

X 

( DERINT ( 1 , 

3 > *R ) « 

(DERINT (2. 

3 ) 

* RDOT ) 

X 

( DERINT ( l , 

4). PHI ). 

(DERINT (2 . 

A ) 

•PHIO ) 

X 

(OERINT ( 1 . 

5 > .THETA )« 

(DERINT (2. 

5) 

♦ THETOT ) 

X 

(OERINTU , 

6 ) « PS I ). 

(DERINT (2 . 

6 ) 

.PSIDT ) 

X 

(DERINT ( 1 , 

7 > • U ). 

(DERINT (2 . 

7) 

tUOOT ) 

X 

(DERINT 1 I , 

8 > . V ) . 

(DERINT (2. 

8 ) 

.VDOT ) 

X 

(DERINT! 1 , 

9 ) .W ) ♦ 

(DERINT (2. 

9) 

» WDOT ) 

X 

(OERINTU , 

1 0) .N ). 

(DERINT (2. 

10 ) 

• NDOT ) 

X 

(OERINT ( I * 

1 l ) «E ) . 

(DERINT (2 . 

1 I ) 

.EDOT ) 

X 

( DERINT ( I , 

1 2 ) .H ). 

(DERINT (2. 

12) 

.HOOT ) 


EQUIVALENCE ( AOSO, 

TARLE (40 ) I . 

(A1S0 , 

TABLE (41 ) ) 

X 

(01 so. 

TAflLE (42 ) ) , 

(THTRO . 

TABLE (43 ) > 

X 

(UO . 

TAnLE (44 > ) . 

(VO . 

TABLE (45 ) ) 

X 

( wo . 

TArLE ( 4fi > ) . 

(PSIO . 

TABLE (47 ) ) 

X 

(PHIO. 

TARLE (40 > ) * 

(THETAO. 

TABLE (49) ) 

X 

(UKNOTS. TABLE < 50 > ) ♦ 

(HDOTDO . 

TABLE (51 ) ) 

X 

( HO . 

TARLE (52 > J . 

(OMEGO , 

T ABLE ( 53 ) ) 

X 

(CGO . 

TARLE (54 ) ) . 

(GRWTO , 

TABLE (55 > ) 

X 

(PO , 

TARLE (56 ) I « 

(OO . 

TABLE <57) ) 

X 

( RO , 

TARLE (50 ) ) 




EQUIVALENCE < A 1 SSTEP. 

TABLE (60) ). 

(01 SSTEP 

TABLE (61 ) ) . 

X 

( THTRSTP . 

TABLE (62) >« 

(AtSGN 

TABLE (63) ) . 

X 

(01SGN . 

TABLE (64 ) ) . 

(ThTRGN 

TABLE <65 ) ) . 

X 

(AOSGN . 

T ABLE ( 66 ) )« 

(PGN 

T ABLE ( 67 ) ) . 

X 

(QGN • 

TABLE ( 68 t ) • 

f RGN 

TABLE (69 ) ) . 

X 

(UGN * 

T ABLE ( 70 ) ). 

( VGN 

TABLE (71 > ) . 

X 

(WGN . 

TABLE ( 72 ) ) . 

(HGN 

TABLE (73 ) ) . 

X 

(PSIGN . 

TABLE ( 74 ) ) . 

(PHIGN 

TABLE (75 ) ) » 

X 

(THEGN . 

T ABLF < 76 ) ) . 

( A1 SSGN 

TABLE (77 ) ) . 

X 

(B1SSGN . 

TABLE! 78) ). 

(ALFFGN 

TABLE (79 ) ) . 

X 

(BETFGN » 

T ABLF ( 80 ) ) * 

(BETSGN 

TABLE (81 ) ) . 

X 

( BDTGN . 

TABLE ( 82 ) > . 

( DT CON 

TABLE (83 ) ) . 

X 

(XIO * 

TABLE (84 ) > . 

(X20 

TABLE (85) ) . 

X 

( X 30 * 

TABLE (86 ) >« 

(X40 

T ABLE (B7 ) ) » 

X 

( T1 0 . 

TABLE (88 ) > * 

IT20 

T ABLE ( B9 I ) . 

X 

( SF 1 0 , 

T ABLE ! 90 ) ) . 

(5F20 

TABLE (91 ) ) . 

X 

(SF30 . 

TABLE (92) ) . 

(SF40 

T ABLE (93 ) ) . 

X 

!SF50 . 

T ABLF ( 9A ) ) ♦ 

(SF60 

TABLE (95) ). 

X 

(SF70 , 

T ABLE ( 96 ) > 




C** 

c 


DATA ME SAGE/ 2HON. 3hOFF / 

INITIAL SET-UP OF PFAL TIME SYSTEM 

CALL NM21 a ( 5LTYPER ) 

CAI *. CYCLE ( 90006S ) 

CALL LOST (ME (77777S I 


C** CALL READOUT SETS UP THE REAL TIME RECORDING FILE 


NT « 
CALL 

32 

READOUT (4 .NT . T 


ASPO . 

GRWT 

. CG 

CALL 

READOUT ( 4 .NT • H 


HOOT . 

OMEG 

. OMEGDT 

CALL 

REAOOUT ( 4 .NT . U 


V . 

w 

. ALFF 

CALL 

PEAOOUT (4 ,NT. UOOT 


VDOT 

WOOT 

. BETF 

CALL 

READOUT ( 4 .NT . P 


0 . 

R 

. DELE 

CALL 

READOUT ( 4 .NT . PoOT 


QDOT 

RDOT 

. LMR 

CALL 

READOUT ( 4 , NT . XA 


Y A * 

ZA 

* OMR 

CALL 

READOUT ( 4 .NT ♦ LA 


MA . 

NA 

* FXR 

call 

REAOOUT (4 .NT . Ph I 


THETA , 

PS I 

. FYR 

CALL 

READOUT (4 .NT , Ph I D 


THETOT, 

PSIDT 

. QE 

CALL 

READOUT ( 4 .NT , AqSS 


Al SS 

BISS 

. WIM 

CALL 

READOUT (4 , NT , AOS 


A 1 S 

BIS 

, THTR 

CALL 

READOUT (4 .NT . XAOS 


YCS 

XCS 

. XTR 
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CALL READOUT U,NT, N 
CALL READOUT (4. NT. QvAX 
CALL READOUT (2.NT. Qv 


EDOT J 
VT ) 


C** FORMAT NO. SUPPLIED TO SUPERVISOR FOR PRINTING JOS CARD FROM ME FILE 
C 

CALL RTROUTE < MF I ,900s4S I 
C 

C** CALL 1NOUT SETS UP I/O TO REQUESTED EQUIPMENT 
C 

CALL I NOUT ( ADC , 8 » DAf,48« LDISI.60. LDISO.10O) 

CALL DATA0LX < TABLE . 1 MO. INTEG.9. LOGIC. Ol. ADC.0. DAc. 48 . 

X LDIST. *0. LOT SO. 180) 

CALL XDSPLAT (LDIS I . LDISO* VARCHNG. I TYPE * IVARBUF. InTABLS) 

C 

C** ASSIGN STATEMENTS SET uP MODE CONTROL RETURNS TO PROGRAM 

C 

ASSIGN 90001 TO NOPEP 
ASSIGN 90002 TO NHOLn 
ASSIGN 90003 TO NRESfT 
ASSIGN 90004 TO NTERm 
ASSIGN 90014 TO NPRInT 
ASSIGN 9001 5 TO NRE Ar> 

C 

C** LOGICAL discretes set TO false 

c 

00 in I NO * 1.0 
10 ADC (I NO) ■ 0. 

DO 20 IND - 1.40 
20 DAC f IND) ■ 0. 

DO 30 IND * 1.60 
30 LDIS I { I ND)« *F. 

DO 40 IND ■ l . J 00 
40 LOISOI I NO ) ■ .F, 


READ 90. FOO0 
READ 90. E009 

READ 90. FOtO 
READ 98. FOl 1 
READ 90. FOl 3 
READ 98, FOl 4 
REAO 90. FOl 0 
READ 90, FO 1 6 
READ 98. FOl 7 
READ 9B * EO 1 8 
98 EORMATCTEIO.O ) 

* CONSTANTS ANO INITIAL PARAMFTERS 

* INITIAL LOGICAL VARIABLES 

LOGIC'Ol l a ,F . 

VAPCHNG » • F « 

MRESeT ■ .T* 

MHOLO « .F* 

MOPER ■ .E. 

* INITIAL INTEGER PARAMETERS 

I COUNT ■ 1 
INTSCME - 5 


INITIAL DOUBLET INPUTS 
AISSTEP - O. 

BlSSTEP » 0. 

THTRsTP ■ 0. 


SAS constants 
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AXlTN * 0 . 

AX 2 TN * 0 . 

BXITn * o. 

BX 2 TN - 0 . 

TX 1 TW ■ 0 , 

TX 2 Tn ■ 0 . 

C 

c* programmed control STICk- DEFLECTIONS 
XI 0 * 0.0 

X 20 ■ 0.0 

X 30 » 0.0 

X «0 ■ 0.0 

T I 0 * 0.0 

T 20 * 0.0 

C 

C** I.C* S ON INTEGRATED VARIABLES 
PO « 0. 

00 * 0 . 

RO » 0 . 

PSIO * 0. 

UO » 0. 

vo * o • 

wo -o. 

NO 3 0 . 

EP - o. 

HO ■ 3000 . 

REVTRAD = 30 . /PIE 

OMEGO = 33 * 9 *REVTRaD 

c 

C** IC’S ON DERIVATIVES used in TRIM 
PDOT ■ O. 


ODOT 

■ 0 . 

ROOT 

* 0, 

PH ID 

* 0 , 

thetdt 

* 0, 

PSIDT 

* 0. 

UPOT 

* o , 

VPOT 

* 0 . 

WDOT 

* 0. 

NDOT 

* 0. 

FOOT 

* 0 . 

HOOT 

« 0, 

recorder gains 

AOSGN 

* 1 . 

A 1 SGN 

■ 1 0 . 

R 1 SGN 

■ 10. 

THTRr,N 

* 3. 

X AOSGN 

= . 066666666* 

XCSGN 

= .10 

YCSGN 

= .10 

xtrgn 

* .to 

PGN 

= 2.3 

OGN 

* 2.3 

RGN 

* 2.3 

PH I GN 

* .023 

THEGN 

- .03 

PS IGN 

■ .03 

UGN 

* .003123 

VGN 

■ .023 

WGN 

■ .023 

HGN 

- .0001 

alffgn 

■ .03 

retfgn 

» .05 

A 1 SSGN 

- 10. 

B 1 SSGN 

* 20. 

bf t sgn 

* 10. 

bdtgn 

» .333333 

;OCKP!7 constants 

NILAgO 

* 0 

NST TOKO 

* 0 

H7 

* 32./1024. 

tlag 

= 1.0 

AC I 

* 1.0 - (TLAG/H7 

BC I 

* ( TLAG/H7 5 * ( I .0 

DC I 

- EXP I-H7/TLAG) 

CSKA 

* PIE/300. 

CSKft 

* -10. 

CSK7 

« 20. 

XPSM 1 

■ 0,0 

YPBM ] 

« 0.0 

XRCM t 

- 0.0 

YRCM 1 

= 0.0 

xrtm i 

« 3. I 416 

YRTM I 

■ 3. I 4 16 

SE10 

■ 16.625 

SF20 

- 8.00 


EXP (-H7/TLAG) ) 

*■ H7/TLAG )*EXP (-H7/TLAG I > 
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SF30 

-8.00 

SF40 

B* 00 

SF50 

B* 00 

SF60 

-B.OO 

5F70 

4.00 

F05< 1 ) 

0.0 

FOB ( 21 

B.21 70 

FOB ( 3) 

10. BOO 

FOB ( 4) 

t 3 # 7 34 

FOB ( B) 

19.183 

FOB ( 6 > 

23.303 

FOB ( 7) 

27.833 

FOB ( 0) 

33.B67 

FOB ( 9 ) 

37.350 

FOB ( 1 0 ) 

4 t .750 

F OB ( 1 1 ) 

44.350 

FOB ( 1 2 > 

61 .990 

FOB ( 1 3 > 

77,930 

FOB < 1 4 ) 

83.834 

FOB ( 1 B ) 

106.784 

F OB < tft) 

1 16.600 

FOB( 17) 

128,336 

F OB < 1 B ) 

128.336 

AOSS 

,00Ol 

At SS 

0 . 

B 1 5S 

0 . 

BETS! 

0 . 

CGO 

1 93.85 

CNTERR 

.01 

OPR 

57.2957795 

DTCON 

! • 

DTO 

.03125 

GRWTn 

8822.8 

HDOTnO 

0 . 

lmr 

9000. 

QF 

0 . 

PPD 

0.01 745329 

WIM 

o. 

DO 75 J» 

1.11 


75 XM AT f; f J ) - .001 

C 

C** COMPUTE ROTOR BLADE PARAMETERS 

C 

00 50 J • I * NR AO 
YPE(J) • Y(J> + PFH 
YTW(J) ■ TWIST#YPE< J )/RB 

W I F l J ) * 3,685*Y( J)/f ( 0MEG0/RFVTRAO ) *RB*RB ) 

BO CONTINUE 

C 

c*+ COWPUTF sines ANO COSINES around AZIMUTH 

c 

DO 60 J ■ 1 « NAZ 

PS I R ■ 2* *P| E*FLOAT { J-l ) /FLOAT (NAZ I 

COSSIR(J) - COS ( PS f R ) 

SINSIR(J) * SINtPSjR) 

60 CONTINUE 
c 

BOVN * FLOAT ( NB ) /PLOAT (NAZ ) 

C 

C** BEGINNING OF REAL T I MF LOOP 
C 

CALL READY 


C** BEGINNING OF RESET LOOP 
C 

90003 CONTINUE 

1F(MhOLD> GO TO 9000? 
t F ( MOPER I GO TO 9000 | 

C 

C** SELECT INITIAL CONDITIONS 

GO To < 90! .902 .903. 904 .905 . 90ft ) « 


901 CONTINUE 

C** HOVER INITIAL CONDITIONS 

C 


AOSO ■ 
A 1 S 0 a 
B1 SO * 
THTRO ■ 
PH 10 * 

THFTaO* 
UXNOt S * 


0. 2739429 
-0*03121 639* 
-0 * 063409593 
0* 1 392309 
-0*01 929372ft 
-0*0621 521 22 
0 . 


vo = o. 

wo « o* 

GO TO 906 


KNOTS 


902 CONTINUE 


C** 33 KNOTS INITIAL CONDITIONS 
C 

AOSO * *2462616 
A1SO --.024696768 
01 SO --*044236309 
THTRo - *0809095015 
PMIO --.014007214 
TMET AO--, 073933757 
UKNOTS- 33.15571344 
VO ■ 0*00 

WO --4.14B259 

GO TO 906 
C 

903 CONTINUE 

C** 80 KNOTS INITIAL CONDITIONS 
C 

AOSO * .2396654 

AISO *-*016621463 
8 1 SO --*011030726 
THTRo - .029313655 
PHIO --.017073581 
THET AO--. 09730 7773 
UKNOTS- 79.928952 
VO • 0.00 

WO *-13.18009 

GO TO 906 
C 

904 CONTINUE 

C** 105 KNOTS INITIAL CONDITIONS 

C 

AOSO • .2591912 

A ISO *-.010925199 
01 SO ■ .021399861 
THTRO * .031807306 

PHIO --,026224227 
THCTaO--, ] 1 58838 
UKNOTS- 105.1509769 
VO - 0,00 

WO --20. 68070 

GO To 906 
C 

906 CONTINUE 

0** 120 KNOTS INITIAL CONDITIONS 

C 

AOSO - .2791526 
AISO *-.021712153 
01 SO » .045708954 

THTRo - ,038308344 
PHIO --.034026789 
THETAO--. 1 291323 
UKNOTS- 120.0 
VO - 0.00 

WO --26.33423 

C 

906 CONTINUE 

UO * UKNOTS* 1.689 
KNOTS • 6 

C 

AOSC - AOSO 

A 1 SC • AISO 

01 SC - 01 SO 

THTRC - THTRO 

C 

C»* TRIM CIRCUIT 

C** IF FUNCTION SWITCH 5 I? ON. ITERAT E CONTROL DEFLECT I ONS 
C** UNTIL STATE VARIABLE DERIVATIVES ARE ZERO. 

C 

IFILDISI (37) ) GO TO 80 

I COUNT ■ 1 

5C0UNT ■ I 

C2 - .01 

NT 1 * 1 

NITER - 0 

LD I So ( 6 1 ) - .F. 

LDI SO Mil) * .F. 

GO TO 100 
80 CONTINUE 

IF (MCOUNT.GT.l ) GO TO 1 00 
1 000 CALL HALT 

LD 1 50(61 ) * *f. 

• 9 l SO ( 11 n - .F. 

- 10 


"N + 1 



r 

CALL SETUP A (XTRIM,THFTAO.PHlO.AOSO*4lSO,01SOiTHTRO.WrM, aoss. 
1 A1 SS «B 1 «S *LMR * 1 > 

1005 CALL XMATRJX(XTRIM) 

C 

C** CALCULATE THE EMATRIX 

r 

00 1010 NCOL * I.NEOol 
CALL XVAL (XTRIM* NCOL ) 

CALL SETUPA ( XTR I M , Th F T A . PH I t AOS . A 1 S .B 1 S . ThTR . W I M , AOSS . A 1 SS • 

1 BISS «LMo .2 ) 

CALL TR1H0ER 

CALL SETUP2 (FTP I M * UDO T « VOOT * WDOT »PDOT .QDOT ,PDOT * EOW I M, EOAOSS 
l Apl SS .BDl SS.EQLMR) 

CALL EVAL (FTP1M.NCOL ) 
lOlO CONTINUE 
C 

C*» FIND THE MINIMUM SUM OF F-MATRIX COLUMNS 
C 

SUMMIN * TRSUM ( 1 ) 

LL * 1 

DO 1020 HT*2,NF0P1 

IF ( SUMM I N* L T • TRSUM {•<■)) GO TO 1020 
LL ■ K 

SUMMIN * TRSUM <K| 

1020 CONTINUE 
C 

C** FIND THE MAXIMUM SUM 

C 

1030 SUMMAX » TRSUM ( 1 ) 

XX * 1 

DO 1040 U»2*NEQP| 

IF ( SUMMAX .GT • TRSUM ( J ) 1 GO TO 1040 
XX * J 

SUMMAX ■ TRSUM(J) 

1040 CONTINUE 
C 

C** SET UP MATRIX EQUATION - EMAT*QVEC = EMAT + 1 

C 

EMAT < 1 ,NE0P2 1 - t *0 
1045 DO 1050 J * 2 » NE OP 1 

1050 EMAT ( J.NEQP2 ) ■ ( 1 • -r 2 ) *EM A T ( J • LL 1 

C 

C** SOLVE MATRIX EQUATION FOR QvEC AND COMPUTE NEW TRIM VECTOR 
C 

CALL S I MEQA (EMAT .NEQP1 . QVEC . XM AT . XTR I M ) 

CALL SETUPA (XTRIM, ThETA .PHI ♦ AOS . A 1 S .B 1 S , ThTR « W I M . AOSS . A t SS . 
1 BISS . LmR .2) 

CALL TRIMDER 

CALL SETUP2 (FTRIM.UDoT* VDOT , WDOT , PDOT . QDOT , ROOT . EO W I M . EG AO SS 
1 A01 SS »BOl SS.EQLMR J 

SUMNfW * 0. 

DO 1070 J ■ 1 * NE ON 

1070 SUMNEW * SUMNEW + FTOtM(J>**2 

IFfSuMNEW *LT. TRMTOL » GO TO 1120 
MCOUNT ■ MCOUNT + l 
IFfMrOUNT • GT • 75 1 GO TO 1135 
IF ( SUMNEW. LT.SUMMAX > GO TO 1090 
1080 C2 ■ . 5*C 2 
NT 1 « 2 

IF (C2. GT. 1 .E-51 GO TO 1045 
GO Tn 1140 

1090 IF (NT1.EQ.21 GO TO l 095 
C2 • AMIN! (2.*C2. 1 .O) 
l 095 NT 1 ■ 1 

C 

C*4 REPLACE WORST VALUES WJTH NEW FUNCTIONS AND ITERATE 
C 

TRSUM (XX) • SUMNEW 
DO 1 1 00 I ■ 1 .NEON 

XMAT < I .XX > » XTRIM ( I ) 

1100 EMATM + l.XX) » FTRT M ( 1) 

!F (SUMNEW. GT.SUMM1N ) GO TO 1030 
LL * XX 

SUMMIN « SUMNEW 
GO TO 1030 
t 120 CONTINUE 
C 

C»* WHITE LIGHT t ON INDICATES VARIABLES ARE TRIMMED 

C 

LDISO (61 1 - *T. 

1130 GO To 1160 

C 

C** RED LIGHT 1 ON INDICATES UNSUCCESSFUL TRIM 

c 

1 l 35 LD I SO ( III) » • T . 

GO TO 1160 
1 140 CONTINUE 


C2 - *01 

NT 1 - 1 

CALL XVAL (XTR IM.LL > 

NITER « NITER + 1 
IF (NITER, LE. 2) GO TO 1005 
LDrSO( 111) - , T • 

1160 CONTINUE 

XTRIm (11) * LMR 

CALL SETUP A < XTR I M, THpT AO »PH I 0 * A OSO # A 1 SO ,81 SO » THTRO , W ] M , AOSS , A 1 SS , 
1 BISS ,LmR * 2 ) 

C 

WO - w 

OELT - DT 

CALL READY 
C 

100 CONTINUE 



XAOSt 

* 

( AOS0 

• DPR - 

AOSCO 

) /XAOSG 


XCST 

» 

( B I SO 

• DPR - 

B 1 SCO 

>/XCSG 


YCST 

* 

(A1SO 

• DPR - 

A 1 SCO 

)/YCSG 

C 

XTRT 

* 

(THTR0*DPR - 

THTRC 0 ) /XTRG 

c 

c*« 

IF FUNCT 

ION SWITCH 13 

IS ON ACCEPT 


C** CONTROL STICK INPUTS FROM COCKPIT 
C 


c 

c 

c 

c 


IF (LOl SI (45 ) ) 

AOSC » 

AOSO 

A 1 SC * 

A 1 SO 

01 SC 

B 1 SO 

THTRC * 

THTRO 

I I 0 continue 



AOS * AOSC 
A 1 S * A 1 SC 
BIS - 0 1 SC 
THTR ■ THTRC 


AUTOH « *F. 

NFpR * 0 

1NT -0 

I PR ImT « 1 


c 


C** ISCHEMES*. I -2ND RK, ?«4TH RE . 3-2NO AM, 4«4TH AM, 5 ■ 1 ST A0 
I SCMwMf * INTSCME 
IF ( INTSCME,EQ,5 ) I SCHEME = 3 
C 


CG 

■ CGO 

DT 

* dto/dtcon 

OELT 

» OT 

DX 

■ (CG - 200 # )/\: 

GRWT 

■ GRWTO 

HDOTn 

“ H00TO0/60* 

MASS 

» GRWT /G 

ONEGO 

■ OMEGO/REVTRAD 

OMEG 

- OMEGD 

T 

• 0, 

VS ■ 1100, 

p 

■ PO 

0 

■ 00 

R 

« RO 

PHI 

- PH I 0 

theta 

* THETAO 

PS I 

* PS 10 

u 

■ uo 

V 

* VO 

w 

■ wo 

N 

- NO 

E 

« EO 

H 

* HO 

NILAO 

* NILAGO 

NSTIcK 

* NS T I CKO 

Xt 

* X l 0 

X? 

« X20 

X3 

* X30 

X4 

* X 40 

T 1 

- T I 0 

T2 

* T20 

SF1 

- SF1 0 

SF2 

* SE 20 

SF3 

* SF 30 

SF4 

■ SF 40 

SF5 

* SF50 
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SF6 


SF 60 

SFT 


SF70 

A1 SCL 


A 1 SC 

Bt SCL 


BlSC 

THTRCL 


THTRC 

PHIL 


PHI 

THETAL 


theta 

PSIL 


PS I 

PHIDL 


0.0 

thetdl 


0.0 

PSIDL 


0.0 


c 

C*» BEGINNING OF OPERATE LOOP 
9000ft CONTINUE 
90002 CONTINUE 
90008 CONTINUE 
C 

C** IF FUNCTION SWITCH 1 I S ON APPLY 
C** DOUBLET INPUT TO SELECTED CONTROL 

C 

IFWnOT.LDIS! (33) > GO TO BOO 
IFMTtGT.l . )*AN0. (T.LT.2. >1 GO TO 501 
IF I (T.GT.2. 1 .AND. T.3. ) ) GO TO 50? 

A 1 SC - A1S0 
B1SC - B 1 SO 
THTRC ■ THTRO 
GO TO 500 

501 CONTINUE 

A 1 SC « AISO - A1SSTEP 
BlSC ■ B 1 SO - SISSTEP 
THTRc ■ ThTRO - THTRsTP 
GO TO 500 

502 CONTINUE 

AtSC - AISO + A 1 S STEP 
B t SC » B1S0 + B1SSTFP 
THTRC ■ THTRO 4 THTRcTP 
500 CONTINUE 
C 

C** CALCULATE CONTROLS FROM COCKPIT INPUTS 
C 

IF I .NOT.LDI S I (45 ) ) GO TO 140 
IFtlNT.LE.il CALL ADCIN 
LD I So ( 8 1 1 * • F . 

LD I SO ( 82 > - .F. 

LD I SO (03 1 » .F. 

LD 1 Sn ( 64 1 « .F. 

C 

C*« IF COCKPIT IS ACTIVATED* FUNCTION SWITCH 13 IS ON - 
C** COMPARE CONTROL STICK POSITIONS WITH TRIM POSITIONS 
C** AND TURN ON WHITE LIGHTS 21 ♦22*23.24 WHEN TRIMMED 

C 

IF t ARS (XAOS-XAOST ) , LF . XAOS*CNT£RR ) LD!S0(8l )».T. 

IF(AnS(YCS — YCS T ).LF. VCS*CNTERR 1 LO I SO ( 02 ) = • T • 

IF(ArS(XCS -XCST j.LF* XCS*CNTERR) LD I SO < 83 1 » • T • 

IF(ABS(XTR -XTRT ! «LF • XTR*CNTERR> LD t SO ( 84 1 = • T . 

GO TO 150 
C 

140 CONTINUE 

X AOS * (A0SC*DPR - AOSC0 1/XAOSG 
XCS ■ (B1SC*DPR - Bi SCO > /XCSG 
YCS ■ < A 1 SC*DPR - A i SCO ) / YC SG 

XTR * (THTRC*DPR - THTRCO l/XTRG 

150 CONTINUE 
C 

C** IF FUNCTION SWITCH 11 IS ON* SEND MODE CONTROL TO COCKPIT 

C 

1F(.nOT.LDISM 43) >G0 TO 151 
IF(INT.LE.I) CALL rPMODE 
GO TO 152 

151 CONTINUE 
mRESpT»lDT S 1(19) 

MHOLD *LD1SI(18) 

MOPER * LD 1 S lit?) 

152 CONTINUE 
C 

C 

C** IF FUNCTION SWITCH 9 I c ON* SAS SYSTEM 15 
C** ACTIVATED AND WHITE LIGHT 5 IS ON 

C 

JMES * 2 

IF(LDIS1 (41 l.OR. (LDlKl ( 44 ) . A ND«LD I S ! (55) » ) UMES = 1 

1 F ( INT.LE. 1 ) CALL cAS 
LDISOI65) = .F, 

IF (LDl S 1 <4 1 ) .OP. (LDI SI (44 ) .AND.LDIS I (55 > > 1 LDIS0(65) * .T. 

IF(LDISI f4I ).0R. (LDlSl (44 ) .AND.LDISI (55) » ) GO TO 160 
A 1 SC S ' 0. 

BlSCc ■ 0. 
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TMTPcS • 0. 

1 F C .NOT • L 01 SM 37 M GO TO 160 


AOSC 

m 

AOSO 


A! SC 

m 

A1S0 


B1 SC 


B1 SO 


THTRr 

• 

THTRO 


160 AOS * 

AOSC 


A 1 S 


A1 SC + 

A t ECS 

81 S 


B1 SC ♦ 

B1SCS 

THTR 


THTRC + 

THTRCS 

PHIL 


PHI 


THETAL 


theta 


PSIL 


PS1 


PHIDL 


PHIO 


THETDL 


thetdt 


PSIOL 


PS IDT 


At SCL 


A| SC 


B1 SCL 


BI SC 


thtrcl 


THTRC 



C 

c** RED LIGHTS 2»3»4* OR 5 ON RESPECTIVELY 
C** INDICATES AOS « Al$, flic* OR THTR IS 
C** BEYOND SET LIMITS 
C 

LD I SO ( 112) » .F. 

LD I SO (113) * .F. 

LD I SO < 114) « • F • 

LD I SO (115) » *F, 

IF(A0S .LT. *1396. OR, AOS .GT..3490) 

I F ( A | S .LT.-.1 746.0R.A1S .GT..1222) 

I F ( B 1 S *LT.-*2305.0R.B1 S .GT..2410) 
IF (ThTR.LT.-. t 222.OR.THTR.GT* .401 o ) 


L D 1 SO ( 112) ■ .T, 

LDISO t 1 l 3) = .T, 

LDI SO ( 1 l 4 > ■ .T . 
LD I SO (115) - .T. 


C 

c** call subroutines 


c 

CALL AEROVAR 

CALL FUS 

CALL TAIL 

CALL MAINROT 

CALL ENGINE 

CALL F AND M 

CALL BOOYOER 

IF fNrPR.EO. 1 ) A UTOH ■ • T. 

CALL EARTH M 

IF(lNT.LE.l) CALL nACOUT 
LDISO (116) ■ .F. 

LDTSOf 117) ■ , F * 

IF(AuTOH) GO TO 170 
IFMNT.GT.l ) 90005,180 

170 CONTINUE 

c 

C** RED LIGHT 6 ON INDICATES PHI OR THETA HAS BECOME TOO LARGE 
C** RED LIGHT 7 ON INDICATES NEGATIVE ALTITUDE 
C 

IF tNFRR.EO* 1 > LD f SO 1116) ■ • T , 

IF IH.LE.0.00 ) L D I SO (117) ■ .T. 

t80 continue 

C 

c** scaling of real time output 


OAC ( 1 ) 


XAOS*XAOSGn 

DAC < 2) 


AOS* A05GN 

DAC ( 3) 


XCS*XCSGN 

OAC ( 4) 


B 1 S*B l SGN 

OAC I 5) 


YCS*YCSGN 

OAC f 6) 


A) S*A 1 SGN 

DAC ( 7) 


XTR*XTRGN 

DAC ( 8) 


thtr*thtrgn 

DAC! 9) 


P*PGN 

DAC ! 1 0 ) 


0*0GN 

DAC ! 1 1 ) 


R*RGN 

DAC ( 1 2 ) 


PH 1 *PH I GN*DPR 

DAC ! 1 3 ) 

* 

theta*thegn*0RR 

DAC! 14 ) 


PS I *PS 1 GN*dPR 

DAC < 1 5 ) 


alff*alf«-gn*dpr 

DAC! 16) 


BETE*BETFGN*DPR 

DAC ! 1 7) 


H*HGN 

DAC ( 1 6 ) 


U* UGN 

DAC ( t9 ) 


V* VGN 

DAC ( 20 ) 


W*WGN 

OAC ( 21 ) 


A I SS*A ) SSGn 

DAC < 22 ) 


BI SS*B 1 SSGN 

DAC ( ?3 ) 


BETS 1*BETSGN 

DAC! 24 ) 


bdts 1*BDTGN 

REAL TIME 

D 

I splay 

TABLP ! 1 

) 

- A1SS 


STABLE! 1 7 ) =PH I 
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T flBL c ( 2 > = AD1SS 

T4BLPO )- BISS 
TABLF(4 )* BDiSS 
TABLF(5 >* BET S I 
TABLES >* BOTS! 
TABLE <7 )■ E 

T ABLF ( 8 )« FOOT 

TABL^(9 )- H 
T ABLF < 1 0 ) = HOOT 
TABLF (11)* N 
TABLF (12)= NDOT 
T ABLE (13)= OMFG 
T A0L F ( 1 4 ) * OHFGOT 
TABLE (15)= P 
T ABLF (16)* POOT 
TABLF ( 33 ) = AOS 
TABLF ( 35 1 =B1 S 
TABLF ( 37) =OMEG0 
TABLF ( I 25 )*CT 
TABLF ( 1 27 ) = CO 
T ABLF ( 1 29 ) = HT I P 
TABLF < 1 31 ) = SHP 
T ABL f ( 133) =YCS 
TABLp (135) *XTR 


*TABLF(1B)=d h IO 
* T ABL F ( t 9)*PSl 
*TABLE(20)=PS10T 
STABLE (2 1 ) *0 
*T ABLE (22 ) *000 T 
STABLE (23 ) *9 
ST ABLE (2 A ) = RDOT 
ST ABLE (25 ) =THE T A 
STABLE (26) *THETDT 
STABLE (27)»U 
ST ABLE (20 ) =UOOT 
S T ABL E (29 >=V 
STABLE (30) = VDOT 
STABLE (31 ) = W 
STABLE! 32 ) = WOOT 
STABLE (34 ) *A [ S 
ST ABLE (36)*THTR 
STABLE! 38 ) *T 
STABLE ( 1 26 )*CP 
STABLE! 1 28 )*MU 
ST ABLE ( 1 30 >»ASPD 
STABLE! 1 32 )*XAOS 
STABLE! 1 34 )=XCS 


IF (LOI S 1 (22 ) > CALL SCANNER (KSC An ) 

CALL DSPLAY 

IF (Lnr SI ( 17 ) ) GO To 200 
IF(VARCHNG) call typfvar 
I F (ENABLE.AND.TYPESW ) CALL TYPEVAR 

90050 CONTINUE 

ENABLE ■ .NOT. TYPES* 

IF ( AuTOTYP.AND.LO! S 1 ( 1 4 ) ) CALL TYPEVAR 
200 CONTINUE 

CALL RTMODE 
90001 CONTINUE 

C 

C** RFC ORO ON RFAL T J he filf 

C 

CALL RFC ORO 

C 

c 

TF(AuTOH) CALL RECYCLE 

r 

c## TIC HARKS RECORDED ON CHART OUR 1 NG EACH SECOND OF OPERATION 

C 

LDI SO ( 1 02 )= .F. 

L D I S n ( 1 1 0 >».F. 

IF(AmCD(T.1 • ) *NE.O, ) GO TO 230 
LD1S0C I 02 )*.T. 

L D t SO ( 1 l 0 ) * * T • 

23^ CONTINUE 
C 

90005 CALL IGPATEl 

C 

C** IF INTSCHE « 5 . USE AoAMS-B A SHFORTH INTEGRATION 
IF ( InTSCWE.EG.5 ) InT - 1 
IF (H,LE«0» 000 ) AUTOH - .T. 

PSI * AMODfPS! .6.2831 8 ) 


INTEGRATION FLOW CONTROL 

C 

IFUnT.LE.I) CALL RECYCLE 
GO TO 90008 
C 

90014 CONTINUE 

WR1 Tp (MF I .240 > I RUN 

IF< IDRINT.GT.I) GO TO 260 

CALL DAYTIM(TIME) 

WRITf ( MF 1 .250 ) I RUN. T I HE (1 ) . T t HE ( 2 ) 
WR I Tf ( HF I « 255 ) MESAGF(JHES) 

IPRINT * 2 
260 IRON ■ 1 RUN + 1 

C 

90030 CALL PL AY0AK ( 90032S ) 


WRITf ( HF I . 300 ) 


X 

T 

• A SPD 


GRWT 


CG 

X 

H 

. HOOT 


OHEG 


omegdt 

X 

U 

• V 


W 


ALFF 

X 

UDOT 

. VOOT 


WDOT 


BETF 

X 

P 

• 0 


R 


DELE 

X 

POOT 

. QDOT 


RDOT 


LHR 

X 

XA 

. YA 


ZA 


OHR 

X 

LA 

♦ HA 


NA 


FXR 

X 

PHI 

. THETA 


PSI 


FYR 

X 

PH ID 

.thctot 


PSIOT 


QE 

X 

AOSS 

. A I SS 


BISS 


WIH 

X 

AOS 

• A 1 S 


B 1 S 


THTR 
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r 

C ** FUNCTION SWITCH 6 ON SELECTS ALTERNATE PRINTED OUTPUT 

C 

IFILOISI (381 ) GO To 90029 
WR I Tc- ( MF I *30j ) 

X XAOS • YCS . XCS 4 XTR . 

X N . NDOT , E 4 FDOT , 

X OMAX t SHP t CT 4 VT » 

X QV 4 OL 

GO TO 90030 

r 

90029 CONTINUE 

WR1TF l WF! ,302 I 

X XF , YF , ZF 4 VTR . 

X LF , WF , NF 4 YTR . 

X IKH«; , Y VS 4 ZHS 4 WIWM . 

X XW , ZW 

GO TO 90030 
POO 32 CONTINUE 
C 

C** FORMAT STATEMENTS 

c 

2A0 FORMATM 2SXI6 ) 

250 FORMAT (4X8HRUN NO * I6.8X12HTHE DATE IS A1048X12H THE TIME IS *10/ 
1 ) 

255 F ORM a T ( A X 7HS AS IS A J n/ ) 

300 FORMAT ( 1 H21 0X1 HT , E20 . 8 4 8X AHA SPO « E 20 • 8 4 0X4HGRWT .020.8. 1 0X2HCG . £20 . 8 
X/ 1 l XI HH4E2O.840X4HHDOT ,£20 . 8 . 8X4 HOMEG 4 E20 . 8 4 6X6H0MEGDT . E20 . 8/ 

X 1 1 XIHU,E20 .04 I | XI hV 4 E20 484 1 1 XI HW 4 020 ♦ 8 4 8X4HALFF 4020.8/ 

X SxAhUDOT 4 E20 .8 .8X4HVOOT.E20.6 4 0X4HWDOT , £2 0 . 8 4 8X4HBETF t E20 • 8/ 

X 1 1 XI HP. E20.e . 1 ! XlH0 4E2r-8 4 1 1 XI HR, E20. 8, 8X4HDELE ,E20.8/ 

X 8xAHPD0T,E20#8 48 X 4 H 000 T 4 E 2 C. 8 .8X4HRD0T ,E 20 . 849 X 3 HLMR, E 2048 / 

X 1 0X2HXA 4020.8, 1 0X2HV A 4 020 . 6 , 1 0X2HZA 4 E 20 4 8 4 9X3H0MR . E20 4 8 / 

X 1 0X2 HL A 4 020.8, 1 0X2HM A 4020.04 1 0X2HNA 4 E20 . 8 4 9X3HFXR .E20 . 8/ 

X 9X3HPHI 4E2O.8,TX5hTHETA4E2O.0 4 9X3HPS I 4 £20 . 8 4 9X 3 HF YR 4 E2 0 . 8/ 

X 8x4 HPH 1 D 4 E 20 #6 . ftX6H T HETDT .E20 .8. 7X5HPS I DTi E20 * 8 4 1 0X2HGE 4E20. 8/ 

X 8x4 HAOSS 4 020.6 4 B X4 H A 1 SS4E2C.8 ,0X4H6! SS . 02 0 * 8 . 9 X 3HW I M . 0 2 0 . 8/ 

X 9x3HAOS 4 020.8* 9X3hA ! S 4020 .8 4 9X3HB 1 S 4E20.8 4 0 X 4 H T H TR 4 E 2 0 . 8 ) 

30 1 FORMA T ( 1 H2 

x 7 x 4 HXAOS 4 E 2 C *8 4 9X3H YCS 4 E20 . 8 4 9X3HXCS 4 E20 . 8 4 9X 3HXTR , E20 . 8/ 

X 1 I Xl HN«E20 .8 *0 x4HnDOT 4E20 .8 . 1 1 XI HE 4E20#8 . 8X 4H0O OT , E20 • 8/ 

X 8x4HQm AX 4E20.8.9X3HSHP. 020*8. I 0X2HCT 4 E20 4 8 4 l 0X2 HVT , 020 . 8/ 

X 1 OX2HQV1 E20. 8. 1 0X?HQL 4020.8// } 

30? FORMAT (1H2 

X 09X2 HXF 4 E2 0.8 4 1 0X2HYF .020.04 10X2HZF . 02 0 . 8 4 9X3HV TR 4 £20 . 8/ 

X I OX2HLF 4 020.8. 1 OX2HMF .020.8. 1 0X2HNF ♦ 0 2 0 • 6 . 9X3H Y TR , 020 . 8/ 

X 9x3HWHS4E2O.B,9X3hYVS4E2O.8.9X3HZHS4 02O, 8, 8X4 HU I WM » E 20 .8/ 

X 1 OX2HXW. E20.8. 1 0X?MZW . 020 .6// ) 

C 

CALL APRINT 
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ADCIN 


Subroutine ADCIN contains the equations for the helicopter collective and cyclic 
stick inputs and the tail rotor pedal inputs (ref. 1 and appendix A). The subroutine allows 
two modes of operation. Primary control is from the control sticks and pedals in the 
cockpit, and alternate control is from the program control station (fig. 1(a)) through the 
use of potentiometers (located on the program control console (fig. 1(b)) representing the 
control sticks and pedals. The alternate mode of operation allows the analyst to conduct 
preselected flight operations to be used, for example, for comparison with actual flight 
data. 
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SUBROUTINE. ADC I N 


C. 

c 

c 

c 

c 

c 


c 

c 

c 


c 

c 


c 

c 


c 


r 


c 


PROGRAM PROVIDES ADC INPUTS THROUGH 
POTENTIOMETERS OR COC<PIT 

TYPE STATEMENTS AND DIMENSION :'. yAOI ALLI S 


RPAL LIM 
LOGICAL LOIS I 


INTEGRATION COMMUNICATION 


COMMON 


X 

/ I ntcomm/ 

T 

DT 

I NT 

. NEO 

. 


X 



I SCHEME i 

DE RINT <2 

.12 1 




X 

/REAL TIM/ 

ADC ( 3? ) « 

DAC ( 64 ) ♦ 

L D I S 1 

{ 1 OH 1 . LD I 

SO < 1 76 ) . 


X 



NO PER 

NH OLD . 

NWfcT.E 

T , N T E 

PM , 


X 



npr I n t * 

NRt Au 








SI. J 

IHROUT 1 NE 

COMMON 1 CAT I ON 



COMMON 









X 

/cntl c 

/ 

AO SC , 

A 1 SC . 

Bt SC 

. THTRC . 

AlSO. . 

H 1 ^Cl 

X 



THTRCL 






X 

/cntl p 

/ 

AO SCO 

A I SCO . 

R 1 SCO 

, T H TRC O . 

XA'V-vP- . 

Y •• . NG 

X 



XCGG 

X T RG . 

xAOSW 

. VC SR 

YCSR . 

X TRP 

X 

/L AGSTK 

/ 

NILAC- 

NE T ICK 





X 

/SCALES 

/ 

SF 1 

SF2 . 

SF3 

. 3^4 

: r ^ 1 ■ 

: • 6 

X 



S E 7 






X 

/ST I CKS 

/ 

X A os . 

Y C S . 

xc s 

. x r p . 

DPR 

"•“D 

X 

/ST I NP JT/ 

X 1 « 

X2 

X 3 

« X4 , 

T 1 

TO 


1 nn 


101 


L|V(X*Aifl) = AM I N 1 ( R . AMAX 1 ( X . 
I F ( LD I s I ( 44 M GO TO 1 n;> 

I F ( N S T I C< • NF • 1 ) GO TO I 0 0 


J F ( T.LT . T 1 

) (iO TO 1 

OO 

IF(T«C»T«TP> r » r TO I 

~,r 

X 1 1 

X 1 


X22 

X? 


X 3 3 = 

X 3 


X44 = 

X 4 


GO TO 101 



CONT I NUF 



XI ! 

0 . 


X22 

0. 


X 3 3 

o . 


X 4 4 = 

0. 




STICK 

C ONT I NUF 



X AO S = 

ADC ( 1 ) *xAOSP 

+ X 1 1 

XCS = 

ADC < 3 ) *XC SR 

+ X 3 3 

YCS = 

ADC ( 2 > * YCSR 

+ X 22 

XTR 

A DC ( 4 ) * X TQR 

+ X 44 

GO TO 103 




A ) ) 


DEFLECT [ ONS 


THROUGH pr.Tf NT 


r,W Trwf, 


COCKPIT STICK DEFLECTIONS 


102 CONTINUE 


XAOSCOL 

= 

L I M ( SF 1 ♦ADC < 1 ) 


« 

0*0, 

XCSCOL 

= 

L I M ( SF3 * ADC < 3 ) 

+ 

SF 6* ADC ( 6 ) * 

“4.8* 

YCSCOL 


L I M ( SF2*AD C (2 1 

+ 

SF5*ADC ( S ) ♦ 

-4 .8 . 

XTRCOL 

= 

L IM ( SF4*ADC (4 ) 

+ 

S p 7 * A D C ( 7 ) . 

-3 • 0 . 

X A O S 

= 

XAOSCOL 




XCS 

- 

XCSCOL +• 4.B 




YCS 


YCSCOL 4 4,8 




XTR 

= 

XTRCOL + 3.C: 





100 CONTINUE 


AOSC 

= (AOSCC 

+ 

XAOSG*X AOS 

) »RPD 

A 1 SC 

= ( A 1 SCO 

+ 

v CSG 

♦ YCS 

) *RPD 

B 1 SC 

= ( B 1 SCO 

+ 

XC E G 

♦ xcs 

>*RPQ 

THTRC 

= (THTRC0+ 

XTRG 

♦ XTR 

) *RPD 

RETURN 






END 







13.3 ) 

4 ♦ B 1 

4.B 1 

3 . C ) 
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AEROVAR 


Subroutine AEROVAR computes the main rotor inflow velocity, the lateral and ver- 
tical dynamic pressures, and the fuselage airspeed (refs. 1 and 2 and appendix A). 
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subroutine aerovar 


real lmr 

LOGICAL 


LDI SI 


COMMON 


/ 1 NT COMM / T 


/REAlT im/ 


PROGRAM COMPUTES ROTOR INFLOW VELOCITY, 
FUSELAGE AIRSPEED* AND DYNAMIC PRESSURES 

TYPE STATEMENTS AND DIMENSIONED v A P I ARLFS 


INTEGRATION COMMUNICATION 


♦ DT , INT « 
ISCHEME, DFR INT ( 2 * 12) 

ADC ( 32 ) ♦ 0AC(64), LDISIC108)* 
NOPER . NHOLD * NRESET * 
NPR1NT , NPEAD 


LDI SOT 196 ) . 
NTERM , 


SUBROUT I NE COMMUNICATION 


COMMON 


X 

/ADV C / 

OMEG 

* W I M , 

VT 

♦ 

QV 

* QL 

* ASPD 

X 

/ADV P / 

ADVP 1 

, A0VP2 






X 

/EUL CS/ 

COSE I 

, S I NF I , 

COSTH 

♦ 

S INTH 

. COSSI 

• SI NS I 

X 

/FTRmPLS/ 

EQLMR 

, EQAOSS. 

EQWI M 

% 

HDOTD 



X 

ZGPARAM/ 

DELT 

* PHO * 

G 

t 

PIE 



X 

/MROT C/ 

FXR 

* FYR * 

lmr 

t 

LRH 

, MRH 

, QMR 

X 


AOSS 

, ALESlY* 

UTS 1 v 

% 

bets r 

* BDTSI 

* Al SS 

X 


ADI SS 

, BISS * 

BD1SS 






EQUIVALENCES 


EQUIVALENCE 


X 


( DER INT ( 1 , 

4 ) , PH 1 ) * 

X 


( DER INT ( 1 , 

5 ) , THETA ) , 

X 


{ DER INT ( 1 , 

6 > *PSI ) * 

X 


( DER INT { 1 * 

7 > *U ) * 

X 


( DER I NT ( ! * 

8 ) * V > » 

X 


C DER I NT ( 1 , 

9 ) * w ) 




COMPUTE 

COSF I 

_ 

COS (PHI) 


S I NF I 

5 

S I N { PH I ) 


COSTH 

= 

COS ( THETA ) 


S I NTH 

= 

SI N { THETA ) 


COSS I 

= 

COS (PS 1 ) 


SINS I 

= 

S I N (PS I ) 



SINES and COSINES OF FuLER ANGLES 


COMPUTE LINEAR VELOCITY ALONG 7-Ay IS 
AND NET AIR VELOCITY PAST Hi iB 

IFT.NOT. LD ISM 37 ) ) GO TO l 00 

W = ( U*S I NTH - V*S I NE I *COSTH - HDOTD >/ (COSTH*COSF I ) 

100 CONTINUE 

WM = W - WIM 


COMPUTE AIRSPEED* TOTAL VELOCITY THROUGH 
MAIN ROTOR*AND MPAN INFLOW VELOCITY 

ASPD = SORT ( U*U + V*V + W*W)/ 1,609 

VT = SORT ( U*U + V*V + WM*WM) 

I F ( V T • LE • 0 • ) VT s ,00001 

WIMP = ADVP 1 ♦L.MR/RHO/VT 

IF(INT,LE*1> WIM = ADVP2*TWIMP - WIM) *DELT + WIM 

LATERAL and vertical dynamic pressures 

QL * 0,5 »RHO* { U*U + V*V> 

QV = 0.5*RH0*(U*U + WM*WM) 


RETURN 

ENO 
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BODYDER 


Subroutine BODYDER contains the equations necessary to compute the total body 
angular and linear accelerations (refs. 1 and 2 and appendix A). The program also con- 
tains a safety feature to prevent program abortion if $ or 9 becomes too large. 
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ooo oooooooo non nr»r> ooo a no 


SUBROUTINE BODYDFR 
C 

C PROGRAM COMPUTES body motion ANGULAR 

C ANO LINEAR ACCELERATIONS AND VELOCITIES 

C 
C 

C TYPE STATEMENTS AND DIMENSIONED VARIABLES 

C 

REAL LA* MA, NA* TXX, I YY * IZZ* MASS 

integration communication 

common 

X /INTCOMM/ T . DT * INT * NEG * 

X ISChEME* DER I NT (2*12) 


SUBROUTINE COMMUNICATION 


COMMON 


X 

/ERROR / 

nerr 






X 

/EUL CS/ 

COS*! 

* SINFI 

* COSTH 

. SINTH 

t COSSI 

. SI NS I 

X 

/F N M C/ 

XA 

* YA 

* ZA 

. LA 

* MA 

, NA 

X 

/GPARAM/ 

DELT 

* RHO 

* G 

* PIE 



X 

/SHIP P/ 

IXX 

* I YY 

♦ IZZ 

* MASS 




EQUI VALENCES 


equivalence 


X 

( DERINT ( 1 f 

1 ) *P 

) * 

X 

( DER I NT ( 1 * 

2 ) .0 

) * 

X 

(DERINT ( 1 * 

3 ) *R 

> * 

X 

(DERINT ( 1 , 

4 ) .PHI 

) • 

X 

(DERINT ( 1 * 

5) • THETA ) * 

X 

(DERINT ( 1 , 

6 ) *PSI 

) * 

X 

(DERINT ( 1 , 

7 ) *U 

) • 

X 

(DERINT ( 1 * 

8 ) • V 

) * 

X 

(DERINT f 1 , 

9 ) . W 

) * 


(DER I NT (2* 
(DERI NT (2* 
(DERI NT (2* 
( DER I NT ( 2 * 
( DER I NT ( 2 * 
(DER I NT (2 * 
(DERI NT <2, 
( DER I NT ( 2 * 
(DER I NT (2* 


1 ) *POOT 

2 ) *ODOT 

3 > * Root 

4 ) .PHI D 

5 ) * THETDT 

6 ) « RS I DT 

7 ) .UDOT 

8 ) * VDOT 

9 ) * WDOT 


COMPUTE ROLL* PITCH* AND YAW ACCELERATIONS 


MOOT 

QOOT 

ROOT 


(LA ♦ ( IYY-I7Z)*Q*R)/IXX 
(MA + ( IZZ-IXX)*R*P)/IYY 
(NA ♦ ( IXX-I YY)4P*Q)/IZZ 


IF PHI OR THETA GETS TOO LARGE PROGRAM 
PUT INTO AUTOMATIC HOLO TO PREVENT ABORTION 
WITH A MODE 4 

IF( ARS(THETA ) *GT* 1 .5.0R*ABS (PHI ) *GT* l #5 ) GO TO 10 

COMPUTE EULER ANGLE DERIVATIVES 

PSIDT * CR*COSFI + Q*SINFI )/C0STH 
THETOT* 0*C0SFI - R*SINFI 
PHID - P ♦ P5IDT*SINTH 

COMPUTE BODY VELOCITY DERIVATIVES 

UDOT * XA/MASS + R*V - Q*W - G*S I NTH 

VDOT ■ YA/MASS + P*w - R*U + G*COSTH*S I NF ! 

WDOT « ZA/MASS + Q*U * P*V + G*COSTH*COSF I 

C 

RETURN 
10 CONTINUE 
NERR • 1 
RETURN 
END 
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CPMODE 


Subroutine CPMODE contains the logic statements necessary to allow mode control 
of the computer program from the cockpit by the pilot. 
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no nn no noon 


SUBROUTINE CPMODE 


C 

C 

c 

c 

c 

c 

c 

c 


PROGRAM PROVIDES MODE CONTROL TO COCKPIT 
TYPE STATEMENTS AND DIMENSIONED VARIABLES 
LOGICAL LDISI. mreset. MHOLD* MOPER 

INTEGRATION COMMUNICATION 


COMMON 

X /REALTIM/ ADC ( 32 ) * DAC(64), LDISI(108>« LDIS0U96). 

X NOPEP • NHOLD , NRESET . NTERM • 

X NPRINT ♦ NREAD 


SUBROUTINE COMMUNI CATION 


COMMON 

X /MODEC / MRESET . MHOLD • MOPER 


I F ( LD I S I (49) ) GO TO 100 
IF (LOIS 1(50)) GO TO 200 
IF ( LO I S I (51 ) ) GO TO 300 
lF(MQESET) GO TO 1 00 

IF(MMOLD) GO TO 2n0 

IF(MOPER) GO TO 3n0 

GO To 400 

RESFT MODE 

100 MRESfT = * T « 

MOPER = • F • 

MHOLD = • F • 

LDISI (1 9) = • T « 

LDISI (17) « #F # 

LDISI (18) * • F * 

GO TO 400 


HOLD MODE 

200 MHOLn = *T# 

MOPER « • F • 

MRESET * • F # 

LDISI ( 1 8 ) = • T « 

LD I S I ( 1 7 ) = .F* 

LDISI (19) = • F * 

GO To 400 


OPFRATE MODE 

300 MOPER = • T • 

MHOLD = • F • 

MRESFT = # F # 

LDISI ( 1 7 ) = • T # 

LDISI (18) = • F * 

LDISI (19) * • F * 

C 

400 CONTINUE 
C 
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RETURN 

END 



D AC OUT 


Subroutine DACOUT contains the equations and logic necessary to drive the cockpit 
instrumentation. The following instruments are included: 

(1) attitude direction indicator 

(2) altimeter 

(3) rate of climb indicator 

(4) two minute turn indicator 

(5) bank angle indicator 

(6) main rotor RPM indicator 

(7) airspeed indicator 
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SUBPOUT INF DACOUT 


PROGRAM provides dacs needfo to 

DRIVE COCKPIT INSTRUMENTS 


TVPF STATEMENTS AND DIMFNSIONFn VARIABLES 


REAL LIM 

LOGICAL LDISI* LDISO, MRESET. mholo. MOPFR 


INTEGRATION COMMUNICATION 


/InTCOMM/ t 

ISCMEME 
/REALTIM/ ADC I 32 ) 
nopep 
nprjnt 


neo 


OT * I NT 

DERI NT ( 2 * t 2 ) 

D A C ( 6A ^ « LDISl ( 108) i LDIS0M96). 
NMOLO « nrfept ■" r '“ 

nread 


NT E RM 


SUBROUTINE COMMUNICATION 


/ACCEL/ 
/AD V C 
/CPI C / 


/GPAPAM/ 
/LAGSTK , 
/MOOEC 
/STICKS x 


QMFG 4 

AC I 4 

YPTMI , 
EOS (IB) 
DELT , 

NILAG i 


ypBMl , 
RHC 4 

mholD 4 

VCS 


A? 

VT 
DC I 
vpgw i 


MOPFR 

yes 


ov 

ypCMi 

CSKA 

P l F 


OL 

VRCM| 

CSK6 


aepd 

yDTM| 


EOU I VALENCE 


EOUI VALENC 

(DERI NT r 1 « A ) »PH] ) • 
(DFRINTM* S).TMFTA), 

( DPR I NT ( 1 ♦ 6 ) 4 PE I ) 4 

( DER I NT { 1 , 3 1 4 P ) 4 

( Deq INT M * 1 2 1 .H ) 4 

( DEP 1 NT ( 2 4 1 2 ) . HOOT > 


DATA PALT / 0.0 / 

L I M < X 4 A 4 B ) * AMlNlfB.AM 

PHJOAL i-AY/A? 

OMEGRPM = 0MEG#DPR/6 .0 


instrument TEET 


I F { .NOT . 
PH I 

theta 
PS I 
HDOT 
H 

PHIBAL 

R 

OMFGRPM 

ASPO 

10 c continue 


LDI SI ( 42 ) ) 

' 0.524 
= 0.7854 
= 0.524 
= 25 . 

= 50. 

= 0. I 

= 0.0262*2.0 

= 324 . 

' 130. 


GO TO 1O0 


DAC (25 ) 
DAC (26 ) 
DAC (27 ) 
DAC (28 ) 
OAC (29 ) 
DAC (30 > 


palt 

AAA 

DAC ( 31 ) 
OAC < 32 ) 


-COS (PHI) 

-S I N ( PH I ) 
-COS ( THETA ) 
-S I N (THETA > 
-COS(PS I ) 

— S I N ( PS 1 ) 


ALT I METER 


PALT + L I M ( H - d A (_ T 4 -15.4 l ^ • ) 

PIE - AMOD (PALT » t 000. ) *CSK4 
COS < AA4 ) 

S I N ( A Aa } 


= 1 .0 


ZSOP 
1 F ( Z SOP ) 

200 CSK5 
GO TO 2C2 

201 CSK5 

202 CONTINUE 

I F < 7 SOP . GE *1000 
XX = 100. 

IDx =0 

GO TO 301 

300 XX * 500. 


L I M (HOOT *60 . . 

200. 201 . 201 

-1 *C 


Rate OF CLIMB 
-3999.994 3999.99) 


• or. zsdp.le. - icon. ) 


GO T Ei 30 n 
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non non onn 


iDy * 8 

301 SSI = ZSDP/XX 

I SI * SSI 

IS2 * IDX*CSK5 

IS * IABSdSl + IS2) + 1 

SS * ABS (SSI - IS 1 ) 

ZSDF = CSK5*((1*0 - SS)*F05(IS) + SS*F05(IS + 1)) 

IF (NILAG.EQ. 1 ) 302, 303 

302 AA5X * ZSDF*RPD 

AA5 * DCI*YRCMt + BC I *XRCM 1 + ACI*AA5X 

XRCM1 * AA5X 
YRCM1 = AA5 
GO TO 304 

303 CONTINUE 

AA5 * ZSDF*RPD 

304 CONTINUE 

DAC(33) * -COS ( A A5 ) 

DAC<34)= SIN(AA5) 

TWO MINUTE TURN 

IF (NILAG.EQ. 1 ) 400, 401 

400 AA6X c LIM((R*CSK6 + PIE), 0.0, 6,2831) 

AA6 * DC I * YRTM 1 + BCI*XRTM1 + ACI*AA6X 

XRTM1 = A A6X 

YRTM I a A A6 
GO TO 402 

401 AA6 * LIM<(R*CSK6 + PIE), 0.0, 6*2831) 

402 CONTINUE 

DAC(35) * LI M (COS ( AA6 ) , -.806, .886) 

DAC ( 36 ) * L I M ( S IN ( AA6 ) , -.500, .500) 

BANK 

IF (NILAG.EQ. 1 } 500, 501 

500 A A 7x “ (LIM(PHIBAl, -.15, .15>)*CSK7 

A A 7 = DC I * YPBM 1 + BCI*XPBM1 + ACI*AA7X 

XPBM1 = A A 7X 
YPBM1 = A A 7 
GO TO 502 

501 A A7 = (LIM(PHIBAl, -.15* . 1 5 ) ) *CSK7 

502 CONTINUE 

DAC ( 37 ) = —COS ( A A7 ) 

DAC ( 38 ) = -S I N ( A A 7 ) 


C 

C 

C 

C 

C 

C 


C 


RPM 

DAC<39> = L I M ( OMEGRPM/ 500., 0.0, .9996) 

indicated AIRSPEED 
DAC (40 )*L IM ( ASPD/1 4 0 . ,0.0,0.9996 ) 

MODE CONTROL DISCRETES 

LDISO ( 1 ) a MRESET 
LD I SO ( 2 ) = MHOLD 

LDISOO) * MOPER 

RETURN 

END 
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EARTH M 


Subroutine EARTH M contains the equations used to calculate the earth rates by 
resolution of body linear velocities through the Euler angles (refs. 1 and 2 and 
appendix A) . 


EARTH M 


Compute Earth Reference 
Velocities and Earth Rates 



SUBROUTINE EARTH M 


PROGRAM COMPUTE FAOTH PAT^S QF^OLUTION 

OF BODY AXIS VELOCITIES THROUGH FULFR ANGLES 


TYPE STATEMENTS AND O I MENS I ONE n VARIABLES 


REAL MOOT 


integration communication 




INCHEMF, OF R 1 NT (?i IP ) 

Si 10 ROUT I NF COMMON I rA' 


COMMON 

X /EUL CS/ COST I 


EQUIVALENCE 


♦ SINFI * CCSTH 
EQU I VALENCES 


CO c ’SI « S I N c I 


X 

( DERI NT ( 1 « 

7 ) , U 

) * 

X 

( DER I NT ( 1 , 

H ) , V 

) t 

X 

( DFR I NT { 1 , 

0 ) , w 

) 4 

X 

( OFR I NT ( E « 

1 0 ) « N D O T 

) « 

X 

( OER I N T (?i 

1 1 ) ,FL)OT 

} % 

X 

( DFR I MT (2, 

1 2 ) ♦ HDO T 

) 


COMPUTE EARTH REFT PF NC:- VELOCITIES 


VP^ 

= V*COSFl - 

X*S I NF I 

WPF 

= V* SINFI + 

W*COSF I 

UPF 

= U*COSTH + 

>;pf. *S 1 NTH 


COMPUTE EARTH RAT 1 -: 


HOOT = U*SINTH - b'PE*COSTH 
NEtOT = UPF* CO SSI - VPE*SINS! 
EDOT = UPE*SINSI + VPE*COSSI 


RFTURN 
f nh 
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ENGINE 


Subroutine ENGINE contains the equations used to calculate required engine torque, 
maximum engine torque available, and the main rotor angular velocity (refs. 1 and 8 and 
appendix A). 
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subroutine ENGINE 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


c 

c 


c 

c 

c 


c 

c 

c 


PROGRAM COMPUTES REQUIRED ENGINE TORQUE* 
MAXIMUM ENGINE TORQUE AVAILABLE* AND 
MAIN ROTOR ANGULAR VELOCITY 


TYPE S t ATEMFNTS AND DIMENSIONED VARIABLE'S 


REAL IROT. MU* mt IP 
LOGICAL LDI S I * LDI SO 

INTEGRATION COMMUNICATION 

COMMON 


X 

/INTCOMM/ 

T ♦ 

DT * 

INT 

* NEO 



X 


1 SCHEME , 

DERINT (2 

.12 1 




X 

/REALT 1 M/ 

ADC (321. 

DA C ( 64 J * 

LO 1 S I ( 1 06 1 « LDI SCH 1961. 


X 


NORER * 

NHOLO , 

nreset 

, NTFRM , 


X 


NPQINT , 

nread 







subroutine 

COMMUN ICAT I ON 



COMMON 








X 

/ADV C / 

OMEG 

W 1 M 

VT . 

QV * 

QL 

ASDD 

X 

/COEF / 

CT 

CP * 

CQ . 

MU » 

MT IP • 

VS 

X 

/CONTRL/ 

AOS 

A 1 S * 

B 1 S , 

THTR . 

omegd 


X 

/ENG C / 

QM A X 

OE 

SHP • 

GRWT , 

OMFGDT 


X 

/ENG P / 

IROT 

CKE ! 





X 

/GpARAM/ 

DELT 

RHO • 

G . 

PIE 



X 

/MROT C/ 

FXR 

, FvR , 

LMR * 

LRH , 

MR H , 

QMO 

X 


AOSS 

ALESI y. 

UTSIY , 

BETS! . 

BDTSI . 

A 1 SS 

X 


ADI SS 

BISS 

001 SS 




X 

/ROT P / 

TWIST 

F MRS « 

CX AO * 

CKABI , 

CYAB? , 

CKAF3 

X 


NP A p 

NA 2 » 

NS . 

PB 

ffh 


X 

/SHIP P/ 

1 XX 

Ivy * 

IZZ , 

MASS 




EQU I VALENCES 


equivalence 

X (DER INT { ! * 7 ) ,u > ♦ 

X f OER I NT ( 1 « 1 2 > * w 1 

compute ENGINE TORQUE REQUIRED 
AND LIMIT TO maximum 


OEO = CKEI * (OMEGO - OMFQ ) + OMR 

IFI.NOT. LDI5K37I1 GO TO 100 
OMEG = OMEGO 
OE = QMR 

100 CONTINUE 

IFMNT.LE.I) QE s OE + (OED - QE J*DELT 
OMAx = ( 202 .631 50* I H - 250C • 1 * 412500*1/ OMEG 

IF! h.LT* 2500 • > QM A X = 605000./ OMEG 

I E (LO l S I C 37 | > GO TO 200 
I F ( OE.GT.OmAx 1 OE = OMAX 

200 CONTINUE 
C 

c RED LIGHT 0 ON INDICATES 

C ENGINE TORQUE GREATER THAN MAXIMUM 

C 

LD I SO ( 1 1 6 ) = .E. 

I F f QE • GE • QM A % 1 LDISOI110J = *T. 

C 

C COMPUTE ROTOR ANGULAR ACCELERATION 

C 

OMEGDT= (OE - QMR )/ IROT 

IFIINT.LEiI I OMEG = OMFG OMEGDT*DELT 


c 

c 

c 

c 

c 

c 

c 


c 

c 

c 


COMPUTE SHAFT HORSEPOWER 


SHP a QMR* OMEG/ 550 • 


COMPUTE ENGINE THRUST, POWER, AND TORQUE 
COEFFICIENTS 


CT 

CP 

CO 


MU 

MT IP 


GRWT / (RHO*PIE*R0*RS* ( ( OMEG *R0 1**2 » ) 

(550.*SHP )/ <RHO*OfE*RB*RB* ( ( OMEG*RB 1**3) ) 

QMR/ I RHO*P I E *R9*R3*R9* ( ( OMFG*RB ) **? 1 1 

COMPUTE TIP SPFEO RATIO ANO MACH NUMBER 
( I • 6 09* A SPO ) / ( 0MEG*R9 1 
<I*609*ASPD + OMEG*RB ) / VS 


RETURN 

END 
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EVAL 


Subroutine EVAL is used by the trim circuit algorithm and sets up an (n + 1) x (n + 1) 
matrix of functional evaluations and computes (n + 1) sums where each sum is the sum of 
the trim equations squared for each column (ref. 6 and appendix D). 


(Tvj 

D 

Set Up (n + 1) x (n + 1) Matrix 
of Functional Evaluations 


L 

Square the Trim Equations 
and Compute Summation 
of Each Column 

("ret 

urnT) 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

r 


c 


c 


c 


SUBROUTINE EVAL(F.NCOL) 


PROGRAM SETS UP AN (N+1JXCN+1) MATRIX OF 
FUNCTIONAL EVALUATIONS AND COMPUTES (N+l ) 
SUMS WHERE EACH SUM IS THE SUM OF THE TR I 
EQUATIONS SQUARED FOR EACH COLUMN 


M 


TYPE STATEMENTS AND DIMENSIONED VARIABLES 


DIMENSION F ( t 1 ) 


SUBROUTINE COMMUNICATION 


COMMON 

X /TRIMAT/ 

X 


XMAT ( 1 1 • 1 2 ) * 
NEQPl « 


EM AT ( I 2 1 1 3 > • TRSUM < 1 2 ) • 
NEQP? ♦ XMATG< 1 1 ) 


NEON « 


EMAT( 1 * NC OL ) * 1 

TRSUM (NCOL) * 


DO 100 J = 2 * NFQP 1 
EMAT ( Jt NCOL ) = F ( J- l ) 

TRSUM (NCOL) = TRSUM (NCOL) 

too Continue 


+ EMAT ( J.NCOL ) **2 


RETURN 

END 
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F AND M 


Subroutine F AND M contains the equations used to calculate the total aerodynamic 


forces and moments (refs. 1 and 2 and appendix A). 


( F AND M ) 



r 

Sum Forces, Compute Moments 
About Center of Gravity, and 
Compute Accelerations 





(RETURN ) 
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non non non non 


SUBROUTINE F AND M 


C 

C 

C 

C 

C 

C 

c 

c 

c 

c 


PROGRAM COMPUTES NET FORCES And MOMENTS 
ACTING ON HELICOPTER 

TYPE STATEMENTS AND DIMENSIONED V A R I ABLF S 
REAL I R • LMR • l.RH, MRH, LA» MA f NA , LF « MF , NF « MASS 

integration communication 

COMMON 

X /INTCOMM/ T ♦ DT , INT 4 NEQ 4 

X I SCHEME # DER I NT (2*12) 

SUBROUTINE COMMUNICATION 


COMMON 


X 

/ACCEL/ 

AX 

t 

ay 

. AZ 







X 

/adv C / 

OMEG 

4 

W I M 

• VT 

« 

ov 

4 

OL 

4 

ASPD 

X 

/ENG C / 

QM A X 

* 

QE 

4 SHP 

« 

GRWT 

4 

OMEGDT 



X 

/F N M P/ 

DX 

• 

DXHS 

« DXTR 

4 

DXVS 

4 

DZ 

4 


X 


DZTR 

• 

DZ VS 

4 DXW 

4 

DZW 

4 

IR 



X 

/F N M C/ 

YA 

% 

YA 

4 ZA 

4 

LA 

4 

MA 

4 

NA 

X 

/FUS C / 

XF 


YF 

4 ZF 

4 

LF 

4 

MF 

4 

NF 

X 


XW 


W I WM 

4 zw 

4 

alff 

4 

EETF 

4 

alfw 

X 

/MROT C/ 

FXR 

« 

fyr 

« LMR 

4 

LRH 

4 

MRH 

4 

OMR 

X 


AOSS 

* 

ALFSI Y 

^ UTS I Y 

4 

BETS I 

4 

BDTSI 

4 

A 1 SS 

X 


AD1SS 

* 

B 1 SS 

* BD1SS 







X 

/SHIP P/ 

IXX 

* 

I YY 

« IZZ 

4 

MASS 





X 

/TAIL C / 

alfhs 

* 

BTVS 

« ZHS 

4 

YVS 

4 

VTR 

4 


X 


WHS 

♦ 

YTR 

4 DELE 








EQUIVALENCES 


EQU I VALENCE 

X ( DER I NT < 1 , 1 ) *P ) * 

X ( DE R I NT ( l , 2 ) ,Q ) 

SUM FORCES 


C 


XA 

= 

XF 

+ 

XW 

+ 

FXR 


YA 

= 

YF 

+ 

YTR 

+ 

FYR + YVS 


ZA 

= 

ZF 

+ 

ZW 

+ 

ZHS - LMR 








COMPUTE MOMENTS ABOUT C*G. 


LA 

= 

LF 

+ 

LRH 

+ 

DZ*FYR + (DZ - DZTR)#YTR + (DZ - DZVS)*YVS 

C 





+ 

I r*omeg*q 


MA 

— 

MF 

+ 

MRH 

- 

DZ*FxR - (DZ - DZW )*XW + (DXHS - 

DX )*ZHS 






+ 

(DXW - DX )*ZW + DX*LMR - IR*OMEG*P 


NA 

“ 

NF 

+ 

QE 

+ 

DX*FYR - ( DXVS - DX)*YVS - (DXTR - 

DX )*YTR 







COMPUTE ACCELERATIONS 


AX 

= 

XA/MASS 




AY 

= 

YA/MASS 




AZ 

= - 

■ ZA/MASS 





C 


RETURN 

tNu 



FUS 


Subroutine FUS contains the equations and functions used in calculating the fuselage 
forces and moments. FUS also calculates the wing lift, drag forces, and downwash veloc- 
ity (refs. 1 and 2 and appendix A). 


CrjD 

I 

Compute Angles of Sideslip and 
Attack, Aerodynamic Functions, 
Fuselage Forces and Moments, 
Wing Forces, and Wing Downwash 


( RETURN ) 
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SUBROUT INF PUS 


r 

c 

c 

c 

c 

r 


r 

c 


r 

c 

c 


c 

c 

c 


r 

c 

c 

c 


c 

c 

c 


PROGRAM COMPUTES FORCES AND MOMENTS FOR 
FUSELAGE AND WING COMBINATION 

type statements and dimensioned variables 

PEAL LF. MF* NFf IW 
LOGIrAL LDJSI, LDI^O 

integration COMMON TCAT l ON 


COMMON 

X /INTCOMM/ T * DT • I NT . 

X ISCmEMF, DFR 1 NT (2*1?) 

X /REALTIM/ ADC ( 32 ) * DAC<64), L01S1I108), 

x NOPeR . NHOLD * NRESFT » 

X NPR I NT . NREAD 


NEO . 

LDI S0< 196 ) » 
NTERM . 



subroutine communication 

/F CNN AME / CLFJ • CMF | « CNF) , CXI , CX2 » CYTR ♦ 

CYF , CZF . CXW « CZW « CYVS * CZHS . 

CDS I Y « CLStY , DW , GFV , GEH 

/FUNCS / F00 l ( 28 ) « FOO2<20>. F003M4). F 00 A ( 1 4 ) ( FPO? I I 1 * 

F006<29>« FOOT ( ] 4 ) * F008I28)* F009 C 1 A ) , F0I0I28), 

FO 1 1 < 35 ) * FO I 3 ( 07 ) • F014I07J* F015(07>« 0018(27), 

DOO 1 ( 09 ) « D002 < 09 ) * 0003(09). D004I09). D005I09). 
DOOft ( 09 ) . 0007(09). 0008(09). 0009(09). 0010(09). 
DO! 1(09). DO 1 3 ( 09 1 . 00 1 A (09). DO 1 5 I 09 ) . D016MB). 

0017(18). FO 1 6 ( 35, t 0 ) • F0l7(35«10). FO J 3 I 7 . 3 . 4 ) 


COMMON 


X 

/AD V 

C 

/ OMEG 

. WI M 

« VT 

* ov 

, OL 

4SPD 

X 

/FUS 

C 

/ yp 

. YF 

, ZF 

4 LF 

, MF , 

NF 

X 



xw 

. WI WM 

. ZW 

4 ALFF 

4 BETF , 

ALFW 

X 

/FU5 

P 

/ CLF? 

. C MF 2 

. CNF2 

. CKRF 

* CKwIWM, 

I W 

X 



VLF 

, VMF 

. VNF 

4 SXF 1 

4 SXP-2 , 

SYP 

X 



SZF 

4 SXW 

. SZW 





EQUI VALENCES 

EQUIVALENCE 

X (DERINTI 1 , 1 >*p )♦ 

X ( DFP I NT I 1 , 2 > .0 ) . 

X (DFR INT ( I f 3 ) .R ) , 

X (OPPINT( l ♦ 7), l) >, 

X ( OER 1 NT ( 1 , 8 ) . v > . 

X (DEPINT ( ] * 9) . W ) 

COMPUTE SIDESLIP AND ANGLES C* ATTACK 

UL * U 

IE(U*EQ.O,) UL * . oOOO 1 
ALFE * ATAN2 ( W~C<PF*W I M,UL I 
BETF a A T AN 2 ( V • UL ) 

ALFW * ALFE + IW 

WHITE LIGHTS 13.14.15 INDICATE ALFF. BETF. 
ALFW RESPECTIVELY BEYOND LT'MT 

LD I So ( 73 ) * .F. 

LD I SO ( 7* ) * .F, 

LD I So ( 75 ) * .F. 

1 F ( AlFF.LT«-3. 1 2.0P. ALFF.GT* 3* 1 2 ) LD1S0(73) « ,T» 

I F (BETF • LT*-3* 1 2 . OR • OE TF *GT » 3* 1 2 j LD I SO ( 74 ) = *T, 

IF(AlFW. LT. -3,12. OR. ALFW.GT. 3,12 ) LDISC(75) = .T. 

COMPUTE CMF1, CZF. CLF1. CNF1. CYF, CZW 

ALFFd ■ ALFF 

bftfp ■ bftf 

ALFWP ■ ALFW 

ALFSlGN *-S IGN( 1 , .ALpFR) 

BFTSfGN *-S I GN ( 1 . .BETFP ) 

ALFWSGN --SIGNM • .ALFWP) 

ALFFP • ALFF 
BFTFP * — AB S ( BFTF > 

ALFWP * — A0S ( ALFW ) 

IF ( AlFFP*LE . - 3* 1 2 ) ALFFP--3 , J l 99 

IFtALFFP.GE. 3,12) ALFFP* 3.1199 

CALL FUNC1 (FOOl .D001 ,ALFKP) 

IF (BFTFP, LE.-3, 1 2 ) BF TFP *-3 . l l 99 

I E (BfTFP « GF #0.00 ) BpTFPs-t *0F- I 0 

IF(ALFWP.LF.-3*12 ) ALFWP — 3.1 t 99 

IF(AlFWP.GF 40400 ) ALFWP*- 1 *0F- I o 

ALFFP * -AB S ( ALFF ) 

IF (AlFFP.LE.-3, 1 2 ) ALFFP«-3. 1 1 99 
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c 

c 

c 


r 

c 

c 


c 

c 

r 


c 


IF ( AlFFP.GE.O ,0 ) ALFFP=- 1 .OF-1 o 

CALL FUNC 1 (F003«D003*ALFFP ) 

CALL FUNC 1 (F004 *0004 ,BFTFP ) 

CALL FUNC 1 (F005 *0005 ,BFTFP ) 

CALL FUNC 1 (F007 *0007 ,BFTFP l 
CALL FUNC 1 ( F009 * 0009 , ALFWP ) 

CMF1 * DOO 1(1) 

CZF = D003 ( 1 )*ALFS fGN 

CLF) ^ 0004 ( 1 ) *BETS IGN 

CNF 1 = DOOB ( M^B^TFIGN 

CYF = D007 M ) *BFTS TON 

CZW = D009 ( 1 ) *ALFW4GN 

COMPUTE CX2* CXI * CXW 

ALFFP * AFS(ALFF) 

BETFP = ABS(BFTF) 

AL^WD = APS ( AL r W ) 

IF ( A(_FFP. LF* 0*000 ) ALFFP r 1 • 0E~ 1 0 

!F (ALFFP. GE • 3 • 160 ) A L C FP= 3 « 1 599 

IF (0f=-TFP*L^*O*OOO ) BFTFP=1 .OE- 1 O 
J F (BrTFP*GE.3« 1 60 ) BETFP=3*1 59 9 
1 F ( Al c WP*LE* 0 *000 ) ALFWP=1 *0E-10 
I F ( AlFWP*GF. 3. 1 60 ) ALFWP = 3.1599 
CALL FUNC I ( F002 *0002 , ALFFP ) 

CALL FUNC I (FOOfe *0006 *BETFP ) 

CALL FUNC 1 ( F008 « 0008 , ALFWP ) 

C X 2 ■= 0002 ( 1 > 

CXI = 0006 < I ) 

CXW 3 0008 ( 1 ) 

CALCULATE FUSELAGE FORCES AND MOMENTS 


LF 

= 

QL*VLF*CLFt 

+ 

U*P*CLF2 

MF 

= 

ov*vmf*cmfi 

+ 

U*Q*CMF2 

NF 

= 

ql*vnp-*cnfi 

+ 

U*R*CNF2 

XF 


QL*SXF1 *CXi 

+ 

QV*SXF2*CX2 

YF 

= 

QL*SYF#CYF 



ZF 

= 

Q V* SZF*CZF 




CALCULATE WING FORCES AND WING DOWNWASH 

ZW = QV*SZW*CZW 

XW = QV*SXW*CXW ♦ TW*ZW 

WIWM * CK W I WM*U*CZW 

return 

END 
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MAINROT 


Subroutine MAINROT contains the equations and functions used to calculate the 
main rotor forces and moments generated due to a given collective and cyclic pitch and 
to given helicopter angular and translational rates. This subroutine includes the calcu- 
lation of lift and drag coefficients, coning angle and flap angle rates, and total main rotor 
forces and moments including ground effect influences (refs. 1 and 2 and appendix A). 
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c 

c 


SUBROUTINE MA 1NR0T 


c 

r 

r 

r 

c 


c 

c 


c 

c 


c 

c 

c 


c 

c 

r 

c 

c 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


PROGRAM COMPUTES FORCES AND TORQUES ON 
MAIN ROTOR FOR GIVEN COLLECTIVE AnD CYCLIC 
PITCHES* ANGULAR AND TRANSLATIONAL RATES. 


TYPE STATEMENTS AND DIMENSIONED VARIABLES 

REAL LMR. LRH, MRH, |_POR * MPDR 
LOGIrAL LDISI. LOISO 


integration communication 


COMMON 


X 

/ I NTCOMM/ 

T , 

OT , (NT 

. NEO 

, 

X 


| SCHFMF ♦ 

DFRINT (2,121 




X 

/REALTIM/ 

AOC ( 32 ) • 

DA0<64>. LDISI (108), LDIS0<196) 

, 

X 


NQPfR » 

N»- D_D , mrfSET 

. nterm 

, 

X 


NPR I NT , 

NR (r AD 






SUBROUTINE COMMUN 

icat I ON 



common 







X 

/FCNNAME/ 

CLF I 

CMF1 ; CNFi 

, CXI 

♦ CX2 

, CVTP 

X 


CYF , 

CZF , CXW 

, CZw 

, CYVS 

, CZHS 

X 


CDS I Y 

CLSIY , OW 

, GEV 

, GEH 


X 

/FUNCS / 

FOO \ ( 28 ) . 

F002<2e). FOO 3 ( 1 4 ) , FOO 4 (14 1, 

F005 (14), 

X 


FOOf, ( 28 I . 

FOO 7(14), F008 ( 20 ) » F009 ( 1 4 ) , 

F 0 1 0(28), 

X 


FO 1 i < 3*5 > « 

FO 1 3 ( 07 ) , FO 1 4 { 0 7 ) • F015(07)» 

0018(27 ) , 

X 


000 i(09). 

0002(09), DOO 3 ( 09 ) , 0004(09)* 

0005(09 ) , 

X 


000^ ( 09 ) , 

0007(09), 0008(09), 0009(09). 

0010(09 > * 

X 


001 1 (09) , 

0013(09), 0014(09), 0015(09), 

0016 (18), 

X 


00 1 7 ( 18), 

^0)6(3*5,10), 

FO l 7 ( 3? , 

,10), FO 

118(7,3,4) 

common 







X 

/adv C / 

OMEg 

WIM , vt 

• ov 

. QL 

, ASPO 

X 

/CONTRL / 

AOS . 

AIS .BIS 

, THTR 

« OMEGD 


X 

/ftrmpls/ 

EGLmR , 

EGA SS , EQW! M 

, HDOTD 



X 

/gparam/ 

DELY , 

QHO « G 

, PIE 



X 

/MROT C/ 

FxR « 

e YR , LMR 

, LRH 

. MRH 

, GMR 

X 


AOSc , 

ALFSIY, UTS1Y 

, 9E T S I 

. BDTS I 

, A l SS 

X 


ADlcS , 

81 SS , BOI ss 




X 

/ROTCON / 

YPE (20) ♦ 

YT W ( 20 ) • W 1 F ( 20 ) 

, 0OVN 



X 

/ROT P / 

TWIST , 

F MRS , CKAO 

, CKAB1 

. CKAB2 

, CK Aft 3 , 

X 


NR AO , 

NAZ , NB 

, RB 

. efh 


X 

/RCOEFF / 

CKL (20 > , 

CXDI20), CKO ( 20 ) , CKDL 

(20), CKQLI20), 

X 


CKM (20) . 

Y < 20 > , WF ( 20 ) 




X 

/SCPSIR / 

SINciR(5C), COSS I R ( 50 ) 





FQUI VALENCES 


EQUIVALENCE 


( DFRINT II , I) ,d 
< DFRINT I 1 , 2 ) ,Q 
( DFRINT f 1 , 7 ) ,U 
( DERI NT ( I , 8 I « V 
(DERINTl l , 9 ) ,* 
f DER INT (I * I 3 ) .H 


SLDR 

SLSD 

SLCD 

SXOR 

SYDR 

SQDR 

SMDR 

SMSD 


a 0.0 
« 0.0 
* 0.0 
- 0.0 
* 0,0 
a 0*0 
■ 0.0 
* 0.0 
« 0.0 


DO 20 I * 1 . NA Z 


INITIALIZE ACCUMULATIONS AROUND AZIMUTH 


MAIN loop around AZIMUTH SUMS FORCES AND 
MOMENTS ON ROTOR, 


BLADE PI TCH-EXCLUDING TWIST 
TEMT - AOS - AIS *C0SSIR<1) - BIS *SINS!R(1) 
FLAP ANGLE 

SETSI ■ AOSS - A 1 SS*C0SS I R I I ) - B 1 SS*S I NS I R < ! > 
FLAP RATE 


( AOSS— • 043 1 


BDTS! - f A l SS*SINSIR( M - B 1 SS*COSS I R ( 1) > *OMEG 
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RESOLVED VELOCITIES AND SPIN RATE 


USVC ■ U*SINS|RM) 4- V*COSSIR(M 
UCVS » U*COSS I R ( I ) - V*S 1 NS I R ( I ) 
OCPS * 0*COSSrR(l) + P*StNSIR(I) 
WTEM = W - BE TS I *UC VS 


INITIALIZE SUMS OVER RADIAL STATIONS 


INNER loop SUMS forces along blade 


DO I O J ■ 1 , NOAD 


CONING ANGLE AND RESOLVED VELOCITIES 

THSIY * TEMT + VTW(J) 

UTS I Y - YPE(J)*OMEG 4- USVC 

WISIY a W I M*WF ( J ) * ( 1 • + UC VS*W I F ( J ) ) 

UPS I Y 3 VPE(J)*OCPS - Y ( J 1 »BDTS I 4- WTEM-wISlY 


IF ( UTS I Y *LE • 0.0) GO TO 10 
PH I S I Y * UPSIY/UTSIY 
AlESIY* THSIY -4 PHISIY 


ALF5IYP * ALESIY 

UTS I vP = UTS l V 

IF (ALES I YP.LE.-.SOfiO* > 

i f ( ales ! yp.ge. .5060«> 

IF (UTS I YP.LE. 1 1 O. ) 

IF (UTS 1 YP.GE. 1 1 00. ) 


ALF5I YP*-. 50604 
ALES I YP* .50604 
U T S|VP=109 #99 9 
UTS I YP« | 099.999 


WHITE LIGHT 10 OR 19 ON INDICATES 
tjTSIY OR ALESIY RESPECTIVELY IS BEYOND 
SET LIMITS 


LDIS0<^9) » . E * 

IEIUTS1Y *LT • 1 1 O.OnO.OR.UTS I Y .GT.11OO.0l 
IE (ALES IV. LT.-. 50605 .OR. ALES I Y.GT. ,50605 1 


LD I SO ( 7B ) * ,T. 

LD I SO < ^9 ) = .T. 


GENERATE LIFT AND DRAG COEFFICIENTS 


CALL FUNC 2 (FO 16.001 6 , A LF S I Y° « UTS ! YP ) 
CALL FUNC 2 (FOl 7 . DO 1 7 , A LE S I YP ♦ UTS l YP > 
CDS I v z DOl 6 ( 1 I 
CL S I v = DOl 7 ( 1 1 


SUM FORCES AND MOMENTS ALONG BLADE 

CKL (J >*CLSI Y#uTS I Y#UTS I Y 4- LPDR 
( CKD ( J ) #CDS I Y*UTS I Y - CKDL ( J 1 *CLS1 Y*UPS I Y ) 
*UT S 1 Y + DPDR 

CKM(J )*CLSI Y*UTSI Y*UTSI Y + MPDR 

(CKO f J )*CDS I Y*UTSI Y - CKOL ( J I *CLSI Y*UPS I Y » 

*UT S 1 v + OR DR 


COMPLETE SUMMATION AROUND AZIMUTH 


SLOP 

* 

LPDR +• SLOP 


SLSO 

3 

LPOR*S I NSTR ( 1) 

4- SL SO 

slcd 

3 

LPOR*C^>SSIR ( I ) 

+ SLCD 

SMOQ 

3 

MPDP + SMOR 


SMSD 

» 

MPDR*S ! N5 I R f I ) 

4- SMSD 

SMCD 

» 

MPDP*CnSSI R ( I > 

4. SMCD 

SXDR 

* 

BETSI *LPDR*COSSTR ( I ) - 

SYDR 


-betsi*lpor*sinsir< ! ) ■ 

SODR 

« 

QPDR + SODR 



20 CONTINUE 


POVN = RHO/NAZ 
rbon 3 RH04R0VN 


COMPUTE CONING ANGLE AND FLAP ANGLE RATES 


EQAOSS * ROVN*SMDR*rKAO/ ( OMEG*OMEG 1 - AOSS 

ADI Sc » CKAR J *SMSD*OOVN/OMEG - CKAB2*Q + C<AB3*B1SS 

BDlSC = -CKAB1 #SMCO*ROVN/OMEG - CKAB2*P - C<AB3*A1SS 

I E (LD ! S I I17j) GO TO J 00 

AOSS * ROVN*SMDR*CKAO/(OMEG*OMEG ) 
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r> n o 


IFUnT.LF.I) AlSS = AlSS + DELT*AD1SS 

IF(!nT*LE.1) BISS = B1S.S + DEL T *BD 1 SS 

too CONTINUE 


C 

C 


c 


c 


c 


COMPUTE ROTOR FORCES AND MOMENTS 
INCLUDING GROUND EFFECT I NFLUFNCF 


HP = H 
UP = U 

IF(HD,LE.0,OO) hp=.oooi 
IF (HP. GE. 50*0 ) HP=4Q.999 

! F ( UP • L E • 0 « 00 ) UP=.000l 

I F ( UP • GE • 200 • ) UP= 1 Q9 • 999 
CALL FUNC1 (F01 5,09! S , HP ) 
CALL p UNC I t FO 1 4 « DO! A «UP ) 
GFH = DO 1 S ( 1 > 

GFV = 0014(1) 


LMR * RBON* fl • + GFH*GEV ) *SLDR 
EOWIm = W T M- *000 328* LMR/(PHO*VT) 

FYR = RBON*SYDR 

FXR = RBON*SXDP 

LRH =-EFH#SLSD*RBON + OM£G*OMEG*FMRS*B 1 SS 

mrh =-efh#slcd*rbon + OMEG*OMEG*FMRS*Al ss 

QMR = RRON*SOOR 


RFTUON 

END 
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PQ 


Subroutine PQ contains the equations used to calculate the constants needed for the 
convolution solution of a linear second-order equation. The subroutine requires the inte- 
gration step size and the roots of the second-order system in order to calculate the con- 
stants (ref. 7). 
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c 

c 

c 

c 

c 

c 


c 


c 

c 

c 


c 

c 

c 


SUB ROUT I NE PQ ( ONE * T WO • THREE * NEEDRTS *P1 1 * P 1 2 * P2 1 • P2 2 • Q 1 1 ♦ Q 1 2 * Q2 1 * Q22 ) 

TYPE STATEMENTS AND DIMENSIONED VARIABLES 
LOGICAL NEEDRTS * LD1SI * LDISO 

integration communication 

COMMON 

X /INTCOMM/ T * DT * I NT * NFC * 

X ISCHEME* DE R I NT (2*12) 

X /REAlTIM / ADC (32)* DAC(64)« LDISI(IOB)* LDISO (196)* 

X NOPER * NHOLD ♦ NRESFT * NTERM • 

X NPRINT , NREAD 

H = DT 

SIGMA = O* 

SIC-MA1 = 0* 

SIGMA2 = 0* 

I F ( T HRFE • NE • 0 • 0 ) GO TO 30 

REAL EQUAL ROOTS 

10 SIGMA = ONE 
20 SH = SIGMA*H 

I F (S I GMA*EQ*0 ) GO TO 50 
PSH = EXP(SH) 

SIGMAS = SIGMA*SIGMA 

PI 2 = PSH*H 

P21 = — S I GMA S*P 1 2 

PI 1 = PSH - S I GMA *P 1 2 

P2p = PSH + SIGMA+P12 

012 = ( H* ( 1 • + PSH) - 2 • * ( PSH - 1 • ) /S I GMA ) /S I GM A S 

021 = P 1 2 

oil = ( l . - PI 1 )/SlGMAS 

022 = Oil 

GO TO 60 

REAL DISTINCT ROOTS 


30 

S I GMA 1 

= 

ONE 


S I GMA 2 

= 

TWO 

40 

SI H 


S I GMA 1 *H 


S2H 

= 

S IGMA2*H 


PS 1 H 

= 

EXP ( S l H ) 


PS2H 


EXP ( S2 H ) 


SI S2 

= 

SIGMA 1 *S I GMA 2 


S I MS2 

= 

SIGMA 1 -S I GMA2 


I F (SI S2 .EQ.O. •OR. SI MS2 .EQ.0 • ) GO TO 50 
P12 = ( PS 1 H - PS2 H ) /S 1 MS2 

P21 =- S 1 S2*P 1 2 

PI 1 = (SIGMA 1 #PS2 H - SI GMA2*PS1 H ) /SI MS2 

P2P = (SIGMAl#PS]H - S I GMA2*PS2H ) /S 1 MS2 

PI Ml = PS 1 H - 1 • 

P2M1 = PS2H - 1 . 

Q ! 2 = (P1M1/SJGMA1 - H)/(S1MS2*SIGMA1 ) - ( P?M 1 / S I GM A2 - H)/ 

X (SI MS2*S IGMA2 ) 

Q2t = PI 2 

Oil = (PI M 1 /S 1 GMA 1 - P2M 1 /SI GMA2 ) /S l MS2 

022 “ Q l 1 

GO TO 60 

50 L D I SO ( 90 ) - • T • 

60 CONTINUE 

RETURN 

END 
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PQ1 


Subroutine PQ1 contains the equations used to calculate the constants needed for the 
convolution solution of a linear first-order equation. The subroutine requires the inte- 
gration step size and the root of the first-order system in order to calculate the con- 
stants (ref. 7). 



SUBPOUT I NE PQ 1 ( W « P « G 1 « 02 ) 

C 

c INTEGRATION COMMUNICATION 

c 

COMMON 

X /INTCOMM/ T * DT « I NT • NEQ 

X I SCHEME * DFRINT(2,12) 

C 

P = EXP ( — W* DT ) 

Q1 = ( 1 #-P ) /W 

Q2 = (DT-01 )/W 

C 

RETUPN 

END 
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SAS 


Subroutine SAS contains the equations used to evaluate the stability augmentation 
system transfer functions and to calculate the corrections made to the control inputs. 
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SUR ROUT INF S AS 


PROGRAM EVALUATES STAHILITv AUGMENTATION 
SYSTEM TRANSFER FUNCTION^ AND CAlCULA T i G 
THE CONTROL SURF AO nfPUTT IONS 


C OMMON 


COMMON 


I NTFGRAT I ON COMMON I CAT I ON 


/INTCOMM/ T 


I SCHEME * DE R I NT ( P % \ P ) 

SUSP OUT I NE COMMUNICATION 


X 

/ANGS as 

/ 

PH I L 

n 

THETAL 


PS IL 

♦ PH I DL 

♦ THPTOL 

* 

I OL 


X 

/C NTS AS 

/ 

A 1 S CS 

♦ 

HI SCS 


THTRCS 






X 

/CNTL C 

/ 

AOSC 

< 

A 1 EC 


R 1 SC 

♦ THTRC 

♦ A 1 SCL 

t 

t." 1 SCL 

♦ 

X 



THTRCL 










X 

/SAS P 

/ 

Aid 

♦ 

A 1 C2 


Aid 

♦ A 1 < 2 

♦ A I ^ ^ 

♦ 

A 1 K4 

♦ 

X 



B 1 C l 

♦ 

R 1 CP 


R1K1 

• Tl< ? 

t P 1 vr ? 

♦ 

P 1 N4 

« 

X 



THC 1 

• 

THC 2 


THK 1 

• T H xp 

« TH^3 

♦ 

TH<4 


X 

/values 

/ 

A I C g 1 

♦ 

A 1 CS2 


R1 CS 1 

♦ R1CS2 

♦ THC S 1 

♦ 

T w C r - 2 


X 

/X I C 

/ 

A 1 X 1 

* 

A 1 X? 


A l X3 

♦ H 1 X 1 

♦ P 1 X 2 

« 

E 1 X? 

* 

X 



T H X 1 


THX2 


THX3 

♦ A X 1 T N 

* axptm 

» 

^ X 1 T N 

• 

X 



HX2 TN 

♦ 

TX 1 TN 


TX2TV 







EQU I VALENCES 


FOUI VALENCE 


( DERINT ( 1 , 4 ) ♦PHI ) 

( DERI NT < 1 , 5 ) * THE T A ) 

( DFR I NT ( 1 , 6 ) *PS I ) 


CALL T RANT N 1 ( A 1 SC « A 1 SCL • A I C 1 ♦ A 1 C2 * A 1 CS 1 , Axi TN. Ay^TN ) 

CALL TRANFNH < PH I * PH I L * PH I DL ♦ A 1 < 1 * A 1 K ? ♦ A 1 < 3 , A 1 K 4 < A 1 C S •? , 

1 A 1 X 1 • A 1 X? ♦ A 1 X~) 

A 1 so S = A 1 CS 1 - A 1 o? 

call TRANFN 1 (HI S C « B 1 SC L * H 1 C 1 « H l C2 ♦ R 1 CS 1 « BX 1 TN « BX2TN ) 
call TRANFN2 ( THETA . THRTALi THE TDL *B1 K I ♦ R 1 K 2 • HI K ? , R ! K 4 ♦ R 1 C S 2 • 
1 Bl X 1 »R I X? * R 1 X “> ) 

H1SCS = B1CS1 + R1CS? 

CALL TRAiMFNl ( THTRC* T HTRCL ♦ T H C 1 »THC?f THC.S1 * TX 1 TNUxPTN ) 

CALL TRANFN2 (PS I • D S I L « PF I OL * T HK 1 ,THK?iTHK3i THK* • THC.S2 , 

1 THX1 * THX2 , THX2 ) 

THTRCS = THCS1 F THC^2 


RETURN 

FNO 
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SETUPA 


Subroutine SETUPA is used by the trim circuit algorithm and initializes XT RIM to 
the starting values of the trim parameters. During the trim iterations, the trim parame- 
ters are set to the XT RIM values (ref. 6 and appendix D). 



SUBROUT INE SETUPA ( X ♦ Y 1 • Y2 * Y 3 ♦ Y4 , Y5 « V6 * Y 7 « Y8 , Y 9 « Y 1 0 * V 1 1 .N ) 


C 

C 

C 

C 

C 

C 

C 

C 

C 


C 


program initializes xtrim to THE STARTING 
VALUES OF THE PARAMETERS* DURING THE TRIM 
ITERATION THE TRIM PARAMETERS ARE SET TO 
THE XTRIM VALUES 

TYRE STATEMENTS AND DIMENSIONED v^RIABl^S 

D I MENS ION X ( 1 1 ) 

GO TO ( 1 00 i 200 ) « N 
100 CONTINUE 


X ( 1 ) = 

X ( 5 ) = 

X ( 9 ) = 

RETURN 

Y1 

Y5 

Y9 

SX< 2) = Y2 

$X< 6) = Y6 

$x(10) = Y 1 0 

$X ( 3 ) 

$X ( 7 ) 

$x(in 

= Y3 
= Y7 
= YU 

$X ( 4 ) 

SX ( 8 ) 

= y 4 
= Y8 

CONTINUE 
Y 1 = X ( 

Y5 = X ( 
Y9 = X ( 

1 ) 
5 ) 
9 ) 

SY2 = X( 2) 
$Y6 = X( 6) 
SY 1 0= X < 1 0 ) 

$Y3 = 
SY7 = 
$Y 1 1 = 

X( 3 ) 
X( 7 ) 
X ( I 1 ) 

SY4 = 
$ Y 8 = 

X ( 4 ) 

X ( 8 ) 


C 
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RETURN 

END 





SETUP2 


Subroutine SETUP2 is used by the trim circuit algorithm and equates FTRIM to the 
trim equations (ref. 6 and appendix D). 



SUBROUT INE SETUP2 (F *Z1 . Z2 . Z3 . Z4 , Z5 . Z6 « Z7, Z8 . Z9 . Z1 0 « Z1 1 ) 

C 

C PROGRAM EQUATES FTRIM TO THE TRIM EQUATIONS 

C 

C TYPE STATEMENTS AND DIMENSIONED VARIABLES 

C 

D I MENS ION F ( I 1 ) 


F ( 

1 ) = 

Z1 

$F ( 

2 ) 

= Z2 

$F ( 3 ) 

= Z 3 

SF ( 

4 ) 

= 74 

F ( 

5 ) = 

Z5 

SF ( 

6 ) 

= Z6 

*F { 7 ) 

= Z7 

SF ( 

8) 

* Z8 

F ( 

9) = 

Z9 

SF (10) 

= Z 10 

SF ( 1 I ) 

= Z 1 I 





C 

RETURN 

END 
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SIMEQA 


Subroutine SIMEQA is used by the trim circuit algorithm and solves a system of 
linear equations by the Gauss- Jordan reduction method (ref. 6 and appendix D). 
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SUBROUTINE S I MEOA ( A , N , CC , XMA T , XTR I M ) 


C 

C PROGRAM SOLVES A SYSTEM OF LINEAR EQUATIONS 

C BY THE GUASS-JORDAN REDUCTION METHOD IN THE 

C TRIM ALGORITHM 

C 
C 

C TYPE STATEMENTS AND DIMENSIONED VARIABLES 

C 

DIMENSION A (12,13)* CC ( 12), XMAT(11«I2)« XT«IM(11> 

DOUBLE PRECISION B<1?«13)« C<12), DXMAT ( 11,12), DXTRIM(ll), D, TOL 

C 

TOL * 5 • D— 20 

KA * N 4- 1 

K2 * 2 

K 3 * 1 

DO in I *1 ,N 
DO 10 J«1 , K A 
C 

C MOVE MATRIX TO WORKING AREA 

10 B< I • J) * OBLE (A(l,Jn 

20 I F ( DABS (B(K3tK3)>- TOL) 30* 1 00 , ! 00 

30 I F (Kp - N) 40, 40, pOO 
40 DO 50 I =1 ,KA 

D * B <K3, I ) 

B(K3, I ) * B(K2, I > 

B(K2, I > * D 

50 CONTINUE 

K2 ■ K2 4 1 

GO TO 20 
100 K 1 * K3 

K3 * K3 4 1 

K2 * K3 4 1 

I F (K 3 - N) 120, 1 20 i 200 

120 DO 130 I *K3 * KA 

B(K1,I) * B CK1 , I )/B <*l ,K1 ) 

130 CONTINUE 

DO 140 I*K3,N 
DO 140 J*K3,KA 

B(!,J) * B ( I * J ) - B f I «K1 )*B (K1 • J) 

140 CONTINUE 
GO TO 20 

200 C(N) ■ B(N*N41 )/B(N,N) 

K3 * N - 1 

210 K 1 * K3 4 1 

D - O, 

DO 2?0 I *K 1 « N 

220 D ■ D 4 8 { K3 • I )*C { I > 

C (K3 ) « B ( K3tKA ) - D 

K3 ■ K3 - 1 

IF(K3) 350, 350, 2lo 

300 TOL * TOL/ l 0 , 

I F t TOL - 1*0-30) 320,310,310 
310 K2 * K3 4 1 

GO TO 20 
320 CONTINUE 
350 CONTINUE 

NM1 * N - 1 

DO 3ft0 1*1, Nf » 1 
DO 380 J*1 ,N 

360 DXMAT ( I , J ) * DBLE(XMaT( I , J) ) 

DO 370 1*1 , NM1 

DXTR I M ( I ) * O.ODO 

DO 3^0 J*| , N 

DXTRlM(I) * DXTRIM(I) 4 DXMA T ( I , J ) *C ( J ) 

370 CONTINUE 

DO 380 J*1 ,NM1 

380 XTRIm(J) * SNGL IDXTRJM { J) ) 

DO 390 J* l , N 
390 CC(J 1 * SNGL (C ( J ) ) 

C 

return 

ENO 
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TAIL 


Subroutine TAIL contains the equations and functions used to calculate the forces 
and moments generated due to the tail rotor and the horizontal and vertical stabilizers. 
General tail velocities are also computed to account for downwash velocities and angular 
rates (refs. 1 and 2 and appendix A). 
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subroutine tail 


Real ! HSf IVS 
LOGICAL LDIS1. 


PROGRAM COMPUTES FORCES AND MOMENTS DUE TO 
TAIL ROTOR « HORIZONTALt AND VERTICAL STA- 
BILIZERS, GENERAL TAIL VELOCITIES ARE 

also computed 


type STATEMENTS and DIMENSIONED variables 


integration communication 


/intcomm/ t , dt , INT , neo , 

1 SCHEME • DE R I NT (2*12) 

/REAlTJM/ ADC { 32 ) « DACI64), L0ISl<108), LDISOtlRft), 
NOPER , NHOLD , NRESET ♦ nTERM , 

NPR1NT . NREAD 


subroutine communication 


/fcnname/ 


/EUNCS / 


/ADV C / 

/cntl p / 

/CONTRL/ 
/FUS C / 


CLE 1 * 

CYF , 

CDS I Y , 
FOO 1 (20 j 
F006I28 > 
FO 1 1 ( 35 l 
D00 1 (09 > 
0006 (09) 
DO I 1 (09 ) , 
0017(18), 


CMFl , CNF l 
CZF , CXW 
CLSIY , DW 
i F 002 ( 28 ) < F 
F 007 ( 1 A ) , F 
FO] 3 I 07 ) • F 
D002 ( 09 ) • D 
D007 ( 09 ) « 0 
DO 1 3 ( 09 ) « D 


r l • CXI 

* ♦ czw 

4 GEV 

F003 ( 1 A | , 
FOO0 (?8 ) , 
FO 1 A (07 ) 4 
0003 (09)4 
0008 ( 09 ) , 
0014(09), 


4 CX2 
4 CYVS 
4 GEH 
F 004 (14), 
F 009 (1 U > , 
FO 1 5 (07 ) , 
0004 ( 09 ) , 
0O09 (09), 
DO 1 5 ( 09 ) , 


F 0 1 6 ( 3S * 10), FO l 7 ( 35 , l 0 > * FC 

rfIM 4 VT , Ov , OL 

A 1 SCO • 81 SCO , THTRC04 XAOSG 
<TRG » XAOSR , YC SR , XCSR 



XW 

# W 1 WM 

4 ZW 

4 ALFF 

/F N M P/ 

DX 

4 DXHS 

4 DXTR 

4 Dxvs 


DZTR 

. DZ VS 

4 DXW 

4 OZW 

/STICKS / 

XAOS 

4 YCS 

4 XCS 

4 XTR 

/tail C / 

ALFHS 

4 8 TVS 

4 ZHS 

. Y VS 


WHS 

, YTR 

4 DELE 


/TAIL P/ 

CK TR1 

t CKTR2 

» THS 

4 I VS 


SZHS , S Y VS 


ECU I VALENCES 


EQUIVALENCE 


( DERINT ( 1 
( OER I NT ( 1 
(DERINT ( 1 
(OER INT ( 1 
(OERlNT ( I 


COMPUTE DOWNWASH AS A FUNCTION OF 
forward VELOCITY 


UP = U 

I F ( l)P,LE ,0*00 ) UP = 0,0001 

IF (UP.GE, 152,0 ) UP =151, 999 

CALL FUNC 1 (F01 3,D01 3,uP) 

DW = 0013(1 ) 


COMPUTE TAIL area velocities 


VTR = V - ( DXTR - DX)*R 

WHS = W + (DXHS - DX)*Q - DW*WIM - wIWM 

COMPUTE STABILIZER ANGLES OF ATTACK 
(HOR* STAB, MODIFIED BY ELEVATOR) 

IF(LDISJ (37 ) ) XCS - (B1S*DPR - SlSCOl/XCSG 

DELE = (10.6052 - 2.0405*XCS ♦ 0 . 24 45* XC S*XCS 1 *RPD 

UL = U 

IF(u*EQ,0.) uL = .00001 

ALFHS = ATAN2 ( WHS ,UL ) + CKTRt#OELE + tHS 
BTVS = A TAN2 ( VTR , UL ) + IVS 


WHITE LIGHT 16 OR 17 ON INDICATES 
ALFHS OR BT VS RESPECTIVELY IS BEYOND 
SET LIMITS 

LO I SO ( 76 ) = ,F, 

LD I SO ( 77 ) = ,F . 

IF (ALFHS. LT. -3. 12. OR. ALFHS. GT. 3. IP) LDISO<76) = .T. 
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n r'i n 


IF (RTVS.LT.-3.06.0P.«TVS.C,T.3.06 ) 


LOT SO ( 77 ) = .T. 


COM D UTE CYVS. CZHS 


c 

c 

c 


c 

c 

c 

c 


c 


c 


r 

c 


AL^HSP= ALFHS 

BTVSP = B TVS 

ALFSGN= s I On ( 1 • * ALFHSP ) 

BT VSGN = S !GN ( 1 . « 9TVSP ) 


ALFH5P= ABS( ALFHS) 
BTVSP = ABS(BTVS) 

IF (BTVSP.LE.O .00 ) 

IF (3TVSP.GE. 3.06 ) 

IF { ALEHSP . LE .0 *00 ) 
IF ( ALFHSP.GF.3. 12 ) 
CALL FUNC I ( FO 1 0 *D0 
CALL FUNC 1 (F01 l .30 
CYVS = DC 1 0 ( 1 ) *BT 
CZHS = DO 11(1) *ALI 


9TVSP=0.000 1 
BTVSP=3.0599 
AlfHSP=0.000 1 
A L fhSP=3. 1 1 99 
O « BTVSP > 

1 ♦ ALFHSP ) 

'SGN 

SGN 


COMPUTE STAR1LIZFQ FORCES 


Y VS = QL*SYVS*CYVS 

ZHS = QV*SZHS*CZHS 


WHITE LIGHT 20 ON INDICATES U. VTR 
THTR BEYOND LIMITS WHEN COMPijTINC- 


LD I SO ( 80 > = .F , 

IF(VTR .LT.-20.0. 
IFOj # L T • OO.O. 

IF ( THTR.LT.-.353. 


OR.VTR .GT.20.O) 
OR.U • GT « 360 . ) 

OR • THTR.GT . .432 ) 

COMPUTE TAIL 


LD I SO < 80 > = • t . 

LD I SO ( 80 ) = .T. 

LDI SO ( 80 ) = . T . 

RO^OR SIDE FORCE 


UP = u 
V T RP = VTR 
THTRP = THTR 
1 F (UP.LE.0.00 ) 

I F ( UP .GE. 360.0) 

I F ( VTRP .LE • -20 • ) 

I F ( VTRP .GE . 20 • ) 

IFCTHTRP.LE.-.353) 
I F (THTRP. GE • .432 ) 


UP= .000 1 
UP= 359 . 999 
VTRP = — 1 9 . 999 
VTRP = 19.999 

THTRP=- • 3529 
THTRP= .4319 


CALL FUNC 3 ( FC 1 8 » DC 1 8, T H T RP , VTRP .UP ) 
C YTR = DO 1 8 ( 1 ) 

YTR = CKTR2+CYTR 


RETURN 

END 


, OR 
C ytc? 
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u u 


TRANFN1 


Subroutine TRANFN1 contains the equations used to compute the values of the first 
of two transfer functions which make up the total SAS transfer function. 



SUBROUT I ME TRAnFnI ( CNT C , CNTCL « C 1 ,C2,CNTS1 ,Xl TN«X2TN) 

SUBROUTINE COMPUTES VALUES OF THE 
C FIRST TRANSFER FUNCTION 

C 

C TVPE STATEMENTS AND DIMENSIONED VARIABLES 

C 

LOGICAL NEEDRTS 
C 

C INTEGRATION COMMUNICATION 

C 

COMMON 

X /INTCOMM/ T * DT , I NT * NEO , 

X I SCHEME , DE R I NT ( 2 i 12) 

C 

CNTDL = (CNTC-CNTCL ) /DT 
NEEDRTS = • E • 

RTYPE = 1.0 

ROOT 1 = -0.4 
R00T2 = -1 */C2 

CALL PQ (ROOT 1 , ROO T2.RTYPF, NEEDRTS iPl 1 *P 12 *P2 1 .P22 t Oil ,012.021 ,022 ) 
X1TL = XITN 
X2TL = X2TN 

XITN = P11*X1TL + P 12*X2TL + Q11*CNTCL + 0!2*CNTDL 
X2TN = P21*X1TL + P22*X2TL + 02 1 +CNTCL + Q22*CNTDL 
CNTS1 * Cl /(2.5*C2 )*X2TN 
C 

RETURN 

END 
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TRANFN2 


Subroutine TRANFN2 contains the equations used to compute the values of the sec- 
ond of two transfer functions which make up the total SAS transfer function. 



SUB POUT 1 NE TPANFN2 ( AnG . AnGL * ANGDL * Cl #C2 * C3* C4 * CNTS2 • XI * X2 • X3 ) 
C 

C SUBROUTINE COMPUTES VALUES OF THF 

C SECOND TRANSFER FUNCTION 

C 

c 

c INTEGRATION COMMUN I c at I on 

c 

C 0 WON 

X /INTCOMM/ T * DT * INT * NEO * 

X [SCHEME# DER I NT ( 2 # 1 2 ) 

C 

XI (XX#XP#XQ1 #XQ2) = XP*XX + XQ1*ANGL + XQ2*AnGDL 

YI(TC#XJ#AI) = (ANG - l./TC*XJ)*AI 

AP = 2.5 

A 1 = (Cl* (C2-AP) )/( AP*< AP-C3 ) * { AP-C4 > ) 

A2 = (C 1* < C2-C3 ) )/ (C3* ( C3-AP )* (C3-C4 > ) 

A3 = <C 1 * < C2-C4 > >/<C4* ( C4-C3 >* (C4-AP ) ) 

C 

CALL PQ 1 (1 */AP,P«Ql #Q2) 

XI = XI (XI * P • Q 1 #02) 

CALL PQ1 (1 •/C3#P*01 #G2) 

X2 = XI (X2#R#Q1 #02 ) 

CALL PQ1 ( 1 •/C4,Pt01 *02 ) 

c 

X3 = XI ( X 3 * P * Q 1 * Q2 } 

Y 1 = Y I ( AP * X 1 * A 1 ) 

Y2 = Y I ( C 3 • X 2 • A2 ) 

Y3 = Y I (C4*X3*A3 ) 

CNTS2 = YI + Y2 + Y3 
C 

RETURN 

END 

98 




TRIMDER 


Subroutine TRIMDER is used by the trim circuit algorithm and executes all subrou- 
tines necessary to compute the trim equations (ref. 6 and appendix D). 



subroutine: trimder 

c 

C PROGRAM EXECUTES ALL SUBROUTINES NECESSARY 

c TO COMPUTE THE TRIM EQUATIONS 

C 

CALL AEROVAR 
CALL FUS 
CALL TAIL 
CALL MA INROT 
CALL ENGINE 
CALL F ANO M 
CALL BODYDER 
CALL EARTH M 
C 

RETURN 

END 
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XMATRIX 


Subroutine XMATRIX is used by the trim circuit algorithm and sets up a matrix of 
independent trim parameters with the first column being XTRIM, the second column being 
XTRIM with a perturbation on the first element, etc. (ref. 6 and appendix D). 

^ XMATRIX (XT) ) 

V _ 

Set Up Matrix of Independent Trim Parameters 

Let Column 1 = XTRIM, Column 2 = XTRIM + Perturbation 

on First Element, etc. 


( RETURN ) 


c 

c 

c 

c. 

c 

c 

c 

c 

c 

c 

c 


r 


1 


2 

r 


SURROUT I NF XMATR IX(XT) 


PROGRAM SFTS UP A MATRIX OF ino-pfndfnt 
TRIM PARAMFTFRS WITH THF FIRST COLUMN 
BFING XT R I v , THF SECOND UF I NG XTRIM WITH 
A PFRTURBATION ON THF FIRST EL rV,c N T ' F TC * 

TYPE STATEMENTS AND DIMENSIONED VARIABLES 


D I MENS I ON XT ( 1 1 ) 


SUUROUT I NF COWlJN [CATION 


COMMON 

X /TRIMAT/ 

X 


XM A T ( 1 1 * 1 2 ) * t. M A T ( 1 2 ♦ 13) • 
NFQP1 « NE OP 2 « 


TWSIJM ( I 2 ) « NEON * 
XMATG (11) 


DO 1 00 I = 1 * NF ON 

DO 1 00 J = 1 * NF OP 1 
XMAT ( I « J ) =i XT ( I ) 

0 CONTINUE 


no ? 00 J = I ,NFON 
XMAT( J+l ) = XMATG ( J ) *XT ( J ) 

IF (XT (J) • F Q . O • 0 ) X M A T ( J * J + I 

CONTINUE 


• or i 


RFTl IQN 
FNn 
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n r> 


XVAL 


Subroutine XVAL is used by the trim circuit algorithm and sets XTRIM to the val- 
ues in XMATRIX to generate n + 1 equations where n = number of trim equations 
(ref. 6 and appendix D). 


C XVALIXT, NCOL) ) 

Set XTRIM to Values 
in XMATRIX 


( RETURN ) 


SUBPOUT INE XVAL (XT » NCOL ) 

C 

c PROGRAM SETS XTRIM TO THE VALUES IN 

C XMATRIX TO GENERATE N + 1 EQUATIONS 

C 

C TYPE STATEMENTS AND DIMENSIONED VARIABLES 

C 

D! MENSI ON XT ( t 1 ) 


SUBROUTINE COMMUNICATION 


C 

COMMON 

X /TRIMAT/ xMAT (11*12)* EMAT<12*13)« TRSUM ( I 2 ) ♦ 

X NEOP1 * NEQP2 ♦ XMATG(II) 

C 

DO 1 00 J = 1 •NEQN 

XT ( J ) = XMAT(J«NCOL) 

100 CONTINUE 
C 

RETURN 

END 


NEON ♦ 
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PROGRAM USAGE 


This section of the report describes the input data (BLOCK DATA and function data), 
the program setup, and the procedures for obtaining a static check, dynamic check, and 
cockpit checkout. Also included are console and cockpit running procedures. 

Input Data Description 

The function data are loaded into the program on data cards using a FORTRAN IV 
READ statement and a 7F10.0 format. All other input parameters are loaded in the 
FORTRAN IV BLOCK DATA subprogram. This subprogram assigns constant values to 
variables by use of the FORTRAN IV DATA statement. These values are then transferred 
to other subroutines through labeled COMMON. The BLOCK DATA input parameters are 
listed as follows with their appropriate COMMON label and typical value. Input parame- 
ters for function generation are also stored in the BLOCK DATA subprogram. 


COMMON 

label 

FORTRAN 

variable 

Input value 

/ADV P/ 

ADV PI 

0.000328 


ADVP2 

3.0 

/CNTL P/ 

AOSCO 

8.5 


A1SC0 

-10.8 


B1SC0 

-13.2 


THTRCO 

23.0 


XAOSG 

0.92 


YCSG 

1.85 


XCSG 

2.81 


XTRG 

-5.00 


XAOSR 

13.3 
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COMMON 

label 

FORTRAN 

variable 

Input value 

/CNTL P/ 

YCSR 

9.6 


XCSR 

9.6 


XTRR 

6.0 

/ENG P/ 

IROT 

2 820.0 


CKE1 

10 529.0 

/F N M P/ 

DX 

-0.5125 


DXHS 

16.55 


DXTR 

26.75 


DXVS 

25.08 


DZ 

7.584 


DZTR 

2.863 


DZVS 

4.40 


DXW 

- 0.67 


DZW 

7.554 


IR 

0.0 

/FUS P/ 

CLF2 

-20.0 


CMF2 

0.0 


CNF2 

-30.0 


CKRF 

0.50 
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COMMON 

label 

FORTRAN 

variable 

Input value 

/FUS P/ 

CKWIWM 

0.0 


IW 

0.244 


VLF 

100.0 


VMF 

100.0 


VNF 

100.0 


SXF1 

10.0 


SXF2 

10.0 


SYF 

10.0 


SZF 

10.0 


sxw 

27.8 


szw 

27.8 

/GPARAM/ 

DELT 

0.03125 


RHO 

0.00238 


G 

32.2 


PIE 

3.14159 

/RCOEFF/ 

CKL(l) 

5.32 


CKL(2) 

4.95 


CKL(3) 

4.35 
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COMMON 

label 

FORTRAN 

variable 

Input value 

/RCOEFF/ 

CKD(l) 

5.32 


CKD(2) 

4.95 


CKD(3) 

4.35 


CKQ(l) 

52.68 


CKQ(2) 

70.79 


CKQ(3) 

92.57 


CKDL(l) 

5.32 


CKDL(2) 

4.95 


CKDL(3) 

4.35 


CKQL(l) 

52.68 


CKQL(2) 

70.79 


CKQL(3) 

92.57 


CKM(l) 

52.68 


CKM(2) 

70.79 


CKM(3) 

92.57 


Y(l) 

9.90 


Y(2) 

14.30 


Y(3) 

18.70 


COMMON 

label 

FORTRAN 

variable 

Input value 

/RCOEFF/ 

WF(1) 

1.0 


WF(2) 

1.0 


WF(3) 

1.0 

/ROT P/ 

TWIST 

-0.175 


FMRS 

0.0 


CKAO 

0.000704 


CKAB1 

0.000704 


CKAB2 

1.0 


CKAB3 

0.0 


NRAD 

3 


NAZ 

4 


NB 

2 


RB 

22.0 


EFH 

0.0 

/SAS P/ 

A1C1 



A1C2 



A1K1 



A1K2 
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COMMON 

label 


FORTRAN 

variable 


Input value 


/SAS P/ A1K3 

A1K4 
B1C1 
B1C2 
B1K1 
B1K2 
B1K3 
B1K4 
THC 1 
THC2 
THK1 
THK2 
THK3 
THK4 


/SHIP P/ IXX 2 530.0 

IYY 11 716.0 

IZZ 10 164.0 

MASS 274.0 
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COMMON 

label 

FORTRAN 

variable 

Input value 

/TAIL P/ 

CKTRl 

1.0 


CKTR2 

4 000.0 


IHS 

0.0 


IVS 

-0.0785 


SZHS 

13.5 


SYVS 

22.5 
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BLOC* DATA 


BLOCK DATA SUBPROGRAM ASSIGNS VALUES 
FOP PARAMETFRS IN LABELLED COMMON 


C 

c 

r 

c 

c 


c 

c 

r 



TYPE STATEMENTS 

and dimensioned variables 

RFAL ir. tw. 

IMS. TVS. JROT. IXX, 

I YY . 

TZZ, MASS 


SUBROUTINE 

COMMUN I CATION 


Common 


X 

X 

X 

X 

X 

X 

X 

X 

X 

X 


COMMON 


X 

X 

X 

X 

X 

X 

X 

X 

X 


/adv p / 

ADVP 1 

ADVP2 





/CNTL P / 

aosco 

At SCO . 

B 1 SCO . 

thtrco. 

XAOSG . 

YCSG 


XCSG 

XTRG « 

XAOSR 

YCSR . 

XCSR 

XTRR 

/ENG P / 

TROT 

CKEl 





/F N M P/ 

DX 

DX HS « 

CXTR . 

oxvs . 

DZ . 



DZTP 

DZ VS . 

Dxw . 

DZW • 

1R 


/FLIS P / 

CLF? 

C ME 2 

CNF2 

CKRF . 

CKW! WM . 

IW 


VLF 

VMF * 

vnf . 

SXF t . 

SXF2 . 

SYF 


SZF 

sxw . 

SZW 




/GP ARAM/ 

DELT 

RHO ♦ 

G . 

PIE 



/RCOEFF / 

CKL (20 ) 

CKD(20 ) 

. CKQ ( 20 ) 

, C < DL (20). CKOL(20). 


CKM (20 ) 

Y ( 20 ) 

. WFI20) 




/ROT P / 

TWlsT 

FMRS . 

CK AO . 

CK AB L 

CKAB2 « 

CK AB 3 


np An 

NAZ . 

N0 . 

R9 

EFH 


/SAS P / 

A 1 C | 

A I C2 

A 1 K 1 . 

A IK? 

A 1 K3 . 

A 1 K4 


81 C 1 

B I C2 ♦ 

8 1 K I . 

8 I K 2 

R I K3 * 

8 ] K4 


THC 1 

THC? ♦ 

THK 1 . 

THK? 

THK 3 • 

THK4 

/SHIP P/ 

I XX 

, I YY . 

1 77 . 

MASS 



/TAIL P/ 

CKTR I 

, CKTR2 , 

1 H S 

I VS . 

SZHS 

S Y VS 


common 

X /FUNCS / 

X 

X 

X 

X 

X 

X 


FOO l 128)* 
FOOft ( 28 1 . 
FOI| ( 35 ) • 
000 l I 09 > • 
DOOis ( 09 > * 
DO I 1 (09). 
DO 1 7 ( 1 8 > • 


F 002 (28 ) 
F 007 (14 ) 
FOl 3(07) 
0002 (09 ) 
0007(09 ) 
001 3(09 ) 
F0| 6( 35* 


, F003U4). F 00 A 

♦ F008 (28 1. F 009 

♦ F014 (07 | * FOl 5 
« 0003(09). 000 A 
» 0008(09). 0009 

♦ DO 1 4 (09 > . 0015 
10). FO I 7 ( 35 . 1 0 


(14). F 005 (14). 
(14). FOl 0(28 ) . 
(07 ) * 0018 (27) * 

(09 ) . D005 (09 ) . 
( C9 ) . OOl 0 (0 9 1 * 
( 09 ) . 0016(18). 

) . FO 1 8 ( 7.3 *4 ) 


ThF following data is FOR COBRA 


C 

C 

c 


C 

c 

c 


c 

c 

c 


c 

c 

c 


DATA SOURCE IS TECHNICAL REPORT 
ECOW-0387-F2A 


/ADVP / — AERODYNAMIC PARAMETERS 

DATA 

1 ADVP 1 / 3.28E-4 X, ADVP2 / 3.0 / 


/CNTl P / 

— CONTROL 

STICK parameters 






DATA 

1 AOSCO / 

8.6 

/. aisco / 

-10.8 

/ . 

R I SCO 

/ 

-1 3.2 

/ 

2 THTRCO/ 

23. 

/. XAOSG / 

0.92 

/. 

YC5G 

/ 

1 .85 

/ 

3 XCSG / 

2.81 

/. XTRG / 

-5.00 

/ . 

XAOSR 

/ 

1 3.3 

/ 

4 YCSR / 

9.6 

/. XCSR / 

9.6 

/♦ 

XTRR 

/ 

6.0 

/ 


/ eng p 

/ 

— ENGINE 

PARAMETERS 





DATA 










1 I POT 

/ 

2820.0 

/ . 

CKE l 

/ 10529.0 

/ 




/F N M 

P/ 

— PARAMETERS 

FOR 

CAL CUL A T I ON 

OF 

FORCES 

AND MOMENTS 

DATA 










1 DX 

/ 

-.5125 

/ • 

OXHS 

/ 16.55 

/, 

OxTR / 

26.75 

/ 

2 OXVS 

/ 

25.08 

/ • 

OZ 

/ 7,584 

/. 

OZ TR / 

2.063 

/ 

3 DZvS 

/ 

4 .40 

/ . 

DXW 

/ -0.67 

/ * 

DZW / 

7.554 

/ 

4 IR 

/ 

0.00 

/ 








/ FUS 

P /- 

-FUSELAGE 

PARAMETERS 






' 

DATA 
1 CLF •? 

/ 

-20.0 

/ * 

CMF2 / 

0.0 

/ . 

CNF2 

/ 

-30.0 

/• 

2 CKRF 

/ 

0.50 

/* 

CKW! WM/ 

0*0 

/ • 

IW 

/ 

0.244 

/♦ 

3 VLF 

/ 

1 00. 

/. 

VMF / 

mo. 

/♦ 

VNF 

/ 

100. 

/. 

4 SXF) 

/ 

10,0 

/* 

SXF? / 

10.0 

/♦ 

SYF 

/ 

10.0 

/. 

5 SZF 

/ 

10.0 

/* 

SXW / 

27.6 

/• 

SZW 

/ 

27.8 

/ 


/GPArAM / — GENERAL PARAMETERS 
DATA 




OELT / 0.03125 /. RHO / 0.00238 /. G 

/ 

32.2 

/ 


2 PIE / 3. 1 4 1 69 / 



c 

c 


c 

c 

c 


c 

c 


/RCOFFF/ MAIN ROTOR ELEMENTAL LIFT AND DRAG COEFFICIENTS 


DATA IS FOR 3 ELEMENT rotor 

data 


1 

CKL ( 1 

)/ 

5.32 

/♦ CKL <21/ 

4.95 

/« CKL ( 3 1 / 

4.35 

/ 

2 

CKD< 1 

)/ 

5.32 

/. CKD<2)/ 

4.95 

/. CKO ( 3 1 / 

4.35 

/ 

3 

CKO ( 1 

)/ 

52.68 

/* CKQ ( 2 1/ 

70.79 

/. CKO ( 3 )/ 

92.57 

/ 

4CKDL ( 1 

)/ 

5. 32 

/•C<OL (2 )/ 

4.95 

/ « CKDL (3 1/ 

4.35 

/ 

5CKQL ( 1 

)/ 

52.68 

/*CKOi-<2)/ 

70.79 

/♦CKOL ( 3 )/ 

92.57 

/ 

6 

CKM( 1 

)/ 

52.68 

/. CKM f 2 1 / 

70.79 

/« CKM (3)/ 

92.57 

/ 

7 

Y( 1 

1/ 

9.90 

/♦ Y ( 2 1 / 

14.30 

/. v (3 )/ 

18.70 

/ 

8 

WE ( 1 

)/ 

1 .0 

/ . WF ( 2 ) / 

1 .O 

/. WF(3)/ 

1 .0 

/ 


/ ROT p / — MAIN ROTOR PARAMETERS 
DATA 

1 TWIST / — 0 • I 75 /. FMRS /0.0 /. 

2 CKAO / 7.04E-4 /• CKAB1 / 7.04E-4 /« CK AB2 / 1 .0 /. 

3 C«A*3 / 0*0 /» NR AD / 3 /* NAZ / 4 /. 

4 NB / 2 /. RB / 22.0 /. EFH / 0.0 / 

/SAS P / — SAS GAINS AND TIME CONSTANTS 

data 

1 A 1 C 1 / /• A1C2 / /♦ AIK1 / /. 

2 AIK? / /* A I K3 / /• A 1 K4 / /. 

3 B 1 C 1 / /. 91C2 / /♦ BIK1 / /. 

4 B1K? / /. B1K3 / /. BtK4 / /. 

5 TMCl / /. THC2 / /. THK 1 / /» 

6 THK 2 / /» THK3 / /* THK4 / / 


/SHIP P / INERTIA AnD MASS PARAMETERS 

data 

1 I XX / 2530.0 /. ivy / 11716.0 /. IZZ / 10164.0 /. 

2 MAES / 274.0 / 


/TAIL p / — TAIL ROTOR AND STABILIZER PARAMETERS 
DATA 

1 CKTR1 / 1.0 /. CKTR2 / 4000* /♦ IHS / 0.0 /. 

2 IVS / - . 0785 /♦ SZH5 / 13.5 /. SWS / 22.5 / 


INTERNAL DATA F OR FUNCTION GENERATION 


DATA IDOO 1 ( 1 ONE 1 * 1 0Nr = 2 .9 ) / 

-2.6O000E-01 , 3.l200nE+00« 3.I2000E+00. 2.60000E-OI, 3.12000E+00, 

2.60000E-01 .0. ,0. / 

DATA ( D002 ( I ONE ) * I0Nr = 2*9 1 / 

7.84000E-01, I • OOOOOE+OO « -O » ♦ 4.90000E-02. 2.672008+0^. 

2 . 1 ftOOOE-O I * 3.5900OE + 00. 2 . 70000E -01/ 

DATA ( D003 ( I ONE > * IONcr = 2.9 1 / 

-2.60000E-01 • 0* , 3 • 1 2000E + C0 , 2 .60000E-0 1 » 3. 1 200PE+0P . 

2 . 6O000E-0 t * 0* ,0. / 

DATA ( DOO 4 < 1 ONE ) * I ONc-* 2 ♦ 9 l / 

-2.60000E-01 * 0. * 3. 1 2000E+00, 2.60000E-01. 3.12000E+0P. 

2 • 60 OOOE -O 1 * 0 . .0. / 

DATA ( D005 < I ONE ) , I0Nr*2 .9 ) / 

-P.60000E-0 1 * 0. ♦ 3. 1 2000E + 00 . 2.60000E-01. 3.12000E+0C. 

2.60000E-01 . 0. .0. / 

DATA (0006 ( I ONE ) « I 0Nr = 2 * 9 ) / 

7.04OOOE—OJ * I .OOOOrE+OO .-0. , 4.90000E-02. 2 . 67200 E+ 00 . 

2 • I 6000E-0 l * 3.59000E+-00 . 2.70000E-01/ 

DATA (D007 ( I ONE ) . I ONc*2 *9 ) / 

-2 .6 OOOOE-O 1 • 0. , 3.1 2000E+00 • 2.60000E-01. 3. 1 2000E+00 . 

2 • 60 OOOE - 0 1 * 0* .0. / 

DATA (DOO0 f ! ONE ) t I ONr-2 .9 ) / 

7 . 84 OOOE— 0 1 * 1 .OOOOoE+OO * - 0 . 1 4 « 90000E-02 * 2.67200E + 00* 

2.t60f0r-01* 3.5900PE+00 . ? . 70000E-0 1 / 

DATA ( D009 ( I ONE ) * I ONf*2 » 9 ) / 

-2.60000E-0! . 0. , 3.12000E+00. 2.60000E-01. 3. 1 2000E+0 0 . 

2.6O000E-01 .0* .0, / 

DATA ( D 0 1 0 ( I ONF ) . I ONE = 2 ♦ 9 l / 

I .6O00OE4-00 * 2.2200PE + 00 .-0. , 1.00000E-01. 8.32000E+00. 

6.20000E-01 . 1 • 35000E+00 * 2 • 1 OOOOE-Ol / 
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DATA (DO] 1 < I ONE ) • I0Ne = 2*9]/ 

1 4*7] OOOE — 0 1 ♦ 5 • 1 OOOr>E-0 1,-0* • 4 . 7 1 OOOE -0 1 . - 4 • 32 OOOE-O 1 * 

1 3.90000E-02 *-3.30000E-0 1 , 9.00000E-02/ 

C 

DATA ( DO 1 3 < I ONE ) * I ONr *2 • 9 ) / 

1 1 • 40000E+0 1 * 3.4000OE+0 1 , -0* ♦ 7.00000E+00* 2.60000E+01* 

1 2 • OOOOOE+ 0 1 * 3 • ?000np + 02 * 1 *18000E + P2/ 

C 

DATA (DOl 4 ( I ONE ) . IONP = 2,9 )/ 

1 3 • 3000 OE + O 1 « l .OOOOnE+02 , -0* . 3.30000E+01* 3*40000E+01* 

1 6*7f>OOOE + Ol « 1 • 0000nE+02 * 1.00000E+02/ 

r 

DATA ( DO 1 5 ( I ONE ) « t ONr= 2 « 9 ) / 

1 1 • OOOOOE+0 1 « 5*0000r>E + 01 *-0* « 1 • OOOOOE+O 1 * l.OOOOOE+OI* 

1 2 • OAOOOE + O 1 « 0. #0. / 

C 

DATA (001 6 ( I TWO ) ♦ ITW0*2* 18 )/ 

1 5 • Oft 050E-0 1 « 5* 0605^E — 0 I « 5.06050E-01* 3,49000E-02« 0* • 

10# 4 0* ,0* , 3 • EOOOOE+O 1 « 1 • 1 OOOOE+C 3 ♦ 

1 1 « 1 OOOOF+O.7,- 1 • 1 0009E+C2, 1«10000F+02* 0. . 0. ♦ 

10 * « 0 • / 

c 

DATA ( DO 1 7 ( I TWO ) * 1 TW0*2 • 1 8 )/ 

1 5 • 06050E-0 1 « 5 • 06050E-0 l * 5.06050E-01* 3.49000E-02* 0* ♦ 

10* « 0. * 0* * 3*5 OOOOE+C l * 1#10000E+C3* 

1 1 • 1 0000E + 03*- l • 1 OOOoE + 02 , 1 *1 0000E+02 « 0* ,0. * 

10 . ♦ 0 • / 

r 

DATA ( DO I 8 ( ITHQEE)* 1 THREE *2« 27 ) / 

1 -9.20000E-02 * 8.4000OE-02. 3.53000E-OU 2.61000E-01, 1.80000E-01* 

1 8#0OOOOE-O2 . 9 • 6000nE “0 ! « 3.48000E-01* 7.00000E+00* 2.00000F+01* 
t 2 • OOOOOE+O 1 * 2 • OOOOnE+O 1 « 2.00000E+01* 0. ,0. , 

10* • 0* ♦ 3 • OOOOOE+OO « 3 • 60000E+02 * 3*60000E+C2« 

1 — 0 * • 1 • 2 0 0 0 o E + 0 2 * 0 * * 0 • * 0 » * 

10* / 

c 

END 
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Function Data Inputs 


-.020305 
-.69404 
.73487 
.53809 
- 1 .0400 
- 1 *3121 
- 1 .3045 
l .0900 
.23077 

6.0000 

.071 154 
3.9221 
- 0.1 8846 
- 8.2750 
0.0000 
-.57121 
-.59000 
.37654 
•23846 
7.9973 
-.00500 
.29918 
.61144 
-.50846 
•05000 
I .3058 
0.0000 
-.93000 
-.76479 
- . 1 4240 
0.0000 
-1 .354 0 
- 1 . 229 l 
- 1 .3530 
- 0.7651 1 
t .0 
t .0 
1 .0 
.622 
.082 
.008 
.01 3 

• 563 

• 622 
.062 
.000 

• 01 3 
.563 
.622 
.062 
.008 

• 0 1 3 
.564 

• 623 

• 1 12 
.008 

• 044 
. 565 
.626 

• 21 2 

• ooe 

• 1 44 

.567 

.632 

• 312 

• ooe 

• 244 
.571 
.640 
.340 
.008 

• 340 

• 570 

• 653 

• 340 
.041 
.340 
.588 
.672 
.340 

• 1 41 
. 340 
.602 
.699 
.340 
. 340 

• 340 

• 624 


-.26538 
-.66849 
.85767 
.34481 
- I .0540 
- 1 • 3180 
- 1 .3023 
1.3200 
3.2308 

6.0000 
*99615 
3.8179 
- 2.6385 
- 7.0627 
-.023333 
-.99000 
-.59000 
. 59000 
3. 3395 
7.6468 
.0066667 
.34846 
.65528 
-.70000 
. 70000 
1 .4530 
-.24333 
-.93000 
-.73239 

- 1.4100 
- 1 .3300 
- 1 • 2482 
- 1 .3762 
- 0.52766 
.7 
. 4 
.9 

.563 
. 0 ! 3 
.008 
.062 
,622 
.563 
.01 3 
. 008 
.062 
.622 
.564 
*01 3 
.008 
.062 
. 622 
.565 
.044 
.008 
.112 
. 62 ? 

.567 
. 1 44 
.008 
.21 2 
.626 
. 57 ] 

.244 
.008 
.312 
.632 
.578 
.340 
.006 
.340 
,640 
.588 
. 340 
. 04 1 
. 34 0 
.65 3 
.60 3 
,340 
.14 1 
. 340 
.672 
.624 
, 340 
, 340 
,340 
. 69*4 


-.54407 
-.62618 
.96585 
.17769 
-1 .0680 
- 1.3158 
- 1 • 3000 
1 , 320 o 

6.0000 

6.0000 
1 . 935 ? 
2 . 979 ] 
- 5.0296 
- 6.6058 
-.046667 
- * 59007 
-.59000 
.59000 
6.2590 
6.7574 
• 01 B 333 
.39229 
. 6991 l 
-.60000 
1 . 318 ^ 

1 . 545 ? 
-.48667 
-.92676 
-.70000 

- 1 . 4064 

- 1 • 3060 

- I . 2*74 
- 1 . 342 ? 
- 0 . 290 ? 1 
.34 

0 . 

• 1 25 
.509 
.01 I 
.008 

• 1 30 

.506 

• 0 l l 

,00 8 

. 1 30 

.500 
.01 1 
.008 
. 1 30 

.509 
.01 1 
,008 

• 1 SO 

.51 l 
. 076 
.008 
.280 

• 51 * 

• 176 
.008 

.340 

.51 9 
.276 
.008 
,340 

.527 

.340 

• 1 07 

.340 

.538 

.340 

.207 

.340 

.554 

.340 

.340 

.340 


-.72308 
-.42655 
.96040 
* 0! 269 
- 1.1010 
- I , 3! 35 
- 1 .0809 
1 .2007 

6.0000 
5,9793 
3.0327 
I .0379 
- 6.6927 
- 4.9015 
1 0439 
-.59000 
-.44234 
• 34484 
7.0238 
5.0648 
. 044098 
.43612 
.61087 
-. 1 3508 
1 . 5462 
1 .2948 
-.73000 
-.89437 
-.93000 

- 1 • 3980 
- 1 .2820 
- ] .2965 
- l .3081 
- .052766 
■ 1 6 
0 . 

O a 

.455 

.010 

• ooe 

• 361 

.455 

.01 0 

.008 

• 361 

.456 

.0 10 
.008 
.361 

.456 

.010 

.008 

• 3*1 

.457 

.010 

.008 

. 3*1 

.459 

• ioe 

.008 

• 3*2 

.463 

.208 

.074 

• 362 

• 46 e 

• 300 

• 1 74 

• 363 

.4 76 

• 340 

.274 
. 365 

.487 

• 340 

• 34 0 

. 36 T 


-.79846 
- . 1 2774 
.95042 

-l • 1 371 
-1.3113 
- 0 . 8! 88 
1 .0841 
6 ■ 0000 
4 .1*79 
3.Q529 
.79389 
-7.7777 
-2.431 3 
-.16705 
-.59000 
-.231 12 
.09968 
7.6349 
2.8779 
,0722 1 3 
.47995 
.36240 
.32984 
1 .4442 
. 8435 ! 
-.79250 
-.86197 
-.93000 

- 1 . 3896 
-1 .2580 
-1 . 3056 
-1 .274 1 

0 . 


.4 06 

.009 

.009 

.406 

• 406 

• 009 
.009 

.4 06 

.406 

.009 
.009 
.4 06 

.407 

.009 

.009 

.407 

• 407 

.009 

.009 

.407 

.409 
.04 1 

• 04 ] 
.409 

• 41 1 

• 141 

• 14 ] 

• 41 1 

.414 

.24 1 

.24 I 
.414 

.418 

.340 

.340 

• 418 

.425 

.340 

• 340 
.425 


-.75846 
. 1*600 
.84925 

- 1 .1836 
- 1 .3090 

- 0 , 4277 © 

6.0000 

,00 

3.9279 

.00 

- 8.2377 

.00 

-.26436 
- ,59000 
-.026556 


1 .2962 

C .0 

- . 85500 
- . 82958 
- .89000 

-1 .3812 
- 1.2340 
-1 .3247 
- 3 .2400 


• 361 
.008 
.0 10 

• 4*5 

.361 

.008 
.010 
.4 55 

. 361 
.008 
.010 

.4 56 

. 361 
.000 
.0 1 0 

.456 

• 361 
.008 
.0 10 

.457 

. 362 
.008 

• 1 08 
.4 89 

• 362 
.0 74 
. 208 
.463 

. 363 

. 1 74 

• 308 

. 4*8 

• 365 

.274 
. 340 

.4 76 

. 367 

• 340 

• 340 
.487 


-.71904 

•46371 

.74405 

- 1 .2479 
- 1 • 3068 
0.8200 

6.0000 

3.9971 

- 8.6650 

-. 4 l 779 
-.59000 
. 1 7508 

8. 1927 

.21082 

.56761 

-.27692 

1 .2030 

-.91750 
-. 7971 8 
-.51620 

- 1 .3728 
- 1.2100 
- l • 3439 
- 1 .0026 


. 1 30 
.008 
• 01 1 
.508 

• 1 30 
.008 

• on 

. 50 © 

• 1 30 
. 00 © 

• Oil 

.508 

• 100 
.008 

• on 

.509 

.280 
.000 
.076 
.51 1 

.340 

.000 

• 1 76 
*5 1 4 

.340 

.008 

.276 

.519 

• 340 

.107 

.340 

.52 7 

• 34 0 

.207 

.34 0 

.530 

.340 
.34 0 
.340 
.554 


8.0144 

.00 

. 1 2246 
. 52378 
•036667 
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-.995 

-.962 

-.926 

-.885 

-.84 1 

-.794 

-.840 

-.990 

- 1.120 

- 1.183 

-.967 

-.752 

-.537 

-.322 

1 07 

.107 

.322 

.537 

.752 

.967 

1.183 

1.120 

.980 

.840 

• 794 

.041 

.085 

.926 

.062 

,995 






-.995 

-.963 

-.926 

-•885 

-.84 1 

-.794 

-.822 

-.920 

- 1 .034 

- 1.140 

-.982 

-.764 

-.546 

-.327 

-.100 

, 1 09 

.327 

.546 

.764 

• 982 

! • 1 4 

1 .034 

.928 

.822 

.794 

• 04 1 

.805 

.926 

. 95 ? 

.995 






-.996 

-.963 

-.926 

-.805 

-.84 1 

-.794 

-.008 

-.865 

-.963 

- 1.04 

- 1 .009 

-.765 

-.560 

-.3 36 

-.112 

.112 

.336 

.560 

.785 

1 .009 

l .040 

.963 

.885 

.808 

. 794 

.84 l 

.805 

.926 

.963 

.996 






-.997 

-.964 

-.927 

-.886 

-.842 

-.794 

797 

-.851 

-.904 

-.958 

- 1.012 

-.817 

-.503 

-.350 

-.116 

. 116 

.350 

.583 

.817 

1 .012 

.958 

.904 

.851 

.797 

• 794 

• 042 

• 886 

.927 

.964 

,997 






- 1 .000 

-.966 

-.929 

-.807 

-.842 

- . 794 

-.707 

-.022 

-.857 

-.892 

-.927 

-.864 

-.617 

-.370 

-.123 

• 1 23 

.370 

.617 

.864 

.927 

• 092 

.857 

.822 

.787 

• 794 

• 642 

• 887 

*929 

.966 

1 .000 






- 1 .005 

-.970 

-.932 

-.890 

-.844 

- .794 

-.779 

-.798 

-.818 

-.837 

-.856 

-.875 

-.666 

-.40 1 

-.1 33 

.133 

.401 

.669 

.875 

.856 

#837 

.810 

,798 

.779 

.794 

.844 

.090 

• 932 

.970 

1 .005 






-1 .013 

-.977 

-.937 

-.893 

-.846 

- .795 

-.773 

-.779 

-.785 

-.792 

-. 790 

-.804 

-.749 

-.449 

-.149 

.149 

.449 

• 749 

.804 

• 790 

• 792 

• 705 

.779 

• 773 

.795 

.846 

.893 

. 93 T 

. 97 ^ 

1.0 13 






- 1.025 

-.987 

-.945 

-•899 

-.850 

-.796 

-.768 

-.764 

-.760 

-.756 

-.752 

-.740 

-.74 4 

-. 740 

-.283 

.283 

,74 0 

• 744 

.746 

.752 

.756 

.760 

.764 

.768 

• 796 

.050 

• 899 

• 94 0 

.987 

1 .025 






- 1 .042 

- 1.002 

-.957 

-•908 

-.855 

-.796 

-.762 

-.740 

-.733 

-.718 

-.704 

-.689 

-.675 

-.4 66 

-.155 

. 1 55 

.466 

• 675 

.689 

• 704 

.718 

.733 

,748 

.762 

♦ 798 

• e 55 

.900 

.957 

1 .002 

! .042 






-1 .067 

- 1 .023 

-.974 

-.920 

-.862 

-.801 

-.755 

-.727 

-.699 

-.6 70 

-.642 

-.614 

-.500 

-.300 

-. i on 

.100 

. 300 

.500 

.614 

• 642 

.670 

.699 

.727 

.755 

.80 I 

.862 

• 920 

.974 

1 .021 

1 .067 






-.3253 

-.0337 

.0039 

• 0635 

• 4818 



-.3538 

-.0555 

- .0026 

.0477 

.4556 



-.3789 

-.0720 

- .0094 

.0267 

.4264 



-.4091 

-.0651 

• 0452 

• 1 548 

• 6001 



-.4561 

-.1155 

-.0056 

• 1 045 

.5627 



-.5027 

-• 1658 

-.0562 

• 054 1 

.5168 



-.5328 

-.0933 

.0527 

• ! 985 

• 7029 



-.591 9 

-.1532 

-.0074 

. 1 386 

.724 1 



-.6510 

-•2131 

-.0672 

• 0767 

.6657 



-.6802 

1 289 

.0547 

.2302 

• 9725 



-. 7436 

-• 1925 

-.0092 

.1744 

.9091 



-.8072 

-.2566 

-•0730 

• 1105 

.8458 




113 


Program Setup 


The program deck is read into the computer after selecting the RESET mode and 
both RELEASE switches on the program control console (fig. 1(b)). When real-time sta- 
tus is achieved, the READ mode is selected to store variables for printed output followed 
by depression of the lower RELEASE switch. The check cases are then run before further 
use of the helicopter simulation program. 

The static check is initiated by depressing Function Sense Switch 6 (LDISI(38)), fol- 
lowed by selection of the READ mode and the lower RELEASE switch (fig. 1(b)). With an 
untrimmed (148.16 km/hr (80 knots)) condition selected, the simulation is momentarily 
place( in the HOLD mode and then in the OPERATE mode for approximately 1 second. 

Data are then printed by selecting the PRINT mode followed by depression of the lower 
RELEASE switch. The values at time equal zero (t = 0) are used for comparison with 
those listed in table I. With successful completion of the static check, the analyst is 
ready to conduct the dynamic check. 

Before executing the dynamic check, the strip chart recorders should be set up as 
detailed on the overlay plots for each control perturbation selected. (See figs. 2,3, and 
4.) The dynamic response of the simulation is verified by imposing a doublet perturba- 
tion equivalent to a 2.54-centimeter (1 inch) stick deflection on each of three control 
inputs, that is, A lg , B 1Sc , and 

Function Sense Switch 1 (LDISI(33)) is depressed and the desired control deflection 
is input by entering the value in the appropriate TABLE location. The values equivalent 
to a 2.54-centimeter (1 inch) stick deflection are as follows: 

A1SSTEP (TABLE(60)) * 0.032 radian for a perturbation on A^g 

C 

B1SSTEP (TABLE(61)) = 0.049 radian for a perturbation on B^g^, 

THTRSTP (TABLE (62)) = 0.087 radian for a perturbation on 

The simulation is placed in the HOLD mode for a few seconds and then in the 
OPERATE mode for approximately 20 seconds for each of the three control perturba- 
tions. Output of the strip chart recorders can then be compared with the appropriate 
independent check overlay plots. (See figs. 2,3, and 4.) Strip chart recorders are then 
reset for normal recording of data. With successful completion of the dynamic check, 
the analyst is ready to conduct the cockpit checkout. If no cockpit is desired, the analyst 
is ready to conduct test cases from the program control console. 

Cockpit Checkout 

After the static and dynamic checks have been executed, the simulation is now ready 
for the cockpit checkout. Telephone communication between the cockpit and the program 
control console is established and Function Sense Switches 12 and 13 (LDISI(44) and 
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LDISI(4 5)) are depressed to activate inputs to the computer from the cockpit. Function 
Sense Switch 10 (LDISI(42)) is depressed for several minutes to allow the cockpit engineer 
to statically check his instruments. During this time, voltages on all cockpit instruments 
(DAC(25)-DAC(40)) are checked on the digital decimal display unit for comparison with 
those listed in table n. When the cockpit engineer is ready, Function Sense Switch 10 is 
released and voltages for full deflections of the controls are checked on the digital decimal 
display unit for comparison with those listed in table HI. After the cockpit checkout has 
been successfully completed, the simulation is ready for piloted test flights. 

Console Run Procedure 

Upon completion of the static and dynamic checks, the simulation program is ready 
for the analyst to set up preprogramed flights through the use of potentiometers repre- 
senting the control sticks and pedals. The procedure is as follows: 

Step 1: Select the RESET mode. All other switches should be in release positions. 

Step 2: Select a value for KNOTS (INTEG(6)) from 1 to 6 corresponding to the preset 

conditions closest to the desired forward velocity in order to initialize the trim 
circuit. 

Step 3: Input the desired forward velocity (TABLE (50)) and altitude (TABLE (52)). 

Step 4: Activate the trim circuit by depressing Function Sense Switch 5 (LDISI(37)) until 
White Light 1 (LDISO(61)) illuminates. Then release Function Sense Switch 5. 

Step 5: Set NSTICKO (INTEG(l)) equal to 1; this allows control stick inputs from the 
potentiometers located on the program control console. 

Step 6: Depress Function Sense Switch 13 (LDISI(45)) to activate subroutine ADCIN; this 
allows control stick inputs to the main program. 

Step 7: Set potentiometers 1, 2, 3, and 4 until White Lights 21 through 24 (LDISO(81)- 

LDISO(84)) illuminate, respectively; this indicates that the controls are trimmed. 

Step 8: Select one of the following variables and input the desired stick displacement in 
inches to program the desired flight condition: 

XI (TABLE(84)) - collective stick 
X2 (TABLE(85)) - lateral stick 
X3 (TABLE (86)) - longitudinal stick 
X4 (TABLE (87)) - tail rotor pedals 

Then set the time of input T1 (TABLE(88)) and the time of removal of input 
T2 (TABLE(89)). 

Step 9: The SAS system may be activated, if desired, by depressing Function Sense 
Switch 9 (LDISI(41)). 
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Step 10: The simulation is now ready for the analyst to conduct a test run. The HOLD 
mode is depressed momentarily and is followed by the OPERATE mode for the 
desired number of seconds. At the completion of the run, the analyst depresses 
the RESET switch which returns the computer to the conditions at the beginning 
of the run. Desired printout may be obtained at this point by depressing the 
PRINT switch and the lower RELEASE switch. For sample strip chart output, 
see figure 5. 


Cockpit Run Procedure 

Upon completion of all check procedures, the analyst is ready to set up the initial 

conditions for the desired piloted flight cases. The procedure is as follows: 

Step 1: Select the RESET mode. All other switches should be in release positions. 

Step 2: Select a value for KNOTS (INTEG(6)) from 1 to 6 corresponding to the preset 

conditions closest to the desired forward velocity in order to initialize the trim 
circuit. 

Step 3: Input the desired forward velocity (TABLE(50)) and altitude (TABLE (52)). 

Step 4: Activate the trim circuit by depressing Function Sense Switch 5 (LDISI(37)) until 
White Light 1 (LDISO(61)) illuminates. Then release Fmction Sense Switch 5. 

Step 5: Depress Function Sense Switch 13 (LDISI(45)) to activate subroutine ADCIN; this 
allows control stick inputs to the main program. 

Step 6: Depress Function Sense Switch 12 (LDISI(44)) to allow cockpit control of the stick 
inputs . 

Step 7: Pilot trims control sticks and pedals to program computed values. Cockpit trim 

is reached when White Lights 21 through 24 (LDISO(81)-LDISO(84)) are illuminated. 

Step 8: Depress Function Sense Switch 11 (LDISI(43)) to activate cockpit control of com- 
puter control modes, that is, RESET, HOLD, OPERATE. 

Step 9: The SAS system may be activated, if desired, by depressing Function Sense 

Switch 9 (LDISI(41)) or by having the pilot depress the SAS switch in the cockpit. 

Step 10: The simulation is now ready for flight runs with the pilot depressing the HOLD 
switch and then the OPERATE switch. At the completion of a flight, the pilot 
depresses the RESET switch which returns the computer to the conditions at the 
beginning of the flight. Desired printout may be obtained at this point by depress- 
ing the PRINT switch and then the lower RELEASE switch on the program control 
console. 
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PROGRAM VALIDATION 


This section presents a comparison of the Langley simulation program (Cl 152) with 
the simulation program of reference 1. Figure 6 presents data X AO g, X^-g, Y cg , 

X T r, and 9 from the Langley simulation program (Cl 152), the simulation program of 
reference 1, and an unpublished program based on reference 1 provided by the U.S. Army 
Electronics Command, Fort Monmouth, New Jersey. The good match between the Langley 
data and the reference 1 updated (unpublished) program data validates the implementation 
of the simulation program onto the Langley RTS system. 

CONCLUDING REMARKS 

A man-in-the-loop real-time helicopter simulation program has been developed and 
integrated with a suitable fixed-base cockpit by using the Langley Research Center real- 
time simulation system. With minor changes to the mathematical model and detailed air- 
craft data, this program can be used to simulate a variety of single-rotor helicopters. 

This fact has been substantiated by a program presently being used in support of two 
piloted studies of the Sikorsky S-61 helicopter, a commercial passenger-type helicopter. 
One study concerns pilot workload, route structures, and Instrument Flight Rules (EFR) 
procedures in congested areas. The second study concerns development and evaluation 
of computer-generated pilot displays and instruments for approach and landing. In both 
studies, the pilots' comments are favorable; that is, the simulation more than adequately 
represents a general single-rotor helicopter. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va., February 27, 1974. 


117 


APPENDIX A 


AIRCRAFT EQUATIONS 

This appendix contains the equations used in program HELIC. The equations are 
listed according to the subroutines in which they are used along with references for their 
source. 

Control equations: 

Subroutine ADC IN (ref. 1) 

A os c = ( 85 ° + °- 92x aos) rpd 
A 1S C = (- 10 - 8 ° + 1.85Y cs )rPD 
B 1S C = (- 132 ° + 2.81X(-.g^RPD 

0 TRc = ( 23 -°° ' 5 - 0 X Tr) RPD 

where RPD is the factor for converting degrees to radians. 

General aerodynamic variable equations: 

Subroutine AEROVAR (refs. 1 and 2) 

_ \/u 2 + v 2 + w 2 
A SPD- 1.689 


V T =(/u 2 + v2 +(w - W IM )' 


ADVPj = 


2ttR 


B 


ADVP 2 = i 


W IM = W IM + 


J MR 


27rR B pV T 


2 „ " W IMj 


At 

T 
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q L = 2 P ( u2 + v2 ) 

?v - |p[“ 2 ♦ (w - W IM )2 

Body derivative equations: 

Subroutine BODYDER (refs. 1 and 2) 

. _ L A + ( I YY ~ 1 Zz)q r 



r xx 


m a + 

( r ZZ ■ 

I xx) r P 


i yy 


n a + i 

( : XX ■ 



r zz 


r cos 

4> + q 

sin $ 


cos 9 


q cos 

$ - r sin $ 

p + 'I' 

sin 9 



• X A 

u = — + rv - qw - g sin 0 


v = + pw - ru + g cos 9 sin <f> 


• Z A 

w = -jjj- + qu - pv + g cos 9 cos $ 

Earth orientation equations: 

Subroutine EARTHM (refs. 1 and 2) 

N = |u cos 9 + (v sin <i> + w cos $)sin fljcos - (v cos $ - w sin <|>)sin 
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E = |u cos 9 + (v sin <i> + w cos <f>)sin 0]sin ^ + (v cos <i> - w sin 3>)cos 
h = u sin 9 - (v sin $ + w cos 3>)cos 8 


Engine equations: 

Subroutine ENGINE (refs. 1 and 8) 

Q E D = ' fi ) + Q MR 

Q E = Q E + ( q e d ' q e)^ 

_ 202.63158(h - 2500) + 412 500 
^max ~ ft 

A Q E ' q mr 

11 ~ T 

a rot 

12 = 12 + riJAt 


S H p - 


Qmr° 

550 


= 


g rwt 

p(7rR B 2 )(fiR B ) 2 


Cp = 


550S 


UUUOjrp 

P^ R B 2 )( nR B f 


C Q = 


Q 


MR 


pR B ( 7rR B 2 )( nR B ) 


M = 


1.689A SpD 
12 Rq 


m tip _ 


1.689Agpjj + 12Rg 

Vs 
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Force and moment equations: 

Subroutine FANDM (refs. 1 and 2) 


X A 

= x F + x w + f xr 




y a 

= Y f + Y tr + Y vs + F yr 




Z A 

= Z F + z w + Z HS - l mr 




l a 

= l F + + d z F YR + ^d z 

- d z 
^TR 

) y tr + ( d z - d z vs ) Y vs 

+ I R f2q 

m a 

- M F + - d z F XR + d x L MR - J 

d z ■ d Z w ) X W + ( d x w - 

d x) z w 


+ ( dx HS " d x) Z HS “ J R fi P 




n A : 

= N F + d x F yR - (d XTR . d x 

) Y TR ' 

( d x vs ■ d x) Y VS + Qe 




Fuselage equations: 

Subroutine FUS (refs. 1 and 2) 

a p = tan-l( W 1 Vf W 1m| 

l f = q L % c z Fi + up % 2 

M F = ^v v m F c m F + u QC m 

1 r 2 
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n F = ^L V n F C n Fi + urCn F 2 
X F = ^l s x Fi c x 1 + ^v s x F2 c x 2 

y f = q L s Y F c Y F 
z f = q v s z F c z F 
Z W =q v s z w c z w 
x w = q v s x w c x w + % z w 

^IWM “ X W z w 


Main rotor equations: 

Subroutine MAINROT (refs. 1 and 2) 

K A = f" 

A 0 l B 


1 


k ab, = r 

1 l B 


K 


ab 2 


1 + nigRge 

% 


K 


m g Rg e f2 


ABg 21 


B 


F MRS = \ N B m B R B e 


= a 0SS ' 01 1SS cos *R %SS sin ^R 

0* = («1SS sin ^R-^lSS cos 

W^ Y = W M W Fy [l.O + W IFy (u cos * R 

Ut ^Y = Y PE fi + u sin * R + v cos * R 


- v sin 
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% Y = Y PE( q cos *R + p sin *r) ' Y/% + w - /3^u cos * R - v sin * R ) - 

e *Y = Y TW + A OS “ A 1S cos ^R - B 1S sin % - (“OSS ' °- 043 ) 

U P, 


$ 




*Y U T 




" Ay + ^Y 


J PDR = Kt Urp 

PL,1 Vy L y L x^y t ^y 


Dpnp “ Kn Cn U rp - Kht Cj U-n Ut 

pdr ^y u y u ^y t ^y dl y l * y p ^y t *y 


M 


PDR *Y Km Y C U^ Y U Tv^Y 


Q 


PDR vJ-Y ^Y^Y^Y K Q l y Cl ^Y Up ^Y Ut 


^Y 


0! 


088 " u^azA 




M 


PDR^y 


“1SS 


- p 


N 


AZ 


' K AB £ 

~~a) Z Z( Mp DR^ Y Sin *r) " K AB 2 q + k ab 3 Ass 


0 


/K AD \ % Y 

= V AB il 


iss ^n az )\ a 



M PDR^ y C0S *r) - K AB 2 P - K AB 3 a iss 


“iss _ a iss + a iss At 


Ass _ Ass + Ass At 


* 


/pNx 


R Y 


F XR In A zJ I J (^UpDR^Y C0S *R ‘ Dp DR* Y Sin *r) 
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Y 

F YR - (§|) I l (-^ l pdr„ s1 " * R - d pdr* y “ s *h 



sin * R ) + Fmrs^ 2 ^ISS 


'I't 


R Y 

1.0 + G^uGiryj ^ ^L Y 


/PN b 

J MR = (N^Ji i - u + u eh u ev; A a^pdr, 


1<Y 


* 


R Y 


^MR=sn Q PD% 


Tail equations: 

Subroutine TAIL (refs. 1 and 2) 

V TR = v - ( d X TR ' d x) r 

W HS " w + ( d X Hg ' d x) q _ °W W IM ' W IWM 
6 e » 10.6852° - 2.0405X CS + 0.2445X Eg 

a HS = tan 1 ( — u^ ) + ^TR^E + *HS 

^vs = tan ’\“u B ) + ! VS 

Zttc = C y 

HS V ^HS ^HS 

Y. rc = q T S Y C Y 

VS H L Yyg Yyg 

y tr = K TR 2 C Y tr 
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INDEXING TECHNIQUE FOR FUNCTION GENERATION 

By Lawrence E. Barker, Jr., and Kemper S. Kibler 
Langley Research Center 

This appendix gives the mathematical and programing details of the technique used 
at Langley Research Center for the generation of aerodynamic functions for real-time 
simulation (ref. 5). The FORTRAN code for the subroutines is given followed by a 
description of the COMPASS version of this technique. 

Since function generation represents the largest portion of the problem central 
processing unit (CPU) time in real-time digital programs with a large amount of aero- 
dynamic data, prime attention was given to developing a technique to reduce this time. 
The mathematical description upon which this function method is based can be illus- 
trated by the following sketches for an arbitrary function of two variables: 



Sketch (a) - Three-dimensional view. Sketch (b) - Bottom view. 


Contour lines where one of the variables is held constant are drawn, such as f(x,y^ and 
f(x,y 2 ). It is assumed that the values f^x^.y^ corresponding to preselected grid points 
are stored in some array and that the area of interest is the area marked with an asterisk 
The interpolated value of the function f(x,y) is obtained by performing a linear interpo- 
lation between f^x,y^ and f(x>y 2 ) where 


f(x,yi) = f(x 1>yi ) + ^-lL|f(x 2 ,y 1 ) - f(x 1>yi j| 


(Bl) 
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«(*.y 2 ) - f ( x i’ y 2) + ^ll? 512 ’ 5 ' 2 ) ' t ( Xl ’ y2 I <B2) 

This interpolation can then be written 

f(x,y) = f(x,y,) + - f ( x > y l3 (B3) 

Equations (Bl), (B2), and (B3) are sufficient to give the value for f(x,y) in the particular 
area. For these equations to be of use they must be implemented to change as the inde- 
pendent variables move to a new area. 

The method is based on a setup assuming the function for equal intervals of the inde- 
pendent variables. This will include most aerodynamic functions. A more generalized 
technique would be to provide for several areas of equal intervals per independent varia- 
ble. The function technique described provides for three areas of equal intervals per 
independent variable. The areas and the corresponding interval for each area are 
defined by an array. The generated function value Is returned through the first word in 
this array. A FORTRAN CALL STATEMENT is used to pass the array of function val- 
ues, the array of function parameters, and the independent variables to the FUNC 
subroutine . 

The key statement in the coding is the division in fixed point mode of the real 
FORTRAN variable AX (or AY or AZ) by its equal increment. This division allows selec- 
tion of the proper interpolation area without a search routine. For this function genera- 
tion technique the execution time is independent of the number of curves per function or 
the number of points per curve. The function value of interest will be obtained by linear 
interpolation between the stored data values by equations similar to equations (Bl), (B2), 
and (B3). These equations are implemented in another form to be more efficient. The 
implemented FORTRAN equations for function of two variables are 

G1 = TX *FUNC (II , J 1) + SX *FUNC (12 ,J 1) 

G2 = TX *FUNC (II , J2) + SX*FUNC(I2,J2) 

RESULT = G1 + (AY/DELY - JS) *(G2 - Gl) 


(B4) 

(B5) 

(B6) 


where 

SX = AX/DELS - IS 

and 

TX = (1.0 - SX) 
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EXAMPLE : Consider the function of two variables defined in sketch (c). For this exam- 
ple, the independent variable X is divided into three areas. These areas are subdivided 
into equal segments. 



XAREA1 XAREA2 XAREA3 
Sketch (c) - Typical function of two variables. 


The tabular form of the foregoing function is as follows: 


XAREA1 
DELX1 =0.2 
BIASX1 = -0.2 


XAREA2 
DELX2 = 0.1 
BIASX2 = -0.4 


XAREA3 
DELX3 = 0.5 
BIASX3 = 2.8 


X Y 

0 

5 

10 

15 

0.2 

FXY01(1,1) 

FXY01(1,2) 

FXY01(1,3) 

FXY01(1,4) 

.4 

FXY0 1(2,1) 







.6 



" 






.7 









.8 









.9 









1.0 









1.1 









1.2 









1.7 









2.2 

FXY01(11,1) 

FXY01(11,2) 

FXY01(11,3) 

FXY01(11,4) 
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where 

XAREA = (XARE A1 ,XARE A2 ,XARE A3) is the value of the independent variable X at the 
beginning of an area (XAREA point) 

DELX is the desired value of increment which begins at the XAREA point 
NUMPTS is the number of data points previous to the XAREA point 

The BIAS(BIASX1 ,BIASX2 ,BIASX3) are functions of the independent variable X and 
the parameters DELX1, DELX2, and DELX3. The formula for determining the BIAS is 

BIAS = (DELX) (NUMPTS) - XAREA 

For the example shown, 

BIASXl = (0.2)(0) - 0.2 = -0.2 

BIASX2 = (0.1) (2) - 0.6 = -0.4 

BIASX3 = (0.5) (8) - 1.2 = 2.8 

The BIAS for changes in DELY or DELZ are calculated in a similar manner. The 
NUMPTS in the BIAS formula would be replaced by the number of function curves pre- 
vious to the YAREA point and the number of families previous to the ZAREA point, 
respectively. 

For a two- variable function the parameter array (locations 2-18) is 

PXY0 1 (XA RE A2 ,X ARE A3 , BLASX 1 ,DE LX 1 , BIASX2 ,DE LX2 , BLASX3 ,DE LX3 ,XPOINT S , 
YAREA2,YAREA3,BIASY1,DELY1,BLASY2,DE LY2,BLASY3,DELY3) 


Substituting numerical values gives 

DATA(PXY01(ITWO),ITWO=2,18)/.6,1.2,-.2,.2,-.4,.l,2.8,.5,ll.,15.,15.,0.,5.,0.,0.,0.,0.1 

The FORTRAN CALL STATEMENT is 

CALL FUNC2(FXY01,PXY01,X,Y) 


CYP = PXY01(1) 
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C * ONE VARIABLE 

C * THE FOLLOWING IS COOED IN COMPASS. THE 

C * CODE IS SHOWN AS AN EXAMPLE OF WHAT IS 

C 

I F ( X I NOEP-XA RF A 2 ) 1 ,1 , ? 

1 OFLX = DEL X 1 

AX = XINDEP + RIASXl 
GO TO S 

2 IF (X INDEP-XAREA3 ) 3 * 3 , A 

3 DEL. X = DELX2 

AX = X I NOE P + HI A S X 2 
GO TO 5 

A DEL X - DELX3 

AX = XINDEP + BIASX3 

5 I S = A X/DFLX 

SX - AX/DFLX - IS 

RESULT = < 1 #-SX )*FUNC US+1 ) + SX*FUNC < I S +2 ) 


C ** TWO VARIABLF 

c * THE FOLLOWING IS CODED IN COMPASS* THE 

C * CODE IS SHOWN AS AN EXAMPLE OF WHAT IS 

C 

I F ( X I NDEP-XAREA2 ) 1*1*2 

1 DEL X = DEL XI 

AX = XINDEP + BIASXi 
GO TO 5 

2 I F ( X l NDEP-XA RF A3 > 3*3*4 

3 DELX = DELX2 

AX = XINDEP + B I ASX? 

GO TO S 

4 DELX = DFLX3 

AX = XINDFP + PIASX3 
8 IF ( Y INDEP-YAREA2 ) G* 6*7 

6 DELY = DEL Vi 

AY = Y I NOE P + BIASYI 
GO TO 13 

7 I E ( Y 1 NDEP— YAWF A3 > 6*8*9 


DELY 

- 

DEL Y2 


ay 

= 

Y I NDEP + 

B I ASY2 

GO TO 

1 3 


DELY 

= 

DE L Y 3 


AY 

= 

Y I NDEP + 

B I ASY3 

I S 

= 

AX/DELX 


JS 

= 

A Y/DFLY 


1 1 

- 

I S + 1 


J1 

= 

JS + 1 


12 

= 

IS + 2 


J? 

= 

JS + 2 


SX 

= 

AX/DFLX 

- IS 

TX 

= 

X 

LA 

1 

• 


G 1 

= 

T X * F U N C ( 

I 1 ♦ J1 ) + SX*FUNC ( 12* J1 > 

G2 

= 

TX*FUNC ( 

1 I * J2 ) + SX*FUNC ( r 2 , J2 ) 

RESULT 

= G1 + ( AY/DFLY-JS)* (G2-G1 ) 


p ortran 

CODED* 


fortran 

CODED* 
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C * THREE VARIABLE 

C » THE FOLLOWING IS CODED IN COMPASS. T HF FORTRAN 

C * CODE IS SHOWN AS AN EXAMPLE OF WHAT IS CODED, 

r 

I F ( X I NDEP-X ARE A2 ) 1 .1 *2 

1 DEL X = DEL X 1 

AX = XINOEP 4 BIASX1 

GO TO H 

2 I F ( X l NJtP- XA PE A 3 ) 3*3*4 

3 DEL X = DELX2 

AX = XINOEP 4 BIA5X2 

GO TO 5 

4 DELX = DEL X 3 

AX = XINDEP + HIASX3 

B IF ( Y INDEP-YAWFA2 ) 6*6*7 

6 DELY = DEL Y 1 

AY = VINDFP 4 BIASYl 

C-0 TO 13 

7 IF(YINL)EP-YAMEA^) 8 * B * 9 


8 

DELY 

= 

DEL Y2 




ay 

= 

YINDEP 4 B I A S Y2 




GO TO 


1 3 



9 

DELY 

- 

DELY 3 




A Y 

- 

YINDEP 4 BIASY3 



1 3 

I F ( Z I NDEP-ZAHEA2 ) 14,14* 

1 5 


14 

DELZ 

- 

DELZ 1 




AZ 

= 

ZINOFP 4 RIAS71 




GO TO 


1 8 



1 5 

1 F ( Z 1 NDEP-ZAREA3 ) 16*16* 

I 7 


1 6 

DELZ 

= 

DEL Z2 




AZ 

= 

Z1NDEP 4 BIASZ2 




GO TO 


1 8 



1 7 

DELZ 

= 

DEL Z 3 




AZ 

= 

ZINDEP 4 BIASZ3 



1 8 

I S 

= 

AX/DELX 




JS 

= 

A Y/DELY 




KS 

= 

AZ/DELZ 




I 1 

= 

I S 4 1 




J1 

= 

JS 4 1 




K 1 

= 

KS 4 1 




I 2 

- 

IS 4 2 




J2 

= 

JS 4 2 




K2 

= 

KS 4 2 




SX 

= 

A x/DFLX - IS 




TX 


1 . - SX 




5 Y 

= 

AY/DELY - JS 




G 1 1 

- 

TX*FUNC ( 11 * J 1 * K l ) 

4 

SX + FUNC ( I 2 * J 1 * K 1 ) 


G2 1 

= 

TX*FUNC ( I 1 * J2 * K 1 ) 

4 

SX*FUNC ( I 2 * J2 * K 1 ) 


G1 2 

= 

TX*FUNC ( 1 1 * Jl * K 2 ) 

4 

SX* FUNIC ( 12* Jl *<2 ) 


G 22 

= 

TX*FUNC (11 * J? *K2 ) 

4 

SX*FUNC ( I 2 * J2 * K? > 


HI 

= 

Gil v S Y * (G2 1 -Gil 1 




H2 

= 

G 1 2 4 SY*<G22-G12) 




RESULT 

= HI 4 ( AZ/DELZ — KS ) 

* (H2-H1 ) 
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SUBROUTINES FUNC1, FUNC2, FUNC3 


Language : COMPASS 

Purpose : To provide functions of one, two, and three variables. The method for gener- 
ating these functions which is described in the reference is basic for an understanding 
of the material presented below. 

Use : CALL FUNC1(FUNC, RESULT ,XINDEP) for functions of one variable 

CALL FUNC 2 (FUNC, RESULT ,XINDEP,YINDEP) for functions of two variables 
CALL FUNC3(FUNC, RESULT ,XINDEP,YINDEP,ZINDEP) for functions of three 
variables 

FUNC the array which contains the function values. 

RESULT the generated function value which is returned through the first word of 

the array. The remaining locations in the array must contain the func- 
tion parameters (see below). 

XINDEP,YINDEP,ZENDEP the X,Y,Z independent variables for the function. 

For functions of one variable, the RESULT array must be dimensioned 
for nine locations and locations 2 through 9 must contain the following 
parameters (in order): 

XAREA2 ,XARE A3 ,BIASX 1 ,DE LX1 ,BIASX2 ,DELX2 ,BIASX3 ,DE LX3 

For functions of two variables, the RESULT array must be dimensioned 
for 18 locations. Locations 2 through 9 are the same as those for the 
function of one variable. Locations 10 through 18 must contain the 
following parameters: 

XP0INTS,YAREA2,YAREA3,BIASY1,DELY1,BIASY2,DELY2,BIASY3,DELY3 

For functions of three variables, the RESULT array must be dimensioned 
for 27 locations. Locations 2 through 18 are the same as those for the 
function of two variables. Locations 19 through 27 contain the following 
parameters: 

YCURVES,ZAREA2,ZAREA3,BIASZ1,DELZ1,BIASZ2,DELZ2,BIASZ3,DELZ3 
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where 

AREA - that is, XAREA1,XAREA2,...,ZAREA2 are the values of the inde- 
pendent variables at which point the function values are tabulated 
at a different incremental value than the previous area. 

DEL - that is, DELX1,DELX2,...,DELZ3 are the incremental values of the 
independent variable at which the functions are tabulated. 

BIAS - that is, BIASX1,BIASX2,...,BLASZ3 are the values that when added 
on to the independent variable enable the function to be treated as 
an equal increment function. 

XPOINTS - the number of points per function curve. 

YCURVES - the number of function curves per family. 

Restrictions: Limits should be placed on the X, Y, and Z upper and lower bounds of the 

function to avoid calling an address out of range. 

Method: See reference. 

Accuracy: Of the order of linear interpolation. 

Reference: Indexing Technique for Function Generation. (See pp. 125-130 of appendix B.) 


Storage : FUNC1 - 13 g locations. 

FUNC2 - 26g locations. 

FUNC3 - 52g locations. 

Subprograms used : None. 

Other coding information : The timing requirements are as follows: 

FUNC1 - 20 microseconds. 

FUNC2 - 40 microseconds. 

FUNC3 - 60 microseconds. 

Source: See reference. 

Author of subroutine: Dave E. Eckhardt, Jr. 

Subroutine date: December 15, 1970. 
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FUNCTION DATA 


This appendix contains the data input to the function generation scheme described 
in appendix B. 


Angle * 
breakpoint 

Cm F 1 ( a r) 

c z f ( Q! f) 


C n Fl (%) 

c y f (%) 

Cz w( Qf w) 

-3.12 

-0.020385 

0.23077 

0.071154 

-0.18846 

0.23846 

0.0500 

-2.86 

-.28538 

3.2308 

.99615 

-2.6385 

3.3385 

.7000 

-2.60 

-.54407 

6.0000 

1.9352 

-5.0296 

6.2593 

1.3185 

-2.34 

-.72308 

6.0000 

3.0327 

-6.6927 

7.0238 

1.5462 

-2.08 

-.79846 

6.0000 

3.8529 

-7.7777 

7.6349 

1.4442 

-1.82 

-.75846 

6.0000 

3.9279 

-8.2377 

8.0144 

1.2962 

-1.56 

-.71904 

6.0000 

3.9971 

-8.6650 

8.1927 

1.2038 

-1.30 

-.69404 

6.0000 

3.9221 

-8.2750 

7.9973 

1.3058 

-1.04 

-.66849 

6.0000 

3.8179 

-7.8627 

7.6468 

1.4538 

1 

00 

-.62618 

6.0000 

2.9781 

-6.6058 

6.7674 

1.5452 

-.52 

-.42655 

5.9793 

1.8379 

-4.9015 

5.0648 

1.2948 

-.26 

-.12774 

4.1679 

.79389 

-2.4313 

2.8779 

.84351 

.00 

.16600 

.0000 

.0000 

.0000 

.0000 

.0000 

.26 

.46371 

-4.1679 

-.79389 

2.4313 

-2.8779 

-.84351 

.52 

.73487 

-5.9793 

-1.8379 

4.9015 

-5.0648 

-1.2948 

.78 

.85767 

-6.0000 

-2.9781 

6.6058 

-6.7674 

-1.5452 

1.04 

.96585 

-6.0000 

-3.8179 

7.8627 

-7.6468 

-1.4538 

1.30 

.96040 

-6.0000 

-3.9221 

8.2750 

-7.9973 

-1.3058 

1.56 

.95042 

-6.0000 

-3.9971 

8.6650 

-8.1927 

-1.2038 

1.82 

.84925 

-6.0000 

-3.9279 

8.2377 

-8.0144 

-1.2962 

2.08 

.74405 

-6.0000 

-3.8529 

7.7777 

-7.6349 

-1.4442 

2.34 

.53808 

-6.0000 

-3.0327 

6.6927 

-7.0238 

-1.5462 

2.60 

.34481 

-6.0000 

-1.9352 

5.0296 

-6.2593 

-1.3185 

2.86 

.17769 

-3.2308 

-.99615 

2.6385 

-3.3385 

-.7000 

3.12 

.01269 

-.23077 

-.071154 

.18846 

-.23846 

-.05000 


*Angle refers to a F , /3 p , or a 
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Angle* 

breakpoint 

C X 2 ( a F) 

C Xi(%) 

c x w ( Q! w) 

Angle * 
breakpoint 

C X 2 ( a F ) 

O 

X 

t— 1 

'tz> 

c x w («w) 

-3.160 

1.0841 

0.09968 

0.32984 

0.000 

-1.0400 

0.00000 

-0.0050000 

-2.890 

1.2007 

.34484 

-.13508 

.049 

-1.0540 

-.023333 

.0066667 

-2.620 

1.3200 

.59000 

-.60000 

.098 

-1.0680 

-.046667 

.018333 

-2.350 

1.3200 

.59000 

- .70000 

.147 

-1.1010 

-.10439 

.044098 

-2.080 

1.0900 

.37654 

-.58846 

.196 

-1.1371 

-.16785 

.072213 

-1.810 

.8200 

.17508 

-.27692 

.245 

-1.1836 

-.26436 

.12246 

-1.540 

-.42778 

-.026556 

.036667 

.294 

-1.2479 

-.41779 

.21082 

-1.270 

-.8188 

-.23112 

.36240 

.343 

-1.3121 

-.57121 

.29918 

-1.000 

-1.0809 

-.44234 

.61887 

.392 

-1.3180 

-.59000 

.34846 

-.784 

-1.3000 

-.59000 

.69911 

.441 

-1.3158 

-.59000 

.39229 

-.735 

-1.3023 

-.59000 

.65528 

.490 

-1.3135 

-.59000 

.43612 

-.686 

-1.3045 

-.59000 

.61144 

.539 

-1.3113 

-.59000 

.47995 

-.637 

-1.3068 

-.59000 

.56761 

.588 

-1.3090 

-.59000 

.52378 

-.588 

-1.3090 

-.59000 

.52378 

.637 

-1.3068 

-.59000 

.56761 

-.539 

-1.3113 

-.59000 

.47995 

.686 

-1.3045 

-.59000 

.61144 

-.490 

-1.3135 

-.59000 

.43612 

.735 

-1.3023 

-.59000 

.65528 

-.441 

-1.3158 

-.59000 

.39229 

.784 

-1.3000 

-.59000 

.69911 

-.392 

-1.3180 

-.59000 

.34846 

1.000 

j -1.0809 

-.44234 

.61887 

-.343 

-1.3121 

-.57121 

.29918 

1.270 

-.8188 

-.23112 

.36240 

-.294 

-1.2479 

-.41779 

.21082 

1.540 

-.42778 

-.026556 

.036667 

-.245 

-1.1836 

-.26436 

.12246 

1.810 

.8200 

.17508 

-.27692 

-.196 

-1.1371 

-.16785 

.072213 

2.080 

1.0900 

.37654 

-.58846 

-.147 

-1.1010 

-.10439 

.044098 

2.350 

1.3200 

.59000 

-.70000 

-.098 

-1.0680 

-.046667 

.018333 

2.620 

1.3200 

.59000 

-.60000 

-.049 

-1.0540 

-.023333 

.0066667 

2.890 

1.2007 

.34484 

-.13508 





3.160 

1.0841 

.09968 

.32984 


*Angle refers to dp, j3p, or a yj. 
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^VS 

breakpoint 

C Y VS (0 V S) 

'Vs 

breakpoint 

c y vs (0vs) 

-3.06 

0.14240 

0.10 

-0.24333 

-2.85 

.51620 

.20 

-.48667 

-2.64 

.89000 

.30 

-.73000 

-2.43 

.93000 

.40 

-.79250 

-2.22 

.93000 

CJ1 

o 

-.85500 

-1.60 

.70000 

.60 

-.91750 

-1.50 

.73239 

.70 

-.93000 

-1.40 

.76479 

.80 

-.93000 

-1.30 

.79718 

.90 

-.92676 

-1.20 

.82958 

1.00 

-.89437 

-1.10 

.86197 

1.10 

-.86197 

-1.00 

.89437 

1.20 

-.82958 

-.90 

.92676 

1.30 

-.79718 

-.80 

.93000 

1.40 

-.76479 

-.70 

.93000 

1.50 

-.73239 

-.60 

.91750 

1.60 

-.70000 

-.50 

.85500 

2.22 

-.93000 

-.40 

.79250 

2.43 

-.93000 

-.30 

.73000 

2.64 

-.89000 

-.20 

.48667 

2.85 

-.51620 

-.10 

.24333 

3.06 

-.14240 

.00 

.00000 
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“HS 

breakpoint 

- 3.12 
- 3.03 
- 2.94 
- 2.85 
- 2.76 
- 2.67 
- 2.58 
- 2.49 
- 2.40 
- 2.31 
- 2.22 
- 2.13 
- 2.04 
- 1.95 
- 1.86 
- 1.77 
- 1.68 
- 1.59 
- 1.50 
- 1.41 
- 1.32 
- 1.23 
- 1.14 
- 1.05 
-.96 
-.87 
-.78 
-.69 
-.60 
-.51 
-.471 
.00 


C Z hs ( Q! HS) 

0.052766 

.29021 

.52766 

.76511 

1.0026 

1.2400 

1.2741 

1.3081 

1.3422 

1.3762 

1.3630 

1.3439 

1.3247 

1.3056 

1.2865 

1.2674 

1.2482 

1.2291 

1.2100 

1.2340 

1.2580 

1.2820 

1.3060 

1.3300 

1.3540 

1.3728 

1.3812 

1.3896 

1.3980 

1.4064 

1.4100 

.0000 


a HS 

breakpoint 

0.471 

.51 

.60 

.69 

.78 

.87 

.96 

1.05 

1.14 

1.23 

1.32 

1.41 

1.50 

1.59 

1.68 

1.77 

1.86 

1.95 

2.04 

2.13 

2.22 

2.31 

2.40 

2.49 

2.58 

2.67 

2.76 

2.85 

2.94 

3.03 

3.12 


c z hs (<*hs) 

- 1.4100 

- 1.4064 

- 1.3980 

- 1.3896 

- 1.3812 

- 1.3728 

- 1.3540 

- 1.3300 

- 1.3060 

- 1.2820 

- 1.2580 

- 1.2340 

- 1.2100 

- 1.2291 

- 1.2482 

- 1.2674 

- 1.2865 

- 1.3056 

- 1.3247 

- 1.3439 

- 1.3630 

- 1.3762 

- 1.3422 

- 1.3081 

- 1.2741 

- 1.2400 

- 1.0026 

-.76511 

-.52766 

-.29021 

-.052766 
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TRIM CIRCUIT ALGORITHM 

By Gary A. McDaniel 
Electronic Associates, Inc. 

A generalized trim algorithm taken from reference 6 and described in this appendix 
has been used successfully in several RTS programs at the Langley Research Center. 

This method, a generalized secant method, is applicable since it requires no first partial 
derivatives in making the linear approximation to the trim equations. This is a desirable 
feature inasmuch as many RTS programs contain function data rather than analytical 
expressions representing the aerodynamic coefficients and derivatives. 

The generalized secant algorithm is an extension of regula falsi to multivariable 
functions. By letting 0 be an m -dimensional zero vector, the problem can be stated: 

Given an m-dimensional vector x = fxp . . ., x m j of independent elements (called the 
variable vector) and an m-dimensional vector F(x) = (Fj(x), . . ., F m (x)j whose ele- 
ments depend on x (called the function vector), find an m-dimensional solution vector 
x D such that F(x 0 ) = 0. The vector x Q is computed by the following steps. Pick a 
nominal value of x 1 . Obtain m additional variable vectors by successively varying 
only one component of xl. There is now a set of m + 1 variable vectors, that is, 
xl, x2, . . x m+1 . With each member of this set, associate a function vector 

F^, F^, . . F m+ 1. From this set of function vectors, compute an m + 1 dimen- 
sional sum vector S = (f^ • F^, F^ . F^, . . F m+ 1 - F m+ ^). Use this S vector to 

find 

m+1 m+1 

gmax _ max and gmin _ min g_^ 

i=l i=l 

With S max and S m * n , there is associated in an obvious way x max , F max , x m * n , and 
pmin pick a value for K and solve the following linear system for an m+1 — dimensional 
vector q = /qj, q£> • • •> w hich is called the q-vector: 


m+1 

I “j - 1 (D1) 

3=1 

m +l 

J q^j = (1 - KjF™ 111 (D2) 

j=l 
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where 0 < K = 1. Then, compute a new variable vector from 
m+1 

*r - 1 ^ 

j=i 

Associated with this new variable vector will be a new function vector F new and a new 
sum vector s new = F new . If S new > S max , halve K and resolve equation (D3) for 
x new Continue this procedure until S new = S max . When this occurs, replace F max 
by F new in equation (D2). Compute a new S, S max , and S min . Replace F min on 
the right-hand side of equation (D2) by the new F min . Replace x max on the right-hand 
side of equation (D3) by x new and resolve equation (D2) until F new = 0 to within some 
tolerance. Therefore, K is doubled up to a maximum of 1 each time S new = S max 
within three trials. 
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TABLE I.- COMPUTER PROGRAM STATIC CHECK 
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TABLE II.- COCKPIT INSTRUMENT CHECK 


DAC 

. i 

Instrument 

Element 

Variable 

Voltage 

Name 

Reading 

25 

cos(4>) 

-86.54 

^Attitude direction 

+30° 

26 

sin(4>) 

-50.04 

J indicator - roll 


27 

cos(0) 

-70.71 

^Attitude direction 

+45° 

28 

sin(0) 

-70.71 

j indicator - pitch 


29 

30 

COS('I') 

sin('I') 

-86.58 

-50.04 

^Heading indicator 

+30° 

31 

32 

cos(h) 

sin(h) 

-95.11 

+30.90 

^Altimeter 

50 ft 

33 

34 

cos(h) 

sin(h) 

-46.96 

+88.29 

j>Rate of climb indicator 

1500 ft/m in 

35 

36 

cos(r) 

sin(r) 

-86.58 

+50.00 

J Rate of turn indicator 

2 needle widths right 

37 

38 

cos(-a Y /a z ) 

sin(-a Y /a z ) 

+41.62 

-90.93 

j> Bank indicator 

1 ball width right 

39 

Cl 

+64.80 

Main rotor rpm indicator 

323 rpm 

40 

a spd 

+72.22 

1 

Airspeed indicator 

130 knots 
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TABLE ID.- CONTROL DEFLECTION CHECK 


Cockpit control 

ADC 

Name 

Range 

Element 

Voltage 

Variable 

Collective stick 


1 

0 to +80 

X AOS 

Lateral cyclic stick 


2 

-60 to +60 

Y CS 

Longitudinal cyclic stick 

Full forward to 
full back 

3 

-60 to +60 

x cs 

Tail rotor pedals 


4 

-60 to +60 

X TR 

Lateral cyclic trim dial 


5 

-100 to +100 


Longitudinal cyclic trim dial 

Full left to 
full right 

6 

-100 to +100 


Tail rotor pedal trim dial 

Full left to 
full right 

7 

-100 to +100 
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(a) Typical program control station. 
Figure 1.- Operational control features. 


CONTROL STATION 
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(b) Closeup of control panel on the program control console. 



t, sec 


(a) Recorder 1. 

Figure 2.- Dynamic check at 148.16 km/hr (80 knots) with A^g doublet input. 
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(c) Recorder 3. 
Figure 2.- Concluded. 
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(b) Recorder 2. 
Figure 3.- Continued. 
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(b) Recorder 2. 
Figure 4.- Continued. 
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(c) Recorder 3. 
Figure 5.- Concluded. 
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Reference 1 simulation program 


O Unpublished reference 1 simulation program 
□ Langley simulation program 




60 r~ 



full stick 
deflection 


40 L 


I I 1 1 J 

0 15 30 45 60 

Forward velocity, m/sec 

1 I 1 1 1 

0 50 100 150 200 

Forward velocity, ft/sec 
Figure 6.- Comparison of static trim stability. 
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— Reference 1 simulation program 
O Unpublished reference 1 simulation program 
O Langley simulation program 
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Figure 6.- Concluded. 
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